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ON THE ORIGIN OF SOLAR FLARES 
T. Gold and F. Hoyle 
(Received 1959 May 27) 
Summary 


Two major problems in the interpretation of solar flares are discussed: 
firstly, by what means high concentrations of energy can be stored in the 
chromosphere; and secondly, what process can account for the suddenness 
with which this energy is released during a flare? It is shown that the energy 
storage can be accounted for only by a particular class of magnetic fields 
whose lines of force have the general shape of twisted loops protruding above 
the photosphere. In active regions with strong magnetic fields it is likely that 
chromospheric motions are entirely dictated by the flow of the much more 
massive photospheric material. Forces arising from pressure, weight or 
acceleration of chromospheric gas are most of the time small, so that the 
motions there are closely those that retain at each moment a force-free field 
compatible with the magnetic boundary conditions at the photosphere and, as 
the chromosphere is a good conductor, with the past history of the motion. A 
sudden release of energy in such a system can result only in a case where the 
magnetic forces act so as to drive the system away from the force-free con- 
figuration, and the energy associated with chromospheric currents can then be 
dissipated into motion or heat. This will occur when twisted magnetic 
loops of opposite sense and opposite twist meet. Such loops attract each 
other, and the annihilation of the longitudinal component of the field where 
they meet leads to a sudden constriction of the current and through this to a 
dissipation of the energy associated with that current. 





1. Introduction.—Of the phenomena that happen in the solar atmosphere, 
flares are the ones connected with the greatest amounts of energy and with the 
greatest range of associated observable effects. An understanding of the basic 

, physical processes underlying the flare phenomenon is therefore required for 
the interpretation of the great range of solar observations. 

Attempts have been made to explain the phenomenon of solar flares in terms 
of electromagnetic processes. Early work by Giovanelli, by Hoyle, and by 
Dungey recognized the importance of neutral points and neutral lines of the 
magnetic field in this connection. Recent observational and theoretical work by 
Severny has gone far towards demonstrating the correctness of the magnetic 
interpretation in general and of the importance of the neutral points and lines in 
particular (1, 2, 3, 4). 

Quite apart from the cogency of the direct observations reported by Severny, 
there is a compelling reason why a magnetic interpretation of flares is necessary. 
No other way has been found that can explain the concentration of energy in 
the flare region, or for its sudden transport there during the flare phenomenon 
itself. 

The optical emission during an intense flare can be estimated to amount to 
more than 5 x 10’ ergs/cm*/sec. (The chief uncertainty here is the degree of 
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concentration of the source into filaments, which may be seen with a resolution 
that is possibly inadequate.) Thus for a flare with a duration of approximately 
1000 seconds the total optical emission is of the order of 5 x 10 ergs/em*. For 
an intense flare occupying an area of, say, 1 per cent of the solar disk, the total 
integrated emission amounts to 10*' ergs. 

The total energy stored in the responsible regions of the solar atmosphere 
must thus amount to at least 5 x 10 ergs for each cm? of the solar surface. If, 
further, we were to assume that flare emission takes place uniformly over a total 
interval of height as great as the entire height above the photosphere at which 
the phenomenon is known to occur, a height of the order of 10cm, the required 
energy storage per cm* would amount to approximately 5 x 107 ergs. Although 
such an estimate is presumably not accurate in a strict sense, it suffices to show 
that the energy density in flares is much greater than the normal thermal energy 
density, which amounts to only approximately 10 ergs percm®* for material of 
density 10’*atomspercm® at a temperature of 7000 °K. 

The kinetic energy content of the gas in the chromosphere before a flare is 
also quite inadequate. This is most easily seen from the fact that Ha lines of 
very great width are seen in intense flares which, when interpreted as arising 
from mass motions of material, correspond to speeds of 80-250km/sec. This 
is about 10 times greater than the speeds that are normally found in such regions 
of the chromosphere. The dynamical energy per unit mass released in flares 
therefore appears to be about 10? times the normal value. It may also be noted 
that the energy density, as estimated from optical considerations, in material 
with particle density of the order of 10'* atoms per cm’, is similar to the kinetic 
energy density of such material moving at 100km/sec, namely approximately 
10° ergs per cm’. 

We are therefore forced to the conclusion that energy must be stored in the 
solar atmosphere in some way not involving thermal or mass motions, and that 
this energy storage must be built up before the flare can take place. The buildup 
must occur over prolonged periods of time, since any rapid method of transferring 
and concentrating energy would itself lead to visible effects. 

If the energy is stored in the form of a magnetic field, then H?/8z must be of 
the order of 10% ergs percm*, and the magnetic field intensity must hence be a 
few hundred gauss. 

In addition to solving the energy storage problem, a theory of flares must 
account for the extreme suddenness of the event. ‘There are no other processes 
on a comparable scale occurring on the Sun that take place so quickly, and some 
form of catastrophic instability must therefore be involved. The requirements 
of the theory can therefore be stated quite definitely. Magnetic field configura- 
tions must be found that are capable of storing energy densities hundreds of times 
greater than occur in any other form, and that are stable most of the time. A situa- 
tion that occurs only a small fraction of the time must be able to lead to instability 
in which this energy can rapidly be dissipated into heat and mass motion. 

There are two possibilities concerning the suddenness of the phenomenon. 
Either it is in the nature of a breakdown, occurring when the stored energy 
exceeds a certain value, so that the process of charging up the region with more 
energy is itself the triggering mechanism. Or, alternatively, the onset of the 
instability is triggered by some circumstance which occurs only very rarely, so 
that the energy-storing process can usually proceed without hindrance. In this 
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case it would be possible for the storage process to be a very common one, and a 
fertile ground for flares might hence exist over most active solar regions—but 
with the initiation of the discharge a sufficiently improbable process to limit 
flares to the observed rate of occurrence. Seeing that a great variety of com- 
plicated magnetohydrodynamical motions are known to occur over active regions, 
as observed by Ha cinematographic techniques, it is clear that in this second 
case a very special condition must be postulated as the requirement for the 
initiation of the discharge. 

2. Force-free fields.—Strong fields that are to be present for prolonged periods 
in the chromosphere must be largely force-free, otherwise rapid motions with 
velocities of the order of H(47p)~"? where p is the mass density, would immediately 
be developed. This requirement would of course be met by fields that are due 
entirely to currents flowing in the massive photosphere below, but such fields 
are not of interest to us here, for there is no way in which this field energy can be 
released except in the regions where the currents responsible for it are flowing. 
We will thus have to be concerned here with force-free fields in which at least 
some of the currents that are flowing are situated in the chromosphere and there- 
fore in the force-free part. Forces arising from the term j x H must be small 
in the chromosphere, not because j is absent, but because it does not have a 
significant component normal to H. 

The ultimate energy source for any large manifestations of free energy in the 
solar atmosphere must no doubt be sought in the hydrodynamics of the con- 
vective motion in and below the photosphere. We therefore wish to find a 
component of this turbulent motion, and a certain initial field configuration, 
which can result in a gradual buildup of chromospheric currents along magnetic 
lines of force. Whilst the observational evidence is not very definite in defining 
the length of time for which such a buildup may be taking place, it seems certain 
that a time much longer than the observed flare process is involved. No sudden 
processes are seen in the photosphere, and no mechanism could be suggested that 
would supply the energy from below the photosphere at the rate at which it is 
observed to be released. Also it is certain that flares tend to occur along arcs 
that are closely related to chromospheric patterns, and again for this reason a 
sudden photospheric supply, concurrent with the flare radiation, seems ruled 
out. 

We are led to think then of a flare process in two separate stages. Firstly, a 
buildup of energy storage in the chromosphere by means of currents directed 
along magnetic lines of force. These currents must be generated through the’ 
interaction of the magnetic field with a component of photospheric turbulent 
motion. The second process is an explosive release of the stored energy into 
heat, bulk motion, emission line and possibly also continuum light, as well as high 
energy particles. 

3. Dissipative decay of currents.—Before proceeding with the discussion of the 
flare mechanism, it is necessary to establish the circumstances in which dissipative 
decay of electric currents may be significant. It can be seen that force-free 
chromospheric currents, running parallel to the magnetic field, will generally 
decay in time scales that are very long compared with those arising from other 
considerations. In the force-free case the time scale for decay is of the order 
aL?/c?, where o is the conductivity parallel to H, L is the characteristic dimension 
of the system under consideration, and c is the velocity of light. The value of 
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o at temperatures in the range from 5 000 °K to 10000 °K can be taken as 10"* sec 
independent of the density (5). Thus for L ~10*cm, the time scale for decay 
is 10° sec, which is long compared with the times that are normally of importance 
in chromospheric processes. 

The situation is more complicated for a component of j normal to H, 
especially in the lower chromosphere where the conductivity is much decreased 
by the presence of neutral hydrogen, as was first pointed out by Piddington (6). 
An expression for the conductivity in directions normal to H has been given by 
Cowling (7) in a form that is convenient for the present discussion. When the 
great majority of the atoms are neutral, as is the case before the onset of a flare, 
the conductivity is less in directions normal to H than it is in the direction parallel 
to H by the factor, 


2p72 
+ — [My(v+v,)+mvv,J> (1) 


where the symbols have the following meanings: 
e, m are the electron charge and mass respectively, 
H is the magnetic intensity, 
c is the velocity of light, 
M is the proton mass, 
v is the collision frequency between electrons and protons, 
v, is the collision frequency between electrons and neutral hydrogen atoms, 
v; is the collision frequency between protons and neutral hydrogen atoms. 


Formula (1) assumes that the ionization of hydrogen is more important than the 
ionization of other atoms that are present in the gas, such as the metals. If another 
atom, magnesium for example, were the most abundant positive ion, then it would 
be approximately correct to replace protons by magnesium ions in the above 
definitions. 

Since the cross-sections for collisions of electrons and protons with neutral 
atoms are not greatly different, we have v,~(M/m)"»v, (the factor M/m arises 
simply from the greater thermal speeds of the electrons), and hence 


mv, = (Mm)"?v,< My,. 
It follows therefore that (1) can be written to sufficient accuracy as 
2 p72 —1 
re H ad (2) 
mMc Vv; 
The quantities v, v,,v; can be determined from the following equations: 


v Son u,, 
Ve OMe, (3) 
vj = OMu;, 
where 
o is the collision cross-section between electrons and protons, 
o, is the collision cross-section between electrons and neutral atoms, 
go; is the collision cross-section between protons and neutral atoms, 
n, is the ion-density, 
n is the density of neutral atoms, 
u, is the mean thermal speed of the electrons, 
u; is the mean thermal speed of the protons. 
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It is a satisfactory approximation, both in the lower chromosphere and at the 
photosphere, to put 


o=3x10-"cm-, 

o,=0;= 10-* cm, 

u,=5 x 10’ cmsec™, 

u; = 10° cm sec", (4) 


from which », v,, v; can immediately be estimated as soon as n, nm; are specified. 

At the photosphere we may put a~10'’ atomscm~-, a;~10ionscm~*. 
These values, together with (4), enable (2) to be expressed as 1+ ~ 3 x 10-7H?. 
This factor of reduction for the conductivity normal to H is evidently of order 
unity, even for field intensities as high as 10° gauss, and hence is of little importance 
at the photosphere. 

The situation is very different in the regions of the chromosphere that are 
relevant to the flare problem, however. Here the ratio n,/n is probably of order 
1ro~* under normal conditions (m = 10* atoms cm~*, n,~ 10" atoms cin-*), and 
must become considerably larger during flare conditions. This implies that v 
always much exceeds v,, so that v, can be neglected in (2), which accordingly 
reduces to 

2f72 
1+ ee, j (5) 


mMc* vv; 
With x defined by 
n;= xn 
(5) becomes 
10242 
I+ — . (6) 

With n=10%atomscm~*, x~10-*, H?~10*%ergscm-*, the value of (6) is 
~1o0*. Currents normal to H decay faster than currents parallel to H by this 
factor. Remembering that we have already shown that the characteristic time 
scale for the decay of currents parallel to H is of order 10° sec (for a system with 
linear dimensions ~ 10*cm), it follows that the time scale for the decay of the 
current normal to H in such a system is ~ 10‘ sec. 

The time scale for decay becomes still shorter for smaller values of n, At 
n= 10 atomscm~*, x~10~!, H*= 10‘ergscm~, the time scale ~ 10*sec; and 
at n~10"atomscm~*, x~1, H*z10tergscm~*, the time scale is ~ 10? sec. 
These values are evidently of an appropriate order to be of importance in the 
flare problem. 

While a density of 10'* atoms/cm* may be a good estimate for the emitting 
material of a flare, such material has probably experienced some degree of com- 
pression (4). Hence we should think of the initial particle density as being 
somewhat lower during the phase that precedes the flare. A value of n= 10 
atoms/cm? might well be more appropriate in the initial situation. 

4. The energy storage mechanism.—Observations of flares seen in the light of 
Ha against the solar disk show them frequently clearly aligned along a pre-existing 
filamentary pattern. Such filaments, when seen against the solar limb, are thought 
to be the commonly occurring arch structures. (Although these markings may 
sometimes represent high-level coronal prominences—which are known not to 
be directly related to flares—most of the patterns visible in Ha are undoubtedly 
related to arch structures in the chromosphere and lower corona.) Such filaments 
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should most probably be thought of as bundles of lines of force of a magnetic 
field that emerge from sub-photospheric regions at one root and re-enter at 
another root. 

Some such firm anchoring in the massive photospheric material seems to be 
essential if field strengths of the order of 100 gauss are to be maintained in the 
chromosphere for long periods of time. This anchoring then makes possible a 
class of field configurations that are force-free above the photosphere, but not 
below it. Despite their high energy density, such fields will not cause an explo- 
sion above the photosphere, while below the photosphere the energy density 
will in no case exceed that of the local material. 

In fact the energy of mass motion of the highly convective sub-photospheric 
material is of the order of 10° ergs/cm*, which is greater than that required in 
chromospheric filaments by a factor of the order of 100, if these are to be the 
seat of the flare phenomenon. No energy difficulty will therefore be encountered 
if a process can be found whereby magnetic energy can be fed into the chromo- 
sphere as a consequence of events that take place at and below the photosphere. 
If this energy is to be dissipated in the chromosphere, it must be associated with 
currents at that level. Since the strong fields must be force-free in the chromo- 
sphere, these currents must flow along the lines of force. We may therefore ask : 
what are the photospheric motions that generate electric currents along the lines 
of force, in those cases where the lines of force stretch up into the chromosphere 
in the form of an arch? 

In general an electric field will be established between the two ends of a magnetic 
arch by horizontal motions in the photospheric material. The integral of E.ds 
along a line of force joining the two bases of an arch evidently must not vanish. 


This condition will in general be satisfied by any non-uniform motion in the 
photosphere. 


cop 





Fic. 1.—Bundle of lines of force protruding from the photosphere into the chromosphere, with 
a twisting motion at the base. 

It is perhaps easiest to discuss this question in terms of the concept of the 
convection of lines of force with the flow (cf. Figs. 1 and 2). The generation of 
a current along the lines of force of an arch is equivalent to generating a field with 
the structure of a twisted arch. The non-uniform motion in the photosphere 
required to send the current along the arch is therefore the motion that will twist 
up the initial arch. A horizontal rotation in the photosphere around each base 
in the same clock-sense is therefore the basic type of motion that generates force- 
free, but not current-free, fields at a higher level. 
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It may be thought that the motion may not be exactly the one required for a 
force-free field, but of course it should be realized that, for fields of the strengths 
that we are concerned with, the chromospheric material will always immediately 


VxH VxH 





Fic. 2.— Bundle of lines of force after a twisting motion from the base has generated a jrelical 
shape. This is equivalent to considering the establishment of a current, j, along the lines of force, 
due to the V x H terms arising at the base. 


move in a direction in which any force is applied. As photospheric motions 
cause an initially untwisted arch to become twisted, the arch moves at all stages 
into a force-free configuration. 

This in turn means that we are concerned with currents flowing in the chromo- 
sphere along the magnetic lines of force. Hence the currents do not suffer the 
greatly enhanced rate of decay that would occur for chromospheric currents 
flowing in other directions. For currents along the lines of force the decay time 
constants are long enough to make a discussion in terms of the convection of lines 
of force quite applicable for the speeds and dimensions that are in fact observed. 

A torque acts across planes normal to the axis of such a twisted filament. It is 
important to realize that this is the case not only at the photosphere, but also over 


the force-free chromospheric portions of the arch { the force-free condition does 


not require the individual components of the stress tensor p,, to vanish, only 


its divergence 
i) 
3 Pik) 
Pa OX, 


For the purpose of estimating this torque we anticipate a result obtained in the 
next section; namely that the number of turns of the lines of force about the 
filament axis, N say, is related to the length of the filament, / say, and to the 
characteristic dimension d of its cross-sections, by 


Nw~lId. (7) 


This requires the stress acting per unit area on a normal plane to have an 
appreciable component about the axis of the filament. The couple about the 
axis accordingly has an average value, of order H*d/87 per unit area, and a value 
of order H*d*®/8z for the whole cross-section (this estimate will be confirmed by 
a more precise calculation in the next section). 

The same torque acts across all normal planes of the force-free portion of the 
filament. ‘The torque need not be exactly the same, however, at the non-force- 
free roots of the filament, because angular momentum about the axis may be in 
the process of being communicated to the material in the lower parts of the fila- 
ment. But since the twisting is impressed on the filament at its roots, the torque 
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at the roots can differ only in the sense of being greater than the torque in the 
force-free part of the filament, i.e. greater than H*d*/8r. Thus in using the 
latter formula we shall not be underestimating the torque at the roots. 

The work done in N relative rotations at the roots is therefore not less in 
order of magnitude than 


(27N).(H?d*/87). 


Now we have seen in Section 1 that H?/87~10*ergscm~*, so that with 
d~2x10*cm for the cross-sectional dimension of the filament, the work done 
is at least of order 5 x 10%Nergs. ‘The relation (7) requires N to lie between 
10 and 10? for a typical ratio of filament length to diameter, so that the total energy 
storage is not less than 10” ergs, which is of the correct order for a normal flare. 

These considerations therefore demonstrate the effectiveness of twisting as a 
mode of energy storage. A precise assessment of numerical values should be 
deferred, however, until the end of the next section, when the calculation will be 
repeated with greater accuracy. 

5. The structure of a force-free twisted filament.—We shall suppose the curva- 
ture of the axis of the filament to be sufficiently small for a limited segment of 
it to be considered straight. We also make the following assumptions : 


(i) all points of a line of force are at the same distance from the axis, 

(ii) all lines of force have the same number of turns about the axis, 

(iii) the field does not depend on the azimuthal coordinate about the axis, nor 
on the distance coordinate parallel to the axis, at any rate in first 
approximation. 

Assumption (ii) implies that the angular velocity of twisting is uniform over the 
roots of the filament. ‘The most questionable of these simplifications is that of 
axial symmetry in (iii), since in an actual filament the horizontal width is likely 
to exceed the vertical depth. 

We use symbols with the following meanings: 

r, 9, = are cylindrical coordinates with respect to the axis, 

¢ is the angle that a line of force makes with the planes normal to the axis, 

v is the number of turns per unit length of a line of force. 

Because of (i), the radial component of the field vanishes. Because of (iii), 
¢,v and the components H,, H, of the field are all independent of @ and z. 
Because of (ii), and the definitions of v and 4, it follows that 

cot d= 277r. (8) 
The force-free condition 
H x curlH =o (9) 


yields @ and z components that are identically satisfied, while the r-component 
of (9) gives 


dlogH ___cos*¢ 
oo ee 


H=|H|. (10) 


It is most important in relation to the later discussion of the stability of the 
flare phenomenon to consider whether a filament satisfying the above conditions 
can be regarded as an isolated structure. To this end, we first consider the 
possibility of satisfying (9) for r<a, together with 


H=0, r>a. (11) 
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It will be useful to write F for the flux through a normal section of the filament, 


Faan["Hs de. - (12) 
0 


The importance of F lies in the fact that it is independent of the twisting of the 
filament. It also turns out that F can take the place of the constant that appears 
in the integration of (10), as will soon be seen. 

When ¢ is eliminated between (8) and (10), the resulting equation yields 


H=A/V1+4rP, (13) 
where A is an integration constant. Substituting H.=Hsin¢g in (12), with 
sing determined from (8), also gives 

z qF 
~ ma*log(1+4q)’ 


The results (13), (14), together with H, = Hcos¢, H,= Hsin ¢, lead immediately 
to 


q= 47° v*a?. (14) 


ce Fq I 

~ wa* log (1 +q) (1 + gr?/a?)!?” 
ne Fq 2nur 

 ma® log (1 +9) 1+ 9r2/a?’ 

sa Fq I 

* wa® log (1 +9) 1+ 9r?/a?" 











We also note for later reference that the total magnetic energy per unit length 
is given by 
F?y? 


log(1+q) Si 


[re 


The energy per unit length of the axial component is given by 


a se F?,2 44 q I sh 
b [Her d= : SORTER (17) 


while the energy per unit length of the longitudinal component follows immediately 
on subtracting (17) from (16), 


a ei F*y?g 
iter cry ice) _ 

Manifestly the field does not fall to zero at r= a, as had been hoped, in order 
that (11) might be satisfied. For r— oo, H declines to zero as 1/r, H, as 1/r’, 
H, as 1/r. On the axis H,=o. Hence we have an essentially longitudinal field 
near the axis giving place to an axial field at large r. 

The tendency of the axial field to pinch is resisted by the longitudinal field, 
and this situation appears stable against a purely radial perturbation. Suppose 
a,r to be increased by some scale factor 7. Then while H, decreases as r~*, H, 
only decreases as t~'. The pinch thereby becomes the stronger, so that the 
filament is compressed towards the axis. An opposite situation obviously holds 
for a decrease of a,r. 

The above considerations show that in order to confine the force-free field 
within a cylinder of radius a, we must surround the filament by a shielding 
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current. Such a shielding current can be thought of as flowing in a boundary 
sheath outside the cylinder of radius a. Plainly this boundary region cannot 
be force-free in the electromagnetic sense. Indeed, while the magnetic stress 
on the outer surface is zero, there is a magnetic pressure 
2 
a ee (19) 
7 (1+ q) [log (1+ )] 
on its inner surface. It appears therefore that the filament cannot be an isolated 
structure unless the electromagnetic pressure (19) is compensated by an excess 
material pressure on the outer surface of the boundary zone. The next step is to 
examine the implications of this requirement. 
By suitable choices of v, F we can arrange that 





x (H*)_.= = (H?),_., = 10% ergscm™~%, (20) 


We then have a situation giving a storage of magnetic energy near the axis of the 
filament that is adequate for a flare but which falls at r =a to a value comparable 
with the thermal energy density of surrounding chromospheric material 
(cf. Section 1 for numerical values). In these circumstances the filament can be 
confined by the shielding current of the boundary zone provided that the material 
pressure is appreciably less within the filament than it is outside. 

From (15), 

(H*),9= (1+ q)(H*),-4» (21) 
so that (20) implies g~ 10. 

From the definition of gin (15) we then have va~1. If Nis the total number 

of turns over the whole length / of the filament v~ N/I, from which 

Nella. (22) 
Thus a is the characteristic cross-sectional dimension introduced in the previous 
section and used in (7). 

With the help of (15) it is easy to work out the total torque across a normal 
plane more accurately than was done in the previous section. Ignoring a small 
contribution from the boundary zone, the torque is given by 

; a H.H,y dr = a*(H?),_ [es +q) _ | (23) 
0 Srv q 1+q 


With va = 1, q= 10?, the right hand side of (23) can be written to sufficient accuracy 
as 





_! 3578 
ae" (H?*),..9 log q. (24) 


This result is similar in form to the expression H?d*/8z, (d= a), used in the pre- 
vious section. It differs in detail, however, in that the torque is reduced below 
the previous rough formula by the factor g~ log q. 

To end the present section, a more accurate assessment of the energy storage 
in a twisted filament will now be given. From (14) and (15) it follows that 


(H),_.,~ 2vF jalog q, (25) 
the approximation g>1+ q being used. Eliminating vF and (H),_, between 
(16), (20), and (25) we thus have 1oma*logg for the energy storage per unit 


length of the filament. For a filament of length / the total energy storage is 
therefore 107a*/ log q. 
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This expression can be put in a convenient form by noting that a*/~ Ah, 
where A is the area of the filament projected against the solar disk and A represents 
height above the photosphere. Writing A=7R*f, where R is the solar radius, 
our result becomes ~ 107*R*hflog q and this is not sensitive either to q or to h, 
the height being always of order 1o°cm. The energy storage is accordingly 
dependent essentially on f, the fraction of the solar disk covered by the filament. 
Inserting h~10%cm, g~10*, R~6-9x 10cm, gives ~2-10*fergs, so that 
for f~10-* the energy storage ~10”ergs—as obtained in our approximate 
considerations. 

6. The boundary zone.—In place of (11) we now write 


H=0, r>b>a. (11') 
We still regard the field for r <a as being given by (15). 


The zone between the cylinders of radii a, b carries the non-force-free shielding 

current. The nature of this current is easily understood. At . 
r=a, H,=2nvaH,, 

so that for va~ 1, the field is mainly axial. ‘Thus the shielding current must be 
mainly parallel to the axis of the filament, and if F>o must be directed in the 
negative z-direction. ‘The total current required to flow in the boundary can 
readily be shown to be ~ $ac(H),_.,. 

If the current density j is assumed for simplicity to be uniform over a normal 
cross-section then 


_ pep ac(H),~4 
=|il= 2n(b? —a*)’ 


(26) 


and the total heat produced is unit length of the filament is 


2a*c* 
= o(b?—a*) * 82 = (A Yrmar 


where a is the effective ae In the present case the current is nearly 
perpendicular to H, so that (6) of Section 3 must be included in estimating the 
conductivity to be used in (27). In fact 


2~ 15-32 7 * 
a/c?= 10 H),_, (28) 
The next step is to estimate 6. ‘This can be done from the usual skin-depth 
formula for an alternating field of angular frequency w. This gives a depth 
c/(wo)"?, i.e. 


n(b*—a?).L = (27) 


citi 
(wo)'?" 


For w we set 7-!, where T is the time scale for setting up the twisted filament. 
. 1/2 
Thus eae (=) (29) 


o 


b-—az 


We now consider numerical values. With 
(H?),_,./827 =10ergcem-*, n=10%atomscm-*, x=107%, 
(28) gives o/c?=3-10-". With 10°sec as a reasonable estimate for T, we have 
b-—a~ 5 10’ cm, from which it follows that the shielding current can be confined 
within a zone of reasonable thickness. ‘The zone becomes thinner at higher n, 





100 T. Gold and F. Hoyle Vol. 120 


but somewhat thicker at smaller n. Thus, lowering the particle density to 
10" atomscm~*, x= 10~', but keeping 7, (H),_, the same, gives 


b—a=1-°5 x 10° cm. 
Since we have considered a~ 10° cm it is clear that the boundary zone is such 
that b—a~a over the relevant chromospheric portions of the filament. 
If as an example we put 6—a=a in (27) and (29), o can be eliminated to give 


2 (H*),_, 
Sioa a Ge) 


for the rate of heat production per unit length of the filament. Multiplying 
(30) by T gives a*(H?) _ ,/127 for the total heat production per unit length taken 
over the whole time interval of the winding process. Using (25) for (H),_, 
then gives 
F2y2 
a 1 
37 (log q)? (31) 


A comparison with (16), which gives the total magnetic energy storage per unit 
length, shows that for q~= 10? the dissipation amounts to only about 3 per cent 
of the stored energy. ‘The fact that energy dissipation is so small is of course 
already clear from the form of (30). This shows that only the energy content 
of a magnetic field of intensity (H),_, is dissipated by ohmic losses, whereas 
the main energy content of the field within the filament comes from regions near 
the axis, where the magnetic intensity is much greater than it is at r=a. 

It will be recalled that it was assumed in Section 4 that ohmic decay does not 
prevent the filament field from becoming twisted. The present result gives 
strong support to this assumption. 

7. The instability of twisted force-free filaments.—Although it is tempting to 
argue that filaments of the sort considered above become dynamically unstable 
(and that the instability constitutes the flare phenomenon), this point of view 
seems to us untenable. We have stressed at the outset that two quite distinct 
problems are involved in the flare phenomenon, the problem of energy storage, 
and the problem of almost catastrophic dissipation. The process of storage is 
comparatively slow and continuous. If such a process were interrupted by a 
catastrophic instability then one must understand clearly why the instability of 
necessity must arise at a sharply defined moment (it must not happen until there 
has been adequate energy storage!). We have not been able to see any plausible 
grounds for the existence of such a sharply defined moment during the winding 
of a single twisted filament. 

The solution to the problem would seem to lie in considering two filaments. 
(cf. Fig. 3) This we can properly do, since the filaments are distinct structures, not 
parts of a connected magnetic complex. Suppose two parallel filaments touch 
each other along a line. ‘The currents flowing in their boundary regions are 
either substantially parallel or anti-parallel. If the axial components of the 
fields within the filaments have the same sense about their respective axes the 
boundary currents are parallel, and the filaments attract each other at their line 
of contact (the boundary zones are not force-free). Moreover the electromagnetic 
forces are of necessity as large as the thermal pressure gradients in surrounding 
material, and can therefore be adequate to keep the filaments pressed together 


for a considerable time interval. 
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Fic. 3.—Two bundles of lines of force, both twisted up. The initial longitudinal field is in 
opposite directions, and the twisting has occurred in the opposite sense. The arrows indicate the 
teoo components of the field. The electric current responsible for the circumferential component is 
in the same direction in each one. ws 

Now formula (29) for the size of the boundary zone is a skin-depth formula 
with the meaning that lines of force can diffuse relative to matter by a distance 
b—ainatime 7. Thusif the two filaments remain pressed together for this length 
of time they must begin to penetrate into each other. If we now add the second 
requirement, that the longitudinal fields within the filaments be directed in 
opposite senses, the respective current systems must continue to attract each 
other however far the interpenetration proceeds. Thus the axial field com- 
ponents are generated by longitudinal currents directed in some sense—which 
are clearly attractive, while the longitudinal field components are generated by 
axial currents in opposite senses about the respective axes—and at a line of contact 
such axial currents are parallel (not antiparallel!) and therefore are again 
attractive. Furthermore, as the axes of the filaments approach each other the 
neighbouring current elements become stronger, so that forces much greater 
than those in the boundary zones can be developed. 

All this can be stated more clearly in the following way. We are considering 
two filaments with the same sense in the axial field component and with opposite 
senses in the longitudinal component; this requires the lines of force to be 
right-handed spirals in one case and left-handed in the other. Interpenetration 
tends to augment the axial component and to annihilate the longitudinal com- 
ponent, leaving conditions suitable for the pinch effect to operate (4). 

The most important point still remains, however. To produce an almost 
catastrophic effect, the axial parts of the filaments, where the fields are most 
intense, must approach one another at an increasing rate. We have seen that 
penetration through the boundary zones requires a time of the same order as the 
winding time 7, which we took ~10*sec. Unless interpenetration is much 
speeded up beyond this there is clearly no possibility of a marked instability 
arising. 

It is here that the dependence on the magnetic intensity of the conductivity 
for a non-force-free current plays a crucial role. Thus the conductivity at mag- 
netic intensity H is given by 

a as 8° 


a ee (32) 
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which agrees with (28) at H=(H),_,. Reference to the skin-depth formula 
(29) shows that the time required for diffusion to take place through a given 
distance decreases as o, and hence as H~?, if n,x are fixed. It follows that the 
axial parts of the filaments interpenetrate each much more rapidly than the 
boundary zones. Moreover the strong part of the field is confined within a 
distance ag~"* of the axis. We therefore require penetration not through a 
distance b—a which ~a, but through a distance less than this by the factor g"?. 
Taking both these factors into consideration the strong axial parts of the filaments 
interpenetrate each other in a time 


I (H? )y =a 

@° (Hag ted 
assuming n,x to stay fixed. With (H?),_,/(H?),..=q"', q= 107, T= 10'sec, 
(33) gives an interpenetration time of ~ 10seconds. 

It must now be emphasized that our present intention is only to discuss events 
leading to instability. The present considerations are not therefore intended 
to apply to the final phases of extreme catastrophic instability when dynamical 
considerations must be introduced. 

The assumptions of nm, x constant were made in the above investigation. If 
strong attractive forces develop between the two filaments the particle density 
n must rise at their interface until the gas pressure becomes sufficient to withstand 
the compressive magnetic forces. It is therefore possible that in a shallow layer 
at the interface m may rise substantially above the value ~ 10! atomcm~* used 
above. Our result is not affected by any such one-dimensional compressive 
effect, however (the compression occurs in a direction normal to the surface of 
separation of the filaments). To see this, we first note that the interpenetration 
time of the filaments is proportional to oD*, where D is the interpenetration 
distance. With o proportional to n?, the time is thus proportional to (nD)?—and 
for a fixed quantity of gas within the filaments nD is independent of the degree 
of compression, this being one-dimensional. Thus the more the magnetic 
forces cause n to increase the shorter the distance through which the lines of 
force need interpenetrate each other—the two effects simply cancel each other. 

The extreme catastrophic instability will, however, be mitigated by an 
increase in the ionization factor x, but even a rise tox ~ 1 would still give an inter- 
penetration time ~ 10° sec. It is perhaps worth writing down the relevant results 
again. The skin depth formula shows that interpenetration takes place through 
a distance aq~" in a time oa*/qc*, and with o/c? given by (32), this can be written 

_g2 77a*x (34) 

qi * 34 

The most favourable case occurs when n is small and H large. With 
n= 10% atoms cm~*, H = 200 gauss, g= 107, a= 2 x 10° cm, x~ I, the formula (34) 
gives a time ~10* sec. These estimates are exactly of the right order to explain 
the observed rapidity of onset of a flare. A value as low as 10" atomscm~* 
for n is not at all implausible, since the pinch effect that follows the annihilation 
of the longitudinal field must increase the density, so that the density after the 
pinch could well be as high as 10 atoms cm. 

One point remains, however. It will be recalled that the quantity x was 
defined as the ratio of the ion density m; to neutral atom density m. Evidently x is 
not limited to the value unity. Indeed x takes a value >1 if almost all the 
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hydrogen becomes ionized. This would greatly increase the time estimate 
given by (34) and would destroy the catastrophic instability of the process of 
interpenetration. The consideration of this point allows a calculation to be made 
of the total flare energy, and this calculation is nearly independent of all the 
foregoing results. 

Plainly x cannot become very large if the longitudinal fields within the 
filaments are not intense enough for annihilation to provide sufficient energy to 
ionize entirely the whole of the filament material. It is easy to work out what 
this limitation implies. If A is the area of the flare projected against the disk, 
and h is the height range of the flare, then the total material content of the flare 
~ AhMn. The energy required to ionize this mass of material ~ 10 AhMnergs. 
For x <1, we thus require the energy of the longitudinal field component,to be 
less than ~10'*AhMn. And since the total field energy exceeds that of the 
longitudinal component by a factor ~ 10, our present condition therefore requires 
the total magnetic energy to be less than ~10'* AhMn. This gives ~ 10*' ergs, 
when n= 10" atomscm~, h~10%cm, and A = 10” cm? (i.e. about 1 per cent of 
the disk). In view of what was said above n should probably be taken less than 
10'3 atoms cm~*, and the energy limit should be reduced correspondingly. - 

If the energy limit is exceeded, the longitudinal component of the filaments 
do not interpenetrate completely at the first instability. Then the material will 
probably radiate without any dynamical instability developing from the axial 
component of the field, since the axial field in this case can pinch on to a 
longitudinal field in both filaments. Eventually, however, the hot hydrogen 
must cool, and recombination of the ionized atoms decreases the conductivity. 
A further interpenetration of the longitudinal field components must then take 
place, until the interpenetration is again reduced by further heating. This will 
go on until the longitudinal field is annihilated in one or the other of the filaments 
—it is of course unlikely that the intensities of the two filaments will be so exactly 
equal that both longitudinal components become effectively annihilated. Since 
the whole process is slow it seems most unlikely that any strong dynamical 
instability will arise from the axial field. Rather does it seem likely that the axial 
components from both filaments will pinch in a stable way on to the residue 
of the longitudinal field. 

It is a possibility that chromospheric flocculi are a manifestation of some 
such process of steady emission. 

To end this section it is worth noting that (16) requires the total energy of 
the whole length / of the filament to be 

: PA 1 
2 log (1+q) ~ 2/logq’ 


where N = vl is the total number of turns of the lines of force. If we now think 
of / as variable with F, N fixed, it follows that the longer the filament the smaller 
the stored energy (the logarithmic term cannot vary appreciably with J). This 
means that random perturbations will always tend to lengthen a filament. In 
our view this is just why filaments are in fact long. An examination of Ha 
spectroheliograms of areas where flares occur shows the filaments to meander 
apparently aimlessly over the areas in question, in just the way that would be 


expected for systems that lengthen under the influence of perturbations, but 
which are otherwise stable. 
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8. General remarks on the dynamical instability.— Although we have not been 
able to consider the dynamics of the instability in a serious quantitative way, 
there are several features that seem worth mentioning. 

Instability seems likely to start at a particular point of the coupled filaments, 
rather than simultaneously along their whole length. This arises from the 
increasing speed of the interpenetration, which demands some unevenness of 
interpenetration. Plainly one point must run somewhat ahead of the others and 
instability will set in first at this special point. So much is clear. It is then 
plausible to suppose that instability spreads both ways along the filaments from 
this special point. This means that instead of waiting for instability to be reached 
separately at each individual place along the filaments we expect that the first 
point of instability will simply trigger off the whole length over which the filaments 
are pressed together. 

Recent work in the laboratory shows that a powerful pinch is followed by a 
second and even by a third pinch. If we assume this to happen in a similar 
fashion in the solar case we can understand at any rate empirically the emission 
of high energy particles by flares, since such particles are in fact produced in the 
laboratory case. 

We have spoken of the interpenetration of two parallel filaments, but it is 
clear that since the filaments have different roots they can only be parallel over 
a limited segment of their lengths. Thus three portions can be distinguished for 
each filament, a portion that is common with the other filament, and two portions, 
one for each root, that are not common with the other filament. 

The destruction of the longitudinal components of the fields of two filaments 
over their common portion allows a considerable change to take place in the 
linkages of the lines of force at the junctions of the common and non-common 
portions of the filaments. Indeed a line of force emerging from the photosphere 
can change the root to which it returns. Thus for instance a line of force that is 
twisted into a right-handed screw on its way out from the photosphere toa junction 
point can return to the photosphere along a left-handed screw. Many of the 
turns of such a twisted line of force can then simply run up to the junction point 
and untwist themselve.. These considerations also suggest that filaments can 
bifurcate, that not all the lines of force emerging at a particular root need re-enter 
the photosphere at one unique other root. 

g. A remark on the structure of the solar magnetic field.—Throughout the above 
discussion we have thought of a twisted filament as being built up in two stages: a 
first phase in which the filament emerges from the photosphere, and a second 
phase in which the filament becomes twisted by mechanical motions at, and 
below, its roots. Nothing in the main argument would be seriously affected, 
however, if the filament were to emerge from the photosphere in an already twisted 
condition. 

This prompts the thought that the phenomenon of the flare may be an essential 
part of a process whereby internal solar fields become untwisted. Itis well known 
that the Sun possesses a deep and highly active convection zone. (Not only are 
there strong theoretical reasons for expecting such a zone, but many stars not very 
dissimilar from the Sun—dwarfs somewhat further down the main-sequence— 
possess structures that can only be understood on the basis of the presence of a 
convection zone.) The magnetic field within this zone must tend to become 
seriously twisted, so much so that it has always been something of a mystery to 
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understand how the twisting comes to be limited in its complexity. Since it 
seems clear that the twisting does not inhibit convection, one possibility is that 
both the convection and the magnetic field assume highly correlated patterns that 
are specifically arranged to prevent any overall increase in the complexity of the 
field. A second possibility arises from the present discussion, namely that a 
balance is reached through an untwisting of the lines of force during discharges 
in the solar atmosphere. On this basis it would scarcely seem to be by chance 
that the somewhat peculiar physical conditions in the chromosphere happen to 
be favourable to the occurrence of flares. Rather should the chromosphere then 
be thought of as an important factor in the structure and equilibrium of the whole 
Sun. 

The authors understand that somewhat similar considerations to those given 
above are being published by P. A. Sweet. The reader is referred to the following 
papers already published by Sweet (J.4.U. Symposium on Electromagnetic 
Phenomena in Cosmical Physics, No. 6, 123 (1958); and Nuovo Cim., Supp. 8, 
Ser. 10, 188 (1958).) 


Our thanks are due to a referee whose comments enabled the paper to be 
clarified at a number of places. 


Harvard College Observatory, St John’s College, 
Cambridge, Cambridge, 
Mass., U.S.A.: England, 


1959 May. 
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Summary 


In the first part of the paper the distorting effects of the Earth’s atmosphere 
on the distribution of light at the extreme limb of the Sun are considered. An 
analytical expression is found which describes the resultant redistribution 
of light in the focal plane of the telescope; the effects of “‘ seeing’’ are 
represented by a Gaussian term, particle scattering and non-specular reflection 
by the telescope optics are included in a power-law term. This redistribution 
function forms the kernel of the integral equation relating the true and 
observed intensities at any point. A geometrical idealization reduces the 
two-dimensional equation to an equation in one dimension. Methods of 
solving this equation which do not presuppose an analytical form of the true 
distribution are investigated and it is shown that the solutions must be 
regarded as discontinuous. 

In order to apply these methods, observations of the K line of Ca 11 were 
made at points within +10” arc of the limb, on two days. The relevant 


parameters in the redistribution function were determined from the observa- 
tions. A numerical application of the solution of the integral equation then 
gave the true line profiles. The accuracy attainable in the solution is shown 
to depend rather strongly on the width of the kernel and hence on the 
atmospheric conditions. 





Introduction 


During its passage through the terrestrial atmosphere the distribution of light 
over the solar surface is falsified. ‘Time-variation in the refractive index of the 
air along the light path produces a characteristic tremor of the solar image, while 
even in its absence molecular and particle scattering together with non-specular 
reflection by the telescope mirrors transfer light from each point to the surrounding 
ones in the image plane. Distortions in the chromospheric intensity distribution 
near the limb due to these causes are particularly severe, since the first results in a 
blurring over distances comparable with the vertical extent of the chromosphere 
and the second in the superposition of a considerable amount of photospheric 
light on the much fainter chromosphere. This latter effect is absent during an 
eclipse, when the photospheric light is occulted before reaching the Earth’s 
atmosphere. There are, however, many advantages in also securing extra- 
eclipse material, and the practical possibility of observing the chromospheric 
emission lines in full sunlight was demonstrated by Hale (7) at the turn of the cen- 
tury. A great reduction in superposed photospheric intensity can be achieved, 
in this case, by using those emission lines which occur in wide absorption lines 
atthesolarlimb. Nevertheless, relatively little photometric work has been carried 
out, because of the difficulty of applying an accurate correction for the terrestrial 
distortion of the intensity distribution. 
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The mathematical formulation of the basic problem is quite simple. Let 
1(€, n) be the true intensity at the point (£, 7) on the solar image, and let O(x, y) be 
the observed intensity at the point (x,y). If R(r) is the fraction of light from unit 
area at (¢, 7) which is transferred to unit area at (x, y), at a distance 


r=[(x-€P+(y—n)*}" 
from its point of origin, the relation between the true and observed intensities is 
given by 


O(x,y)= i) f I(é,n) R(x—§y—) dé dn. (1) 
I) : 


R satisfies the normalization equation 
JJ [ Renagan~t, @) 


since the light can only be re-distributed. R varies only very slowly ohn wave- 
length and is quite independent of it over a range comparable with the width of a 
spectral line. 

If the photospheric surface is represented by the semi-infinite plane 
— ©<x<o, and all isophotes are supposed parallel to the y-axis, the intensity at 
wavelength A at any point, which may be taken as one on the x-axis, is 


Os,0)= [" 1(6)4& [”_ Ree—E.n)do. 3) 


This represents the intensity distribution along a radius of the image. 
Performing the integration with respect to 7 symbolically by writing 


S(e—e)= [" Ree-E.n)dr, (4) 


where S is the redistribution function for the x-direction, the problem is seen to be 
that of the solution of the integral equation 


Os)= [" 16) Sex-eae (5) 


for J at each wavelength, where O,(x) has been written for O,(x, 0). 

A variety of approximate solutions of this equation has been made. The 
earliest is that due to Elvey and used by Keenan (9) in 1932. Equation (5) may be 
written 


O,(s0)= [" 1,6) SoA ae + [" 116) Ste—E)dE 

=J, +d, (6) 
where J, represents the contribution to the observed intensity from the intensity 
inside the limb and J, that from the chromosphere beyond the limb. Clearly a 
partial correction can be effected by subtracting J,. If the profile of the line 
inside the limb is independent of position (i.e. the limb darkening is independent 
of wavelength), the values of J, at different wavelengths differ only by a scale 
factor; also, outside the emission line J,=o0. The correction then involves 
subtracting, from the observed profile at a given height beyond the limb, the 
absorption line (standard profile) found inside the limb, reduced in the ratio of the 
measured intensity of the continuum at the given height to that at the position 

R* 
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where the standard profile is determined. Owing to its simplicity this method or a 
variant of it has frequently been used (6, 11, 12, 28). In all such cases, quite 
apart from the approximation involved in J,, it is the integral J, in equation (6) 
which is found and not the desired true distribution. 
This point was considered by R. O. Bishop (1). He worked with equation 
(5) in the form 
0 oo 
O,(*)= | TP (€) S(x—€) dé + | T(E) S(x—€) dé 
=J'+dy’. (7) 


J? is the photospheric intensity, which is zero for £ > 0, and hence J,’ is the photo- 
spheric contribution to the observed intensity. J,’ is the contribution from the 
chromospheric intensity J°, which is not zero for £ <0, since the line of sight from 
the observer to a point on the photospheric surface passes through the chromo- 
sphere. J,’ was evaluated by Elvey’s method, the standard profile being deter- 
mined well inside the limb, where the contribution of J° to the observed intensity 
is much reduced. The equation was then solved for /°, taking a Gaussian form 
for S in J,', under the assumption that the chromospheric emission decreases 
exponentially with height above the photosphere. 

The basic weakness of Elvey’s and associated methods lies in the assumption 
that the line profile is independent of position. This is least likely to be so near the 
limb, and the true form of the profile is unknown in this region (11), from which the 
greatest contribution toJ, in equation (6) comes. This was realised by Michard* 
(16, 17), who consequently attempted to solve equation (5) directly without 
assuming any specific form for J. However, in so doing he took a Gaussian form 
of the redistribution function S (as Bishop did in his integral J,’). This is known 
to describe quite accurately the time-averaged effect of the seeing (22, 17); it 
does not take account of the other redistribution processes enumerated above. As 
a result, Michard found erroneously large intensities in the outer chromosphere 
and had to apply an additional correction in this region (17). 

It is clear that no completely satisfactory method of dealing with the correction 
for distortion has been proposed. Further work on the problem has been in 
progress at Oxford since 1956. 


1. The redistribution function and solutions of the integral equation 


1.1. The redistribution function.—The primary need in any attempt at a general 
solution of equation (5) is an adequate form of the redistribution function R and 
hence of S. It will be recalled that the redistribution consists of two independent 
physical processes. The time-averaged effect of the seeing as mentioned above is 
represented quite accurately by a Gaussian term, and this has been shown by 
Wanders (29) to account extremely well for the intensity distribution in the con- 
tinuum across the solar limb to a distance a few seconds of arc outside. The 
inclusion of a second term is necessary to take account of the ‘‘tail’’ of the 
observed distribution beyond this, produced by scattering of light in the Earth’s 
atmosphere and by imperfections in the surfaces of the telescopic system. For 


* This was most lucidly stated in a letter to Mr S. V. M. Clube, to whom I am grateful for the 
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this, a power-law term was used by Zanstra (30) and Kinman (10) in their 
investigations of ‘‘long-range scattering ’’ in a stellar image and in the circum- 
solar aureole respectively. 

These considerations led H. H. Plaskett to suggest the following treatment. 
Suppose the redistribution function to have the form 


R(r) = A exp (—17/2a*) + B(1 +1r?/2a*)™, (8) 


where A, B, a and n are constants. If the solar limb is again taken as the line 
x=o, and all isophotes supposed parallel to it, the observed intensity at any 
wavelength A at a distance x measured radially from the limb is, substituting for R 
in equation (3), 


Os0)= [~ [”_ 1e){Aexp{—[ee- OF +9*)/204 | 
+ Bia+ (xP +9*Y/aat}"| dédy. (9) 


Integrating with respect to 7 yields, again writing O,(x) for O,(x,0), 


O,(x)= { ” 1,(€Xav 2A exp [— (x- €)?/2a*] 


+aV 2m BY (n— 3)[0(n)] [1 + (x—£)2/2a*] "9 dé, (10) 


which is the specific form of equation (5) to be solved for J,. The expression in 
braces, the kernel of the integral equation, is the redistribution function, S, for 
the x-direction. 

One relation between the four parameters of R (and hence of S) is furnished 
by the normalization condition (2), which gives 


1 =2ma* [A + B/(n—1)}. (11) 


The remaining three, most conveniently, a, n and B/A, may be found from a 
comparison of the observed and predicted distributions at some wavelength in the 
continuum. Wanders’ investigation (29) showed that a linear law of limb- 
darkening and a discontinuous drop at the limb were a satisfactory representation 
of the true distribution at these wavelengths, i.e. 


T,()=1g(A) (t-—mE),  E<o; 
=0, £>0, (12) 


where /,(A) is the intensity at the limb, and m (the limb-darkening coefficient) is 
independent of wavelength over a small range. Substituting equations (12) in 
equation (10) gives 


O,(x)= i) " _ 1o(A)(1—mé fav 2a A exp [-(x—£)/244] 


+aV 2m BY (n— 4) (T'(m)} [1 + (x—£)*/2a"]-" +49} dg. (13) 
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Effecting the integration with respect to £, one obtains for the predicted distribution 
O,(x) 


O,(x)/Ip(A) = V22a°A (1 — mx) exp (—u*®) du + 4/2m ma®A exp (— x*/2a*) 
ziaJ/2 
a/2 


+ V2 2a°B (1 —mx)'(n— 4) (T(n)]” cos?" ~3 6 d0 
are tan z/a/2 
+ V 20 ma® BY (n— 4) [(n— 3/2) P(n) (1 + x2/2a”)" 32] -1, (14) 
for n>3/2. ‘The unknown limb intensity, J,, also occurring in equation (14), 
can be eliminated by making the comparison between the logarithms of the 
observed and predicted distributions. 

1.2. Solutions of the integral equation.—Equation (5), of which equation (10) is 
a particular example, is a singular equation of Fredholm’s first type, whose solution 
by a method of successive substitutions (without assuming any analytical form of 
the distribution) has been discussed by Burger and van Cittert (3, 4) (in their 
formulation of a correction for spectroscopic dispersion of spectral lines). After 
conversion of the integral equation to one of the second type, the well-known 
Liouville-Neumann method for that type yields a solution of the equation as an 
infinite series of successive approximations, /,, to the desired intensity distribution, 
I, where the rth predicted distribution, O,, is given in symbolic form by O,=I,.S. 
Successive iterates are formed by adding the difference between the actually 
observed distribution and the last predicted one, i.e. J,,,=1,+(O-—O,). The 
process is started by putting /,=O. For non-singular equations, the Liouville— 
Neumann series starting from this or any other approximation converges as r —> 00 
to the unique solution, J, of the equation, provided that this exists and is 
a continuous function (15). This was verified by Burger and van Cittert for 
their particular example. 

The effect on the solution of observational error in O has been discussed by 
Kreisel (13). He showed that, in these circumstances, the solution of equations 
of type (5) is indeterminate to some extent, since an oscillatory function of any 
amplitude with a period small compared with the width of the kernel produces 
only a small change in O; these solutions are superposable on the solution of 
the accurate equation, since equation (5) is linear. However, he found that the 
possible range of smooth solutions of the equations was limited and further that, if 
the true distribution itself is smooth and continuous, the solutions in this range 
lie close to the true distribution. 

But when an attempt is made by Burger and van Cittert’s method to recover, 
from an observed distribution, a true distribution having a discontinuity in the 
variation of intensity with distance, a distortion of the true distribution similar 
to Gibbs’ phenomenon in Fourier expansions near a point of discontinuity (5) 
occurs. The slope of successive iterates J, increases at the position of the 
discontinuity, but the values of J, at adjacent positions have spuriously high values 
on one side of the discontinuity and correspondingly low values on the other. 

Now the drop to zero in the photospheric intensity actually occurs in a distance 
of the order of 0”-2 arc (14) and the chromospheric intensity may therefore rise 
by a factor of up to two in this distance. On the other hand, the width of the kernel 
in equation (5) is of the order of a, the radius of the tremor disk which appears as a 
scale factor in S, in equation (5). This is 1” or 2” arc even in very good conditions. 
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Consequently, the solutions of equation (5) relevant for the present case must 
be regarded as discontinuous. 

When an intensity distribution having a discontinuity is used, successive 
iterates can no longer be formed by Burger and van Cittert’s scheme, since the 
residuals between the observed and predicted curves are due to a combination of 
errors in quite unrelated distributions. But, since the true distribution when 
inserted in equation (5) must produce the observed distribution, a solution can be 
obtained by trial and error. If there is only one discontinuity in the range of 
integration, the distributions on either side being smooth, the residuals between 
the observed and the predicted distribution resulting from some assumed true 
distribution are insensitive at sufficiently large distances from the discontinuity 
to the precise distribution on the far side of it. In these circumstances the residuals 
are a guide to the next approximation, but there remains a zone of indeterminacy 
near the discontinuity (the limb, in the present case) extending to a distance on 
either side approximately equal to a, the radius of the tremor disk, whichdefines 
the half-half-width of the kernel. ‘ 

On spectra taken in good conditions the chromosphere itself shows a reticulated 
structure, so that the positions of further discontinuous changes in intensity may, 
on occasion, have to be located and taken into account. With the admission of 
discontinuous solutions of equation (5) the uniqueness of the adopted solutions is 
called into question. Spurious solutions are conceivable, though some check on 
their authenticity is provided by the necessity of the solutions forming a set which 
is physically consistent from wavelength to wavelength. A scheme for the 
location of the positions of discontinuities (other than the limb) and an ordered 
sequence for the approach to the true intensity distribution are required. The 
examination of the effects on the solution of observational error also needs to be 
extended, 

The attempt to find a rigorous method of converting the observed to the true 
projected intensity distribution near the limb thus raises a number of interesting 
problems. A considerable amount of further work will be necessary before a 
completely satisfactory solution of the problem can be reached. 


2. Application to the K line of Call 

2.1. Observations and microphotometry.—Observational material for an appli- 
cation of the theory of Section 1 to the H and K lines was obtained by Plaskett in 
the early summer of 1956, using the 19°8m telescope (24), which produces a 
solar image approximately 18cm in diameter. Since absorption by the glass 
increases heavily with decreasing wavelength in this region, only one prism of the 
associated spectroscope (25) of focal length gm was used; consequently, the 
dispersion was approximately 1:8Amm~'. The width of the spectroscope slit 
was 0°025 mm and the wavelength A 3900 A was centred on the plate. 

Definitive plates were obtained on two days, May 27 and June 5. Each plate 
carried three chromospheric spectra taken within a few minutes of one another. 
In addition, for plate calibration, two spectra of the centre of the disk were taken. 
One was with a neutral-tint logarithmic wedge (21, 22) before the slit, the other 
with it removed, this latter (‘‘ blank ’’) showing any non-uniformity of illumination 
over the height of the wedge. All exposure times on a single plate were of identical 
duration (go seconds on May 27, 4 minutes on June 5), suitable photographic 
densities being obtained for the centre spectra by placing calibrated neutral-tint 
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filters in the beam before the slit. The plates (Ilford Ordinary N30) were all 
stored for at least 24 hours before they were brush-developed in Parkhurst’s 
two-solution developer (19), which gives fine grain and little background fog. 

The height of the spectroscope slit for the chromospheric exposures was 
12°5 mm, and the solar image was so placed that the slit, which lies in an E~W 
direction, made an angle of 75° with the image radius, the limb intersecting the 
slit at its mid-point. During the exposures this position was accurately maintained 
by means of a guiding plate (26). The position angle of the Sun’s pole given in 
the Nautical Almanac is — 17° and — 14° for May 27 and June 5 respectively, and 
the corresponding heliographic latitudes of the centre of the disk are — 1° and o’. 
Thus the spectra, which were of the upper limb, lie almost exactly at the North 
pole. 

One plate from each day was selected and microphotometer tracings (23) 
centred at the H line and at the K line were made parallel to the dispersion for two 
chromospheric spectra on each plate. High magnification ( x 48-2) and a narrow 
analysing slit (projected width on the plate= 10-4 microns) were used, the plate 
being driven at the slowest speed (3 cm per hour) past the analysing beam, whose 
height was limited too-5; mm. The tracings were made at some 15 consecutive 
positions, separated by o-5mm, on each spectrum. If the observations are 
regarded as typifying a radial distribution, the interval between the positions is 
17-35 arc. 

Calibration curves for the two plates were constructed, from tracings made 
perpendicular to the dispersion across the wedge and ‘‘blank’’ spectrum at 
several wavelengths free from absorption lines of Rowland intensity greater than 
zero, using Page’s (20) values of the wedge constant. The intensity scales of the 
two plates were chosen so that the intensity in the ‘‘ blank’’ centre spectrum was 
the same ineach. The wavelength scale was set up by interpolation between the 
positions of unblended photospheric lines, a least-squares solution having shown 
that the relation between wavelength and position on the tracing was linear. 

Intensities at an average of 35 wavelengths, between AA 3931°5 and 3935°6A 
and between AA 3967-0 and 3971-2 A, were measured on each tracing. However, 
since no smoothing of the records was made before the measurements, the 
profiles showed considerable irregularities. As the time elapsing between the two 
exposures on either plate did not exceed 10 minutes, an average for each plate was 
taken, the position of the limb on each spectrum having been located (the 
procedure is described in Section 2.2). The results for the K line are displayed 
as a function of radial distance from the limb at a number of typical wavelengths in 
Fig. 1. Numerical values are given in Tables Ila and IIb in the columns headed 
O (= observations). 

Since the absorption by the glass of the spectroscope makes it impossible to 
find the intensity variation with wavelength over a range sufficient to determine 
the required level of the ‘‘ continuum ”’ across the K line at the centre of the disk, 
conversion of the intensity scale to absolute units was made in the following way. 
The intensities at the centre of the disk were measured at the four wavelengths 
listed in Table I which lie more than 1 A from the line centre, so that the intensities 
are not disturbed by emission in the line core, and they fall between the small 
absorption lines in the line wings. The mean intensities for the two plates are 
given in the second column of Table I. They have been compared with the values 
found by Houtgast (8), who made a careful determination of the local continuum. 
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Fic. 1.—Comparison of the true and observed i: ity distributions (I(x) and O(x) respectively) 
as functions of distance (x) from the limb, at two . 
Abscissae: Radial distances from the limb in sec of arc. 
Ordinates: Intensities in units of 1-48 X10" erg cm-* sec per unit solid angle per A 
wavelength interval. 
(a) 439354 A(Ki); (4) 3934-4 A (Ky emission small); (c) 439340 A (Ky emission a 
maximum); (d) A 3933°7 A (K;). 
Unfilled circles and dotted curves refer to the epoch 1956 May 27-493 U.T. Filled circles and 
full curves refer to 195 June 5-579 U.T. 
In order to avoid confusion between the I(x) and O(x) curves in (c) and (d), the curves are 
displaced vertically. The I(x) scale is shown on the division between the two parts of the figure; 
the O(x) scale is the one on the outer borders. 
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This yields a value of 311 on the present scale for the continuum at the centre of the 
disk at K. Interpolation at A 3933-7 A in the table of monochromatic intensities 
of the continuum, at the centre of the disk in the normal direction, given by 
Minnaert (18) yields 4-59 x 10! erg cm~* sec~! per A wavelength interval per 
unit solid angle. Thus multiplication of the present intensities by 1-48 x 10!” 
converts them to these absolute units. 


Tasie | 
Wavelength Intensity 
A 
3928'9 100°9 
3931°5 57°4 
3935°0 47°2 
3936°3 69°6 


Plates for the determination of the horizontal apparatus function of the 
spectroscope in its one-prism form, from the mercury line A 4046 A, were taken 
during the season when the solar plates were obtained, using a '® Hg isotope 
lamp as source. Details of the method are given in an earlier paper (25). The 
width of spectroscope slit, type and processing of the plates and conditions for 
microphotometry were all the same as for the solar spectra. The pattern has a 
half-half-width of o-c41 mm and a well-resolved satellite centred 0-136 mm to the 
red of the central maximum and attaining an intensity 13-5 per cent of it (cf. 2, 25). 
With a mean dispersion 1-77 Amm_" in the region of H and K, the half-half-width 
is 0073 A. Since the H and K lines and their emission cores are intrinsically 
broad and since the interest in the present work was in the correction for terres- 
trial distortion, no correction for spectroscopic pattern has been applied, especially 
as its calculation by the method of Burger and van Cittert (3, 4) (Section 1.2; 
in this case x is the wavelength and S the apparatus function) is well known and 
straightforward. However, there are in the red wing of H a number of medium- 
strong lines of photospheric origin, which, being relatively sharp, will have 
been much distorted by the spectroscopic pattern. The measurements of the 
H line have therefore not been used, except those from the blue wing in the 
determination of the parameters of the scattering function to be described in 
Section 2.2. 

2.2. Determination of the parameters of the redistribution function. —Owing to 
the impossibility, already mentioned (Section 2.1), of locating the continuum 
on the plates, equation (14) for the variation of intensity across the limb had to be 
applied in the wings of the ‘H and K lines. The basic assumption concerning 
the true distribution, equations (12), is presumably applicable provided that 
wavelengths outside the emission core are considered. At wavelengths more 
than 1 A from the centre of the line, in both wings of K and the blue wing of H 
(its red wing is disturbed by emission in the He line), curves of the logarithm of 
the observed intensity distribution, log O,(h), plotted against height, A, in the 
spectrum were parallel. Since the logarithm is increasingly sensitive as the 
value of its argument decreases, it appears that emission beyond the limb is 
absent in these wavelengths and a mean log O,(h)—A curve was formed from 
values at two or three wavelengths for each of the three wings. 

The value of m, the limb-darkening coefficient at wavelengths in the line, 
now appears as an extra parameter, bringing to four (a, n, B/A and m) the number 
to be determined. Preliminary calculations with forms of the redistribution 
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function other than equation (8) indicated a value of m= o, and this was adopted. 
Theoretical curves, log [O(x)/I,]—x, were constructed from equation (14) for 
various values of a with B=o (a Gaussian distribution) and also with A=o 
and some values of n (redistribution following a pure power law). An observed 
distribution lay between the Gaussian distribution and a particular power law 
distribution at some value of a. The values of a and nm so found were quite 
definite, and with these a family of curves having various B/A values was con- 
structed and the value of the best fit with observation noted. (The position 
of the limb, the value of h for which x= 0, was thereby located.) A comparison 
of the parameters showed that the a and nm values were the same in all three wings 
of the two lines, a satisfactory result, and that they were the same for the two 
spectra on the same plate. The B/A values varied a little from wing to wing in a 
given spectrum and by a similar amount between the two spectra on a-plate. 
Since this ratio is sensitive at low intensity levels, where irregularities in,silver 
grain distribution on the plate have their greatest effect, a mean B/A for each 
plate was taken. Values adopted for the plate of 1956 May 27 were a=4” arc, 
n=2°0, B/A=2:25; those for the plate of 1956 June 5 were a= 3” arc, m= 1°75 
and B/A=1-'5. From the normalization equation (2), 4=1/1047 and 1/547, 
respectively. Confidence in the significance of the results is strengthened by 
the fact that the three limb positions determined from the individual wings of 
the H and K lines during the curve fitting were identical for any one spectrum. 

As a numerical solution for the true intensity J in equation (10) was envisaged, 
S (the redistribution function given in braces in the equation) was evaluated at 
intervals of 1” arc from €=0” to = 30” with the parameters found for each plate. 
A smooth curve was drawn and the mean values were found for 1” arc intervals, 
from the first, centred at £=0", out to that at€=10". Beyond this out to 22”'5, 
steps of 2” arc were taken, followed by one step of 5” at each end. ‘The equal 
areas under the curve in the two tails of the distribution beyond 27-5 were 
determined by integration. 

2.3. Solution of the integral equation.—The values of the parameters having 
now been determined, at wavelengths in the absorption wings (K,) where there 
is no emission outside the limb, it only remains to find the values of J,, the limb 
intensity, in equation (14) which give the best distribution of residuals from the 
observed values. This was done for wavelengths AA 3932°0, 3932°4 and 34935°4A. 
The J, values from the two plates at any one of these wavelengths did not differ 
by more than 0-3 unit (3 units are approximately 1 per cent of the continuum 
at the centre of the disk, cf. Section 2.1); the values for A 3935-4 A are given as an 
example in Tables II(a)and 1I(b). This agreement between the plates is not forced 
and is a satisfactory result. The non-systematic nature of the residuals between 
the observed and predicted distributions (illustrated in Figure 1 (a) for A 3935-4 A) 
reflects the adequacy of the analytical form of the redistribution function. Their 
size (which seldom amounts to more than 0-5 unit) indicates the size of the 
residuals which must be tolerated in the solutions at other wavelengths. 

For wavelengths nearer the line centre, AA 3933°0 (0°1) 3934°4A, the observed 
imensities were plotted as a function of distance from the limb (the distributions 
at three typical wavelengths are shown in Fig. 1(5)-(d)), and a curve drawn 
through the points. Intensity values at intervals of 1” arc, corresponding to the 
division adopted for the core of the redistribution function, were tabulated over 
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the radial distance covered by the observations, approximately 10” arc inside 
and outside the limb. These were treated as the observed O,(x) values in the 
ensuing solutions of equation (10). 

A first solution was made by the method of Burger and van Cittert (Section 
1.2), except that, in order to increase the speed of convergence of the solutions, 
the solution at any one wavelength (multiplied by a suitable factor) was used as a 
first approximation at an adjacent wavelength. Thus the distribution at 
A3935°4 (or A3932°0)A, already found, was used at A3934°4 (or A3933°0) A, 
where there can be only a little emission outside the limb. A solution by Burger 
and van Cittert’s method converged rapidly to give an adequate representation 
of the observations. The corresponding true distribution was used as a first 
approximation at A 3934°3 (or A3933°1) A and so on through the line. The true 
distributions so found for both plates at all wavelengths showed an intensity 
minimum near the limb and a maximum some seconds of arc beyond the limb. 
Kreisel’s work shows the impossibility of attributing such a distribution to 
genuine fine structure in the chromosphere, since the widths of the redistribution 
functions at half intensity are 9” and 7” arc on the two plates. The form of the 
distribution was eventually attributed to Gibbs’ phenomenon as described in 
Section 1.2. 

Accordingly, a second solution was made, which took into account the 
discontinuous change in intensity at the limb. As a first approximation, at any 
wavelength, the intensity was taken to be constant inside the limb, and a linear 
function of distance outside, this latter form of variation being suggested by the 
observed distribution at wavelengths near the line centre (see Fig. 1 (d)). The 
value adopted for the intensity inside was approximately that of the point furthest 
from the limb. At most wavelengths the values adopted outside the limb were 
those obtained by extending the straight-line portion of the observed curve. At 
wavelengths well away from the line centre, however, where the intensity outside 
the limb is low compared with that inside, the observed distribution outside is a 
poor indication of the true one and it is much quicker to begin with values indicated 
by the run of the solutions from wavelengths nearer the line centre. 

At a given wavelength the first approximation J, was tabulated for the same 
positions as the observed distribution. The numerical integration according to 
equation (<) was then performed to produce the first predicted solution O,; it 
consists of summing the relevant products of J, and S. O, and O were then com- 
pared. A second approximation /, was then chosen, from the run of the residuals 
at some distance from the limb (see Section 1.2), and integrated and the predicted 
O, compared with O. The process was repeated until O was satisfactorily 
reproduced (within the likely error of observation) by the predicted values. The 
corresponding / values were then adopted as the true intensity distribution. 

One further detail of the computational process must be recorded. Since the 
range of integration is infinite, values of the rth approximation J,(x) have to be 
assumed beyond the range of x-values where the resultant O,(x) can be compared 
with the measured intensities. When this occurred outside the limb, a linear 
extrapolation of the intensity to zero was used. Inside the limb, an extrapolation 
was always necessary, and a constant value, that of /,(x) for x between 8” and 10” arc 
from the limb, was taken; let it be denoted by /,(—5). The contribution to the 
integral in equation (5) for this unobserved region is then simply 


1(—b) [ S(x—£) dé, 
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Fic. 2.—True line profiles of the K line as a function of radial distance from the limb. 
Abscissae: Wavelengths in A. 
Ordinates: Intensities in units of 1-48 X 10"* erg cm=* sec" per unit solid angle per A 
wavelength interval. 
(a) 1956 May 27-493 U.T.; (6) 1956 Fume 5°579 U.T. 
The distance from the limb in seconds of arc is shown on the profiles. Outside the 
limb successive profiles are separated by distances of 1" arc. 
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which can readily be found from the tabulated values of the redistribution function 
and the previously determined area of its tail (Section 2.2). 

2.4. Results—The results of the solutions for the two plates are presented 
graphically in Figs. 2 (a) and 2 (6), in the form of the contours of the lines at a 
number of positions inside and outside the limb. In Fig. 2 (6) the intensity 
inside the limb is seen to be steady to within 6” arc of the limb and the line is fairly 
symmetrical; thereafter the intensity in the line core decreases as the limb is 
approached. Outside the limb, the lines are centrally reversed in the low 
chromosphere, becoming simply peaked some seconds of arc further out; they are 
asymmetrical. ‘This asymmetry would have been accentuated had the correction 
for spectroscopic distortion been applied. In Fig. 2 (a) the lines show the same 
general characteristics, but are fairly symmetrical outside the limb. Inside the 
limb, in this case, the intensity is almost independent of position except for an 
asymmetrical brightening in a very narrow region centred 6”-5 arc inside the limb. 
Possibly this is caused by a small facula; its line profile is shown dotted. The 
numerical values of the true and observed intensities for the two plates are 
compared in Tables II(a) and I1(b) for the typical wavelengths AA 3935-4 (K;), 
3934'4 (K, emission small) 3934:0(K, emission a maximum) and 3933°7 (K;) A. 
Of the two entries for each wavelength at each position the left is the observed 
value (O) and the right the true one (J). 

These results are depicted in Fig. 1 in the form of “‘ profiles of the limb ”’ at 
the four wavelengths. Here the broken lines refer to the plate of Table II(a) and 
Fig. 2 (a), the full lines to the other plate. The smooth curve through the observed 
points (shown as unfilled and filled circles, respectively, for the two plates) are the 
predicted intensities corresponding to the adopted true intensities marked /(x). 
The fit between observation and prediction is seen to be quite satisfactory. 

It is most clearly seen from Fig. 1 that the variation of chromospheric intensity 
outside the limb does not differ significantly from a linear distribution. As 
explained in Section 1.2 the residuals between the observed and predicted 
distributions at any stage in the calculation do not furnish values for the next 
iterate at positions near the limb. Beyond this, however, there was no evidence of 
a departure from the linear distribution, and so the form of successive iterates was 
not substantially changed. Further calculations showed that the insensitivity of 
equations (5) in the present case is such that an exponential variation would yield a 
predicted distribution not measurably different from that given by the linear form. 
Thus the attempt to arrive at the true form of the distribution cannot be considered 
completely successful in the present case. 

A considerable improvement in the situation would be expected if the width of 
the kernel of equation (10) were reduced. This requires a reduction in the 
size of the tremor disk, which appears as a scale factor in the expression for the 
redistribution function, equation (8). In the practical case, though not in the 
theoretical one, this sets a limit to the fineness of the structure which can be 
recovered from the observations. With a sufficiently small tremor disk, possible 
irregular variations of the intensity with position could be detected ; it would then 
be possible to see if differences, such as that in the apparent slope of the emission 
outside the limb shown in Fig. 1, are wholly or partly due to a statistical variation in 
the fine structure. 

A new series of observations is, therefore, planned. Attention will be paid 
to securing spectra in good sky conditions with the 35 m telescope (27), which 
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produces an image having a linear scale almost twice that produced by the telescope 
used in this preliminary investigation. The associated spectroscope of 11m 
focal length has a grating as the dispersing element so that both the dispersion 


TABLE II 
Observed and true K line intensities 


Wavelength d 3933°7 3934°0 3934'4 


Position O I O I O I 
(a) Mean epoch 1956 May 27°493 Duration 1™ 30% 


—10"°0 109 «=«<TI'l , 176 18:4 
—8-7 108 11°3 ° ° 176 184 
—7°3 Ir'O) 6114 . : 17°93 184 
—6°0 1rrO) =OII"4 , 160 18-4 
—4°6 106 11°2 ‘ 15°4 184 
—3°3 10°2 «O«II'l : 1472 18-4 
—1'9 98 110 : 126 18-4 
—0o'6 go 610°9 ‘ 10°8 184 
+08 88 10°71 ° 88 36 
+2°'1 8-6 . . . 70 24 
+3°5 70 . . . 48 1'2 
+48 6:2 . . . 3°5 oo 
+6°2 5‘2 , : ‘ 26 oo 
+7°5 3°83 . . . o’o 

+8*-9 27 . ° oo 


(6) Mean epoch 1956 Fune 5°579 Duration 4™ oo 


—9"'5 18 12°3 15°6 180 6189 28-7 
—81 1270 123 16°6 17°99 189 28°7 
—6:8 116 =612°3 156 166 171 189 26:0 «287 
—5°4 rr8 = 12°2 15°5 16°4 1770 §=6—18'9 258 28-7 
—4'1 114) ««12"1 14°5 16:2 16:5 189 24°5 287 
—2°7 109 =«II'9 138 15°7 154 189 22:2 28-7 
—14 98 11°5 12°0 86150 12°5 189 181 28-7 


: ‘ 1m! ‘ 14°3 : 18-9 ; 28°7 
oo a na {143 1004 io va { tie 


+1°3 92 99 8 10°5 i 99 =«0'0 
+2°7 8-1 9°0 , 8-8 5°4 16 6-2 oo 
+4°0 78 8-0 , 7°2 4°0 12 4°2 oo 
+5°4 6-2 66 . 5°4 28 o8 30 oo 
+6°7 55 sl . 3°6 22 of 24 oo 
+81 39 4806 -3°S ° 1'5 19 oo 20 8=60'0 
+9°"4 28 20 : roe) oe) oo 
Intensity unit: 1°48 x 10'* erg cm~* sec™! per unit solid angle per A wavelength interval. 


and resolving power will also be improved. The width and shape of the 
spectroscopic pattern should be more favourable, with a consequent decrease in the 
spectroscopic distortion of the line profiles. 
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A GENERAL FORMULA FOR THE CALCULATION OF 
ATOMIC PHOTO-IONIZATION CROSS SECTIONS 


A. Burgess and M. }. Seaton 
(Received 1959 May 5) 


Summary 


The general formula is derived by considering the model of a single 
electron moving in a central field. Approximate bound-state radial functions, 
accurate for large radial distances, may be obtained once the effective 
quantum numbers »(=n*) are known. Bates and Damgaard have shown 
that such functions may be used to obtain good estimates for bound—bound "' 
transition integrals. For bound-free transitions use is made of approximate | 
free-state radial functions having exact asymptotic forms, their phases being 
given by 5=7p where yp is the extrapolated quantum defect (u=n-—v). 

The results of extensive numerical calculations are summarized in tables 
which permit the rapid calculation of transition integrals once the energy 
levels are known. Both bound—free and bound—bound transition integrals 
may be obtained. For bound—bound transitions v+v’ good agreement is 
obtained with the Bates and Damgaard tables for ted I'S. 

Comparisons are made with other results for bound—free transitions, as 
obtained both from theory and from experiment. In nearly all cases the 
general formula gives results at least as accurate as those obtained in the best 
alternative methods of calculation. 





Introduction.—In order to calculate the intensities to be expected in a recom- 
bination spectrum one requires many radiative transition probabilities and 
recombination coefficients. Most of the bound-bound transition probabilities 
required may be calculated using the tables of Bates and Damgaard (1). In 
order to calculate the recombination coefficients as functions of electron tempera- 
ture one requires the corresponding photo-ionization cross sections as functions 
of the kinetic energy of the ejected electron. 

Bates (2) has given formulae which may be used for the calculation of photo- 
ionization cross sections when the bound-state wave functions are known and 
when regular Coulomb functions may be used for the state of the ejected electron, 
Although these formulae are very useful for certain special problems they cannot 
be used for many cases of practical importance. Vainshtein and Yavorsky (3) 
have attempted to allow for distortion of the continuum wave function by using 
an effective charge Z,, but they give no general rule which may be used to 
determine the value of this effective charge. Their method may be criticized on 
the grounds that their continuum functions will have incorrect asymptotic form. 

Our method of calculation is similar to that used by Bates and Damgaard 
for bound—bound transitions. The asymptotic forms of the wave functions, for 
both the bound and the continuum states, are determined from the known energy 
levels of the initial atomic system. We use the theory of the Quantum Defect 
Method described in a previous paper (4), which will be referred to as QDM. 


9 
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Our results may be used to calculate probabilities for both bound-free and bound- 
bound transitions. They are found to be consistent with the Bates and Damgaard 
tables and enable similar tables to be obtained for larger values of the effective 
quantum numbers. 

In later papers we shall apply the results of our calculations to the problems 
of interpreting recombination spectra in gaseous nebulae and of developing a 
quantitative theory of ionization equilibria in nebulae. Our results should also 
prove to be useful in stellar opacity calculations. 

2. The quantum theory of atomic photo-ionization.—We consider ionization of 
an N-electron atom* by light quanta of energy Av. The quantal formula for the 
cross section is 
_ Sartety 


2 
= [¥eGnx .--Xyla)RY (X,,X2,...Xyli, E) dr (1) 


where » is the frequency and 


a- 


R= 5 r;. (2) 


j=1 

The initial atom wave function is ‘Y(x,,X,,...Xy|a) where x;=(r;,¢;) is the 
space and spin co-ordinate of electron j. The final wave function for the system 
of ion plus ejected electron is ‘Y’(x,,X,,...Xy|i, Z), E being the electron kinetic 
energy. “The energy conservation condition is hiv = 1+ E where J is the threshold 
ionization energy. The initial wave function is normalized to unity and the 
final wave function to 


| We(i, E’)¥(i, E)dr=3(E-E’), (3) 


8(E—E’) being a Dirac delta-function. The summation >’ in (1) is over all 
initial and final states of fixed energy and w is the statistical weight of the initial 
energy level of the atom. 

Equation (1) may be re-written 


° 


a, = Al (+H) 28 (4) 


where 
2 
$= | We(a)R¥ (i, k*) dr (5) 


and where « is the fine-structure constant and a, is the Bohr radius. In (4), k® 
is the electron kinetic energy and both J and &? are in Rydberg units (13-60 eV 


or 109737cm™~'). ‘The quantity S is in atomic units and the normalization 
condition for ‘Y(i, k*) is 


| W* (i, k’)¥ (i, 2) dr = 08(k2— ke). (6) 


The atom wave function ‘¥(a) is anti-symmetric in the co-ordinates X,,X,,...Xy 
and the ion wave function ¥'(X,, Xq, . ..X,_,|i) is anti-symmetric in the co-ordinates 
X,Xq,---Xy_y- For ¥(i, k*) one should use anti-symmetric functions of the type 


F(X) Xg ...-Xy]i, k*) 


N 
=N-” 2, (— 19 -PP (xy, Key « « -j_gy Hj 41, + oe Ryley (%)). (7) 


* The initial system is referred to as the atom. It may, of course, be an ionized atom. 
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If :, is of the form 
U(X) = Ying (#)8(ol m,’)~Gue(?), (8) 


where Y;,,, is a normalized spherical harmonic and 5(o|m,') a normalized spin 
function, the normalization condition (6) requires that G,, should have asymptotic 
amplitude k-*? (5). In (5), 2’ includes a summation over /'m,'m,’. 

Assuming LS coupling one may take the atom quantum numbers to be 
aSLM,M, and the ion quantum numbers to be «"S"L"M,"M,". For the final 
state it is convenient to use functions ‘V'(a"S"L"ki’ S'L' M,' M,') where S’L’M,'M,’ 
refer to the complete system; such functions may be obtained from (7) and (8) 
using vector-coupling formulae. The radial functions G,, will then depend on 
a” S"*L*kl'S'L’ but will not depend on M,’M,’. The usual selection rules apply: 
S=S', M,=M,', L=L’', L'+1, M,=M,', M,;'+1. For S/w we obtain 


8 : > 
w (28+ 1)(2b+ 1) MoM PLM’ 





“3 
| Y*(2SLM,M,)RY («”S"L’kl'SL'M,M,')dz\ . 
(9) 


We consider the usual approximation in which the atom and ion wave 
functions are built up from one-electron orbitals and we assume yf, to be ortho- 
gonal to all the orbitals for the atom. We then obtain the selection rule that 
only one-electron jumps occur*. Let the nl electron be ejected and let the 
threshold ionization energy be J,,. The angular momentum quantum number 
of the ejected electron is then /’=/+1. From (g) one obtains 


S 


o a E P(r)rGy(r) dr] (10) 


Ww Fell 


where P,,(r) is the radial function for the nl electron. If more than one value 
of L’ is allowed by the selection rules the right hand side of (10) should include 
a summation over L’. Let n"l” be an electron which is not ejected. We assume 
the radial functions P,., to be the same for the atom and for the ion; this will 
introduce only a small error in the calculated cross section (6, 7). With this 
approximation the C, in (10) are algebraic factors obtained from the integrations 
over spin and angular co-ordinates. These may be evaluated using standard 
atomic structure methods (8, 9, 10). We consider the three cases of greatest 
practical importance. 

Firstly we consider an initial configuration consisting of closed shells 
together with an electron ni. For ejection of n/ one obtains for C,, ,: 


C(Lol + 1) = (1+ 1)/(2l+1) 
C(lol— 1) =1/(2l+1). (11) 
Secondly we consider an initial state 2S”L"nlSL where aS"L” specifies a 
parent term. One obtains the selection rule that transitions do not occur to 


states aS"L” of the ion with (S", L")#(S’, L”). For transitions to states 
aS”L” of the ion one obtains for C;: 


C(S"L'ISL+S’L'l' SL’) =(al+1)(2L'+1)WAWLL'; 1L")C(lel’) (12) 


* One may obtain finite probabilities for transitions violating this rule if one allows for differences 
between the orbitals for the atom and the orbitals for the ion. 


9* 
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where C(I’) is given by (11) and where W is a Racah coefficient. Expressions 
for W? are given in Table I. Using a sum rule for W? one obtains 


> C(S*L'ISL+S’L'I' SL’) =C(I-+1). (13) 
- 


In (10) one may therefore use C,=C(/-l') if the radial integral is effectively 
independent of L’ or if only one value of L’ is allowed by the selection rules. 


Tasie I 
L’ W%U+1LL’; tL") 
oe (L" +L —I(L" +L —l-1(L"—L+14+2L"-—L+1l+1) 
4(2L + 1) 2b — 1) L(+ 1)(2l + 1)(2l + 3) 
(L°+L+l+2)—L"+L4+l4+1)(L"—L4+14+1(l"+L-J) 
4L(2L + 1)(L + 1)(2l+ 1) + 1)(2l + 3) 
(L°+L+14+3KL"+L4+l+2)(-L°+L4+1+2-L"+L4+1l4+1) 
4(2L + 3b +1)(2Lh + 1)(2l + 3) + 12l+1) 
Note that W°(ll’LL’; 1L”)=W*I'1L’L; 1L”). 





L 








Thirdly we consider ionization from a configuration n/¢ of equivalent electrons 
With initial state nl@SL and final state nl*-"S"L"ki'SL’ one obtains 


C(#SL+ik“S"L'I' SL’) 
=| (@SLile—(S"L")ISL)PC(S’L*ISL+S"L'I' SL’) (14) 
where (/¢SL{|#-1(S"L" ISL) is a fractional parentage coefficient. This coefficient 
is unity for g=2. Defining 
C(RSLoikS"L'l') = ¥ C(h#SLoir“S’L'I' SL’) (15) 
Fr 


one obtains, using (13), 
C(MSLohS’ Ll’) =q\(USLIl-(S"L’ ISL)PC(1’). (16) 
Therefore 
C(ESL+i“ SLI SL’) 
=(2/+1)(2L'+1)W°4ULL'; 1L")C(M#SL>kS’ Ll). (17) 


Values (6, 11) of the coefficients C(p*SL-+>p?S" LI’) are given in Table II. 


Tase II 
The coefficients C(p*SL+p*—S"L"s) 

S"L" Cc SL S*L’ 
1S 1/3 “ ‘Ss 
? 2/3 3p "D 
os 2/3 7p ~ i 
2p 2/3 Ip 4S 
>? I ip *D 
1p ° ip ?- 
1S ° 1S 4S 
*p 1/2 1S "D 
1D 1/2 1S “4 
1g ° 4 ? 
eg 1/2 sp Ip 
iD 5/18 4 1S 
1S 2/9 1S 2p 


Note that C(p*SL-.p*-1$"L"d)=2C(p*SL-.p*“8" Ls). 


WwWwwwewwewwwnndn 4 
AUuuutk PPE HH HHS 
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3. Use of the Quantum Defect Method.—We consider the radial equation 


E ~~ Vip) +E |y=o (18) 


where E and V are in Rydbergs and where V(r) = —22/r for r large. We put 
E=e2* and p=rz. (19) 
For E<o solutions of (18) satisfying 


| ° Pr)dr=1 (20) 


exist only for the eigenvalues ¢,, of «. We put «,,= —1/* where v is the effective 
quantum number*. The quantum defect is defined at the eigenvalues by 
p(e,)=n—v where n is the usual principal quantum number. A generalized 
quantum defect yu(e) is obtained on interpolation or extrapolation from the 
values of y(e,,). - 
For «> 0 we put « = k?/z? and define G,,(r) to be the solution of (18) satisfying 
G,,(0)=0 and 
Gir) ~ k-*? sin [x +8(e)] (21) 


where 


x=kr—4lr+ i In (2kr) + arg P'(/+ 1 —iz/k). (22) 
We wish to calculate 
g(ols €P)=Luf” Pail Geel) dr (23) 


where /’ =/+1 and where J,,, =2*/»*, «' = (k'/z)*. We assume that the eigenvalue 
spectrum is known but that the form of V(r) may not be known for small r. 
The Quantum Defect Method enables us to obtain the form of the radial functions 
for large r. Our method of calculation depends on the fact that, for typical 
atomic potentials, g(vl; «’l’) is usually insensitive to the exact form of the radial 
functions for small values of r. Approximate values of g may therefore be obtained 
if we make certain reasonable assumptions concerning the radial functions for 
small r. 


It is shown in QDM that 


Piulr)=2"*K(v, JW, ,,,(2e/v) (7 large) (24) 
where W is the Whittaker function. The normalizing factor is 
K(v, /) = [f(v) 40 (v ++ 1) (v—D)]- (25) 


U(v)=1+ 55). 


When » is not small one may usually put {(v)=1. The asymptotic form of the 
Whittaker function is 


W.1+s (ele) = (2) em” § (vs Dp++O(-+) (26) 
where 5, =1 and 


b= = [Kl+1)—(v—t)(>—t+ 1) (¢>1). (27) 


*The symbol n,j® is often used for the effective quantum number. In QDM we used my or r,. 


where 
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An approximate bound-state radial function is defined for all r by 


P(r) =2"*K(», 1) (72)' em ¥ bio Np~. (28) 


For v=n=(l+1), (/+2)... the expansion terminates at t=n—/—1 and is then 
equal to the exact hydrogenic eigenfunction, Bates and Damgaard (1) have 
shown that, with suitable choice of ty, (28) may be used to obtain good estimates 
of the bound—bound transition integrals, 


[© PailrPwe(r) ar 
0 
We shall consider positive energies E’ such that 

exp[—27/e’ 1?] <1. (29) 

The phase 3’ in Gy, is then given by 
8'(<’) = mp '(€’) (30) 

where yp’ is the extrapolated quantum defect for the n'l’ series. The Coulomb 
functions y,'=y,(ie’'-?,l’; p) and y,' =y,(ie’'-"*,l’; p) defined in QDM are 
solutions of (18) in the limit of large r. Using (29) we obtain the asymptotic 
forms 


; 2 ys 

yy nae (<7) sin (x’) (31) 
, A’ v3 , 

yy! ~ - (27a) cos (x’) (32) 


pew \me te 
where 
A’=(1+el?)[1+e'(l'—1)?]... [1 +e]. (33) 
At the origin y,’ behaves as p+! and y,’ as p~" together with logarithmic terms. 
We define 


F(é',l'; p)= (=)" yy (34) 


and 


22 
where 7, is real and positive. Thefunctions F and H behave as p’ *? at the origin 
and have asymptotic forms F~e'-"‘sin(x’) and H~ —«’-"*cos(x’). As an 
approximate positive energy function we use 

Gyr(r)=274{F(@',1'; p)cosmp’'(e')—H(e',!'; p)sinmp’(e’)} (36) 
which behaves as p' +! at the origin and has asymptotic form 
Gyy~ k-™* sin (x’ + mp’). 
Using (28) and (36) the integral (23) is given by 
a(vl; <l)=f(vl; <l')cosmp'(e’)—A(vl; ¢'l') sin m'(e’), 


H(e1's p)= (1 exp(— ree)" (S) 9 (35) 


where 
fot; et)= =E9 (2)" $ aos” pre*eFlet's p)dp 
d v" V} t=O 0 
an 


v 


A(vl; él')= =) (=) $ bv, nf” p’tite-"H(e', I’; p) dp. 
t=0 0 


It is seen that g, f and h do not depend explicitly on z. 
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Choice of t, and of tj. For small r the radial functions we have adopted depend on 
the choice of t, and of r,. The method of calculation depends on the fact that, 
once reasonable values of these parameters have been chosen, the integrals are 
insensitive to small changes in the parameters*. 
The Coulomb functions y,’ and y,’ are solutions of 

PFE 48) 338 G7 

dp ae ai + ? +é€ p =0. (40) 
With increasing p the first point of inflection of y’ occurs at p=p, where 

- 143), 2 eeee. 
Pr Pr 








- 2 — 
Fic. 1.—The regular Coulomb function, y,' and the irregular function y,' for « =0, l’ =2. 


For the values of ¢’ of interest p, differs little from the value p,=/'(I'+1)/2 for 
e’=0. For p>p, the functions y,’ and y,' are both oscillatory, but for p<p,'’, 
y,' is small and y,’ starts to diverge (see Fig. 1). The cut-off factor 


[1 —exp(—typ)]}¥*** 

multiplying y,’ in the function H is therefore chosen to be close to unity for 
p>p, and to be small for p<p,. The condition we adopt is 7,p,~5 giving 
ty=10/I'(l'+1). This condition suggests an infinite value for 79, corresponding 
to omission of the factor [1—exp(—79p)] in the s-wave function H(e’,o; p). 
This is in fact quite satisfactory since H(e’,o; p) has only a weak logarithmic . 
divergence at the origin. Functions G,,(r) calculated with 7, =10/I'(l' +1) 
have been compared with radial functions calculated by numerical integration 
for various atomic potentials. For the values of p of importance satisfactory 
agreement was obtained. Final calculations of the integrals were made for 
T)= 0, 7,=5, Te=1°5 and 2:0. 

The choice of t, has been discussed by Bates and Damgaard. They consider 
two possible criteria. Criterion (a), which is similar to that used for an asymptotic 
expansion, is that the series should be terminated at the minimum term. Criterion 

* One could obviously consider various alternative definitions of the radial functions for small r. 
Instead of using the asymptotic expansion of the Whittaker function, Hylleraas (12) has represented 
this function by a contour integral and has then deformed the path of integration so as to obtain 
convergent bound—bound transition integrals. This method gives results very similar to those 
obtained by Bates and Damgaard. 
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(b) is to neglect in the integral powers of p less than 2; t) is then the integer 
satisfying 
v+l’—1<ty<vtl. 
Bates and Damgaard found these two criteria to give practically identical results. 
This was also found to be the case in most of our calculations but (a) and (b) were 
found to give significantly different results for (/,/’)= (2,1) and (3,2). In these 
cases criterion (a) was adopted. With exact wave functions the integrand would 
behave as p'+"+# for p small and it is therefore evident that criterion (5) will become 
unsatisfactory for large /, I’. 
4. Evaluation of the integrals.—The integrals f and h may be expressed as 
power series in ¢’: 
fl; eT) =fool; V)+ef,(vl; V+... (41) 
A(vl; el’)=h(vl; ')+ehy(vl; I’)+.... (42) 
Using 
A’ =1+$l'('+1)(2l'+1)e'+... (43) 
together with formulae given in Appendix I of QDM we obtain for the leading 
terms in the expansions of F and H: 


F(<,1; p)= (no)! {Juuss( VB) 


+ S[Kl+ t)(2l+ 1 n4s(V 8p) — 3(1+ 1)(2pW a143(V 8p) 


+ (20)"Ja..4(V8)]} +0(2) 


and 
H(6,1; p)=[1 — ee} mp) { Youi(V 8) 


+ < [l+ 1)(al+ 1) Yo, (V8p)—3(1+ 1(2p) Yara V 89) 
+ (2p)*¥uee(VEp)}} + OC) 


The integrals for f, and f, may be evaluated using 


[e-rrpt amex VB) dp 
_ P(m+q+2)2rtisynrers 
3 T'(2m +2) 
(x3) where ® is the confluent hypergeometric function defined by 
ax  a(a+1)x* 
e311 ee+1)a! °°" 





O(m+q+2,2m+2; —2v) (46) 


@(a,c; x)=1+ (47) 


Using Kummer’s transformation, 

O(m+q+2,2m+2; —2v)=e-*'O(m—g,2m+2; 2v), (48) 
and recurrence relations satisfied by the confluent hypergeometric functions 
(14, 15) one obtains 

R(vi, I’) 


folvls)= Facey Soh?) 


WLP) x ROLLY EU + al +1) »__ AOI, I’) 
fiw, P) sal 12 Sol, F) “6 +1) 


and 
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where 


ba 20 12 grew 
RON =| eT (2l’+1)! 


ty 
S(vl, ')= 2 ahvl, l')d{v,1'), 


I(v+l'+3) e-*, 





O(v1, 1) = ¥ pv}, 1')0,(v1, 0), 
t=0 

_ (Kl+1)-(v—t)v—-t+1)) 
os 2t(l’ +v—t+3) G1 (#21), 

P=(+v—t+4)(' +v—t+3)q, 
$v, l’)=O(l'—v+t—1, 2l'+2; 2»), 

oe ny, UU 4+1)-vel—vt+t- 
6,(vl, 1 )=balv,! )+ (U’ + 1)(2/’ +3) 
The functions ¢,(v,/’) may be evaluated using the series expansion (47). For 


given »v, I’ it suffices to calculate ¢, for two values of t; to obtain ¢, for other values 
of t one may then use the recurrence relation 


(l' +v—t+3)b)_1(v, 1’) +2(t-—2)d(v, ')— (Ul —v+t—1)br(v, l') =0. 
Calculations of f, and f, were made for the following values of the arguments: 
l I’ v 
0'2(0°2)5°0,* 1°0(1°0)12°0 





Jo=1l, 





2) 4. (v,U'+1). 


3 


1'0(0°2)5:0, 
1°0(0°2)5-2, 
2°0(0°2)5°0, 
2:0(0°2)5°0, 
3°0(0°2)5°0, 


6°0(1-0)12°0 
6°0(1°0)12°0 
6-0(1-0)12°0 
6°0(1°'0)12°0 
6-0(1-0)12°0 


Analytic expressions were obtained for hk, and h, but were found to be so 
complicated that it was decided to evaluate these functions by direct numerical 
integration. Whenever Aumerical integration was used for hy, h,, the functions 
fo and f, were calculated both by numerical integration and from the analytic 
expressions. Numerical integration was used to calculate A, and h, for the follow- 
ing values of the arguments 

l I’ v 
I ° 1-0(0°2)3°0 
I 2 1'2(0°2)3°0 


For (i,/')=(0,1) the functions f and A were calculated by numerical inte- 
gration for «’=0°00,0°05, 0°10. The continuum Coulomb functions given by 
Bates and Massey (16) were used, the range of p being extended by numerical 
solution of the Coulomb equation. The ranges of v were 


é’ v 


° 1:0(0°2)3°8 
0°05 1-0(0°2)3"0 
O-10 1'0(0°2)3"0 


* In this range we calculated f, but not /,. 
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5. Behaviour of the integrals.—The functions f(vl; <'l') and h(vl; <'l’) were 
found to be oscillatory functions of v with period tending to 2 for large v. Further- 
more, f and h were found to have phase differences of approximately 7/2. This 
behaviour is illustrated in Fig. 2. We therefore put 

eee See ES co 
f(l; «l= iG Ty cos a[v+x(vl; <l’)] (49) 


and 


h(vl; &l’)= ey inal + x <T’)) (50) 


giving 
g(vl; <l')= rip comale to Ce) + xo e'l')}. (51) 


The functions G and y, calculated from 

G=0%(f24h)2, a(v+ x)=aretan (h/f), (52) 
were found to vary smoothly with v. For v not small it was found that x could 
be represented by two or three terms of an expansion in powers of 1/v and it was 
therefore decided that x must tend to a finite limit for large v. 

The behaviour of the integrals may be understood by considering the limit 
of large v. For p not too large the radial functions for small negative energies 
are similar to those for small positive energies. Let v vary from v, to (v,+2) 
where v, is large. There is then only a small change in the energy parameter 
«=-—1/v*, The main change in the radial function is a change in the positions 
of the nodes, this change being similar to that for a positive energy function 
when ,» changes from yp, to y,—2. From the fact that the radial functions are 
similar for v=v, and v=v,+2 one may understand why y depends little on 
v for large v. In a similar way one may understand why f and A are 7/2 out of 
phase in the limit of large v. 

If f is known as a function of v the equation f= Gcosm(v+ x) does not, in 
general, suffice to determine both G and x. However, if one assumes G to be 
non-zero, x is determined at the values of v for which f=o. There is generally 
one such zero between each integer value of v. If one further asumes that x 
varies smoothly between the values determined at the zeros of f one may deter- 
mine x for intermediate values of v by interpolation and for larger values of v 
by extrapolation. Having calculated x from the zeros of f one may calculate 
h from h=ftanz(v+y). The values of h obtained in this way may then be 
compared with those obtained by direct numerical integration. Such comparisons 
have been made for (/, I’) = (0,1), (1,0) and (1,2). Thecomparison for (J, /’) = (1, 2) 
is shown in Fig. 2. The two methods of calculation are found to give good 
agreement for y>(/+1). For /=2 and 3, A was not calculated by numerical 
integration. 

For v> 5, f, and f, were calculated only for integer values of v, and G was 
then calculated using extrapolated values of x. Since G varies slowly with v, 
values for non-integral v greater than 5 may be obtained by interpolation. 
Energy variations. Values of f(nl; <'l’) may be calculated exactly for all <’. 
Such calculations for* n=1, 2, 3, and 4 show that the linear approximation, 
f=f,+«'f;, is valid only for very small values of «’. We therefore sought to 


* These calculations were made by Dr A. Dalgarno. 





No. 2, 1960 Calculation of atomic photo-ionization cross sections 


3 
fo (vi52) 


2 








ho (v1 52) 








Fic. 2.—The functions fv1; 2) and hofv1; 2). The full line curve for hg is obtained from 
ho=fy tan w(v+y,,3) where x,,_ is obtained from the zeros of fy. The curve - ~- - for hy +s 
obtained from numerical integration using t,=1°5 and the curve using T,=2°0. 
Comparison of radial functions suggests that +, = 15 will usually give the most accurate results. 


express our results in a form which would be valid for larger values of «’. We 
first considered negative values of ¢’. 

The positive energy functions F(e’,l'; p) and the bound-state functions 
P,,,(r) are both proportional to the function F defined in QDM, this function 
being an entire analytic function of the energy. Comparison of the proportionality 
factors (QDM, Section 4), shows that the integral 


(ort) = | ° PalPre(r) ar (s3) 
is given by 
2 eee 
Olrbt)==[ ae | (corere(s Sr) Ge) 
or, using (51), 


oMl= Fl ene] F" 


x(x; S311) covn| v—v' + x( ol 53+) |. (55) 
The Bates and Damgaard calculations show that x ( vl; aa , r) andG (» . ai Jt ) 
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vary smoothly but increase rapidly as v tends to v’ and have large but finite values 
for v=v'. Since we are not concerned with the immediate neighbourhood of 
v= v’ we considered functional forms which are singular for y=’. The following 
were finally adopted : 


G(ul; <l’)=(- sass poe (56) 


x(ul; el’) = xiy(v) + .+ = Av) a <Y ——; Bir(v). (57) 


The functions G,y, yy. Xiys %y and OH were iiiiesined from the results of our 
calculations. In all cases it was found that «,, could be taken to be independent 
of v and in nearly all cases that £,, could also be taken to be independent of v. 
The energy variations obtained give close agreement with: (i) the Bates and 
Damgaard tables for (v’ —v)>1°5; (ii) values of g(vo; «'1) for v= 1-0(0'2)3-0 
and «’=0-00, 0°05 and o-10; (iii) values of f(nl; «'l') calculated exactly for 
n=1,2,3 and 4ando<e’' < 1-0 (iv) values of f(nl; ¢'l’) calculated from asymptotic 
expansions (17), for the ranges of n, /, e’ andl’ over which the asymptotic expan- 
sions are valid. Further checks are provided by the comparisons with other 
calculations (Section 7). 

6. Summary of results.—Let I,,, be the threshold ionization energy and k* 
the energy of the ejected electron, both measured in Rydberg units, and let z 
be the charge on the final ion. Put J,,,=2*/v* and k?=2%e’. The photoionization 
cross section is 


on T+? ‘ye 

a;(nl) = 8-559 x 10 0) > Chall; ¢l’)|? cm? (58) 
Tu Ueltl 

where 

py Gl; el 12 

gels et) = AD conav+u'(e)+x0ols PY, (59) 


G(vl; el’) =(-1G,(v) [1 +4], (60) 





ee , ‘2 
x(vl5 eT) = xur(v) + ae + 5 Bl) (61) 


and where p’(e’) is the extrapolated quantum defect for the n'l’ series. Expressions 
for the algebraic factors C, are given in Section 2. The quantity 
2 Ou(e) 

{(v, l= I+ - a ae. ’ (62) 
where p(e) is the quantum defect for the n/ series and «= —1/v*, is usually 
close to unity when v is not small. For vl, G,,(v) behaves as (v—/)"* and {(v, /) 
behaves as (v—/); therefore {-"G,, tends to a constant. 

Results for v>(1+2). One may take B,, to be independent of v and one may 
represent x,,(v) by 

Xu(¥) = Ay + by ]/v + Cyr. (63) 
The coefficients @,,, b,y, Cy, &y and B,, are given in Table III, and the functions 
G,(v) and y,,(v) are given in Tables IV and V. These tables give values of G,, 
and y, for v>/+1; they may be used for v between (/+1) and (/+2) but will 
give results slightly less accurate than those obtained from Tables VI to VIII. 


Whenever Tables IV and V are used x,, should be calculated from (63) and 
Table III. 
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Tables III, IV and V have been published previously in a short report on 
the present work (18). 
Taste III 


Cw by 
+0°2515 —0'078 0°000 
—o'171 0*000 . 0°000 
+0600 0°000 +0°0535 
—0°272 0*000 —o'o19 
+1°170 0000 +0°106 
+0°599 —2°432 +0°050 


Taste IV 
Gy) 
(1, ')=(0, 1) (1, 2) 


2°723 ie 

2°095 2‘840 
1°856 ° 27264 
1°718 ° 2°010 
1°623 : 1°856 
1°553 : 1°749 
1°498 , 1666 
1°452 ° 1601 
1°414 : 1'546 
1°381 ° 1501 
1°352 : 1°461 
1°327 ’ 1°427 


Cn Aut wnreet 


TaBLe V 
; yur) 
(1, l')=(0, 1) (1, 2) (2, 1) (2, 3) 


1°754 ose 

1°605 1'°574 exe oon 

1°§91 ‘ 1°582 1°819 1°447 

1*590 : 1°579 1-771 1°535 

1°591 . 1°582 1°741 1°544 

1°594 1°587 1°722 1°549 

1*596 A 1°593 1°707 1*556 

1°599 A 1*598 1°697 1*564 

1°601 1603 1688 1°573 

1603 1608 1682 1°581 

1-605 , 1°614 1676 1°589 1°926 
1°607 ‘ 1°618 1672 1°596 1'928 


e~x Auhwndvre ct 


Results for v small. The various (I, I’) cases are considered separately. 
ns->k’p 


‘ ev : 
x(vo; €'1)=x—(v)+ Tadly ORO 


Gop You 2d Xo; are given in Table VI for y= 1-0(0-2)2°0. 
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Taste VI 


v Goilv) Yor(v) xo1(¥) 
06 3°259 1°85 +0°143 
08 2976 1°77 +0°085 
1’o 2°739 1°701 +0°043 
1°2 2°527 1°655 +o°o1! 
1*4 2°360 1°632 — 0008 
1°6 2'244 1°620 —0°020 
1°8 2°162 1°612 —0°O3I 
2°0 2*095 1°604 —O'O4I 


np+k’s (v1; ¢’o) is independent* of ¢’ for all <’ and y,9(v) is independent of 
v for v2. Table VII gives Gio(v)/(v—1)** for v=1-0(0-2)1°6, and G(r), 
Yi0(¥) and x49(v)=x(v1; €’0) for v= 1-0(0°2)3-0. 


Tasie VII 


Gil) Y 10) xX 10() 


0°000 1°333 — 0°330 
0°670 1°55 —0°321 
0°826 1°585 —0°313 
o'9gII 1°630 —0°306 
0962 1°655 —0°300 
0'999 1°667 —0°295 
1°029 1°667 —0°290 
1058 1667 — 0°286 
1080 1°667 —0o'281 
1*100 1°667 —0°277 
I°rl7 1°667 —0°273 
np—k'd 


, 
€v 


ey? 
aap x 0°362 + Trevi Baal”) (65) 


x(vt5 €'2)= xy2(v) + 


Table VIIE gives G,.(v)/(v—1)** for v= 1-0(0-2)1-6 and G,,(v), y49(v), X12(v), 
B,2(v) for v= 1-0(0-2)3:0. 


Tasie VIII 
v 

ae G,,(v) Yia(v) Xi2(¥) B,2(v) 
5°69 0°000 2°340 0650 0'079 
5°57 2°489 I°QII “511 “069 
5°02 3°174 1*703 “389 "054 
4°25 3°291 1°625 287 038 
xen 3°O75 1°624 ‘210 "029 
2°757 1658 “164 "035 
2°512 1°675 "1425 053 
2°415 1°635 *131 *068 
2°386 1°593 "IIS 068 
2°340 1°576 "0945 ‘060 
2°25! 1°597 0°073 0°050 


* a39 and 810 being both zero. 
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7. Comparison with results derived by other methods.—For bound—bound 
transitions Bates and Damgaard have compared the results of their calculations 
with many other results, obtained both from theory and from experiment. They 
conclude that their method gives remarkably accurate results for simpler systems 
and, for complex systems, results which are useful and in some cases precise. 
Our calculations are based on a further development of the Bates and Damgaard 
method and enable the Bates and Damgaard tables to be reproduced for | »—v'| 
not too small. Our results for bound-bound and bound-free transitions should 
be of an accuracy not inferior to that of the Bates and Damgaard results. 

In assessing the accuracy of calculated transition probabilities it is necessary 
to consider whether the radial integral is sensitive to small changes in the wave 
functions. When such sensitivity occurs our results will be sensitive to small 
changes in the effective quantum numbers. The sensitivity is usually determined 
by the degree of cancellation which occurs between positive and negative contri- 
butions to the radial integral (23). A further cause of sensitivity may also be 
noted. In transitions vl>e’]+1 with v small (for example, in certain 2p—> e’d 
transitions) the radial integral will be sensitive to P,, for large values of p. For such 
values of p the sensitivity of P,, to v will be determined mainly by the exponential 
factor e~*” in (28). To obtain a measure of the sensitivity we write equation 
(59) as ; el’ 

g(vl; ¢'l’)= oa ode) (66) 
wae. de )=aly tne’) +x(; <T)]. (67) 
We may then say that sensitivity will be great when cos ¢ is small. 

Radiative transition probabilities may be expressed in terms of the matrix 
elements of the dipole length, of the dipole velocity or of the dipole acceleration. 
All three give identical results when exact wave functions are employed but may 
give significantly different results when approximate wave functions are used. 
The dipole length expression gives the greatest weight to large radial distances and 
is therefore to be preferred in our method of calculation. In the following 
discussion we do not consider the relative merits of the alternative expressions 
when other methods of calculation are employed and when quoting results 
obtained by other methods we consider only the use of the dipole length formulae. 

In certain previous calculations accurate bound-state wave functions have 
been employed together with undistorted regular Coulomb functions for the 
continuum states. In such cases the accuracy of our calculations may be tested 
by putting »’ =o in our formulae. 

We consider separately transitions of the following types: (i) ms—e’p, (ii) 
np-—>«'s, (iii) mp—>e'd, (iv) nd-><'f, (v) nd—e'p and nf->«'d and (vi) initial n/* states 
with g>1. 

(i) ns—>e'p transitions 

He 21S. Both Huang (19) and Goldberg (20) have used g-parameter bound- 
state functions together with undistorted Coulomb functions for the ejected 
electron. Their results are in close agreement and, as shown in Fig. 3(a), are 
in good agreement with our results calculated with »’=0. For this case the 
sensitivity is small (cos¢(o)=0-82). In Fig. 3(@) we also give results calculated 
with values of y’ determined from extrapolated quantum defects of the n'P 
series. The correction for distortion is seen to be quite small. 
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He 2*S. The sensitivity is not large (cos¢(o)=0°67). Fig. 3(6) shows that our 
results with 4’ =o are in good agreement with results obtained by Huang (19) 
using a 9-parameter bound-state function together with undistorted Coulomb 
functions for the ejected electron. Goldberg (20) has used a less accurate 3-para- 
meter bound-state function and, for the ejected electron, has used regular Coulomb 
functions calculated with an effective value of z suggested by n*P energy level 
data. Although this procedure makes some allowance for distortion it may be 
criticized on the grounds that the continuum functions will have incorrect asymp- 
totic forms. In Fig. 3(6) Goldberg’s results are compared with our results 
calculated with «’ determined from n*P energy levels. It is seen that the effect 
of distortion is similar in the two cases and it is a good deal greater than that for 
24S. Dalgarno and Kingston (21) have shown that our results, with allowance 


(a) He 2 '$ (b) He 235 


, 


0 








Fic4 3+-(a) He 2'S 
- -— -—-— Calculations of Huang (19) with p’=0 for the «'P continuum. 
{ ‘—.—.-— General formula with p’ =o. 
General formula with yw’ obtained from the n’ 'P series. 
(6) He 2*S 
- --- Calculations of Huang (19) with p’ =0 for the «*P continuum, 
General formula with p’ =0. 
Calculations of Goldberg (ao) with some allowanze for distortion in the 


<®P continuum. 
General formula with x’ obtained from the nP series, 


for distortion, are consistent with the requirement of continuity between bound- 
bound and bound-free matrix elments, the bound—bound matrix elements being 
calculated by the Hartree-Fock method. The results of Huang are not con- 
sistent with this requirement. 

Li 2s. In view of surprisingly large differences between experimental and 
theoretical results we consider both the 2s-np series and the 2s->e’p ionization 
continuum. For the bound—bound transition a> the oscillator strength is 
defined by 


f(b, 2) =4(E, — E,) —8(6, a) (68) 


where the energies E,, E, are in Rydberg units and where the line-strength S is 
in atomic units. For ionizing transitions we define 


a) =42"(I, 3%) S(¢,2) (69) 
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where S(e’, a His defined by (5) and (6). In Fig. 4 we plot 4v’*f(np, 2s) against 
—1/v' and df(e’, 2s)/de’ against +/e’. All points should lie on a single smooth 
curve (as may be shown by the argument leading to equation (54)). We give 
the experimental results of Fillipov (22) for $v’*f and the experimental results 
of Tunstead (23) for df/de’. These independent experimental results are seen 
to be reasonably consistent one with the other. We also give the results of Hartree— 
Fock calculations for $v'*f (Fock and Petrashen, 24) and for df/de’ (Stewart, 
25). Finally we give the results of Bates and Damgaard for transitions to 3p, 
4p, and sp together with our results for transitions to higher excited states and 
to the continuum. 

Since the sensitivity is only moderate (cos (0) = 0°36) the difference between 
the experimental and the theoretical results are larger than would have been 
expected. 











01 


Fic. 4.—Oscillator strengths for Li 2s-»>n'p and 2s-> «'p. 


Experimental results of Tunstead (a3) for (df/de’) and curve through 
experimental points (©) of Fillipou (a2) for (vf y/2). 

Results of Hartree-Fock calculations made by Stewart (ag) for (df/de’) and 
by Fock and Petrashen (24) for (v'*f y/2). 

General formula. For (vf y/2) the points [-) are obtained from the Bates 
and Damgaard tables and the points x from our calculations. 


If the experimental results are correct there must be comparatively large 
errors in the separable wave functions used in the Hartree-Fock calculations and 
in our method of calculation. 

Although it would be of interest to make further calculations using non- 
separable variational functions for the Li ground state it is by no means evident 
that any substantial modification of the calculated cross section would be obtained. 
The Hartree-Fock term value, €,, = — 0-393, is in very good agreement with the 
experimental value, €,,= —0-396. The Hartree-Fock mp term values are also 
in good agreement with experiment. Our calculations, based on experimental 
energies, would not be in such close agreement with the Hartree-Fock calculations 
if there were substantial errors in the Hartree-Fock energies and phases. 

Na 3s. In Fig. 5 our results for the photo-ionization cross section are compared 
with results of Hartree-Fock calculations (Seaton, 26) and of an experimental 
determination (Ditchburn, Jutsum and Marr, 27). Since the sensitivity is 


10 
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considerable (cos¢(o)=0-12) and since the error in the experimental results 
may be as great as +20 per cent, the general agreement between theory and 
experiment may be considered to be satisfactory. ‘There is as yet no satisfactory 
explanation of the non-zero minimum in the measured cross section, the calculated 
cross sections having zero minima corresponding to a change in the sign of the 
matrix elements. 

In Fig. 5 we have also plotted some results obtained by Rudkjobing (28). 
These are discussed further in the next paragraph. 
Na 4s. Agreement with calculations by Rudkjobing (Fig. 6) must be regarded 
as poor even when allowance is made for fairly great sensitivity (cos ¢(0) = 0-20). 
Rudkjobing used a Hartree-Fock 4s function and calculated a continuum function 
Gy, using a central field adjusted to give approximate agreement with experi- 
mental energies for 3s, 3p and 3d. To these functions G,,, he then added linear 
combinations of Hartree-Fock radial functions P;,, and P,, in such a way as to 








0.10 


Fic. 5.—The photo-ionization cross section for Na 3s. 
Experimental results of Ditchburn and Marr (27). 
Hartree-Fock calculations (Seaton, 26). 

Values calculated by Rudkjebing (28). 
General formula. 


obtain «’p radial functions orthogonal to P;,, and P,,. This addition of combina- 
tions of P;,, and P,,, produces an undesirable distortion of the ep radial functions 
at fairly large radial distances. No such distortion appears in the Hartree-Fock 
radial functions (Seaton, 26). Since the dipole matrix elements are large for 
45-3p and particularly for 4s-4p the admixture of even small amounts of 3p and 
4p functions has a considerable influence on the calculated 4s—e’p matrix elements. 
Similar considerations apply to the 3s-<’p transitions and for this case the 
moderately good agreement between our results and those of Rudkjobing appears 
to be largely fortuitous. 

Mg?* 38 and Si+ 4s. Biermann and Liibeck (29) have calculated cross sections 
for ionization from several states of Mg* and of Sit. To the central potentials 
obtained from non-exchange self-consistent field calculations they add polariza- 
tion potentials behaving asymptotically as a/r* where the polarizability « is 
estimated from quantum defects for the f series. The potentials are further 
modified by including a function intended to make approximate allowance for 
exchange interactions. This function contains an adjustable parameter Af which 
is assumed to be constant for the bound and continuum states of each / series, 
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Fic. 6.—Na 4s. 


Values calculated by Rudjebing (28). 
General formula. 


the value of Af being chosen so as to give the correct energy for one of the bound- 
states of the series. This method should give good wave functions for the state 
used for the determination of AB but may give less accurate functions for the 
remaining bound-states and for the continuum. For the p-wave potentials the 
values of Af are determined from the 3p state of Mg+ and from the 4p state of 
Sit. In Fig. 7 our results are compared with those of Biermann and Liibeck. 


(a) Mgt 3s ” (b) Si* 4s 











Fic. 7. 


(a) Mgt 3s, (6) Sit 4s 
---- Biermann and Liibeck (a9) 
General formula 


The sensitivity is considerable (cos $(0)=0-19 for Mg* 3s and cos $(0) = 0°08 for 
Sit+4s) and the differences between the results obtained are therefore no greater 
than would be expected. 

K 4s. For this case the sensitivity is so great that accurate calculation of the cross 
section is scarcely possible. Bates (30) has shown that approximate agreement 
with the measured cross section (27) may be obtained using a Hartree-Fock 4s 
function and ¢’p continuum functions calculated with a suitably modified central 
field. It is of interest to see if the continuum functions which give approximate 
agreement with the experimental cross section have phases consistent with 


10* 





140 A. Burgess and M. }. Seaton Vol. 120 


extrapolated quantum defects. We have therefore made calculations using a 
Hartree-Fock 4s function together with continuum functions of the form (36). 
From np energy levels we obtain p'(e’)=p'(0o)—o-261e’ with p’(o)= 1-711. 
We have made calculations for various values of y’(o). It is seen (Fig. 8(a)) 
that the best agreement with the experimental cross section is obtained with a 
value of p’(o) not far different from the experimental value (y’(0)= 1-690 in 
place of 1-711). 

In the general formula we consider the effect of using the extrapolated 
experimental quantum defect, pw’ = 1-711 —0-261e’, and of varying the effective 
quantum number »v for 4s. With the experimental value, v= 1-770, we obtain 
cos¢(o)=—o-15 and with the Hartree-Fock value, v=1-852, we obtain 
cos ¢(o)= +0°09. Fig. 8(5) shows that the best agreement with the experimental 
cross section is obtained with values of v intermediate between the experimental 
and Hartree-Fock values. 


- ¥=1.852 
(0) =1.711 “4 


1-0 
= y=1.830 


a“ 
ZF = 1-690 
Z ¥=1.770 


1 0.0 — 9" 810 


~ 0.00 005 0.10 0.00 005 gf «(00 














Fic. 8.—Cross sections for photo-ionization of K 4s 
(a) The experimental cross section (— ——) (27) and cross-sections ( ) calculated using a 
Hartree—Fock 4s function and continuum functions (36) with two different values of p'(o); 
p'(o)= 1-711 is the experimental value. 
(6) The experimental cross section (— ——) and cross sections ( ) calculated using the general 
formula with th: extrapolated empirical quantum defect and vurious values of the effective quantum 
number v for 4s. The experimental value is v= 1-770 and the Hartree-Fock value is v= 1-852. 


Comparison of the curve »’(o)=1-711 of Fig. 8(a) with the curve v= 1-852 
of Fig. 8(b) shows, for a sensitive case, the differences in the calculated cross 
section due to differences between the Hartree-Fock bound-state function and 
the function (28) calculated with the Hartree-Fock effective quantum number. 
Cat 4s. Although this is also a sensitive case the sensitivity is less than for K 4s. 
From mp energies we obtain »’=1-435—0-407«’. With the experimental 4s 
effective quantum number, v= 2-141, we obtain cos ¢(o)= +0-094 and with the 
Hartree-Fock value, v = 2-196, we obtain cos ¢(o)= +0°255. The corresponding 
threshold cross sections are 0-4 x 10~! cm?® (experimental v) and 2-8 x 10~"® cm? 
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(Hartree-Fock v). A full Hartree-Fock calculation (16, 31) gives a threshold 
cross-section of 1:2 x 10~!¥cm*, 

The Cat+ 4s cross section is required for the study of the interstellar ionization 
equilibrium. For the iso-electronic case of K 4s, and also for Na 3s, the best 
estimate of the cross section is obtained with a value 7 of v intermediate between 
the experimental and Hartree-Fock values. We may expect that this will also 
be the case for Cat4s. We put 


v = cv(experimental) + (1 —c)v(Hartree—Fock). 
For K 4s we obtain the best results with c~o-39. For Ca* 4s we adapt c=0°39 
to obtain }=2-175. Fig. 9 shows the cross section calculated with this value of 
y and, for comparison, the cross section calculated by the Hartree-Fock method. 








Fic. 9.—Cross sections for photo-ionization of Ca* 4s. 


- — — — Hartree-Fock calculation (16, 31). 
General formula with effective quantum number v = 2-175 for 43 (see text). 


The full line curve is probably the most accurate estimate of the cross section. Note that the 
kinetic energy of the ejected electron, in units of 13°60 eV, is k*=4e’. 


(ii) mp—>e’s transitions 

Na 3p-e’s and 4p->e’s. For these two cases there is no great sensitivity 
(cos¢(o)~o-8). Fig. 10 shows our results to be in good agreement with the 
results of Rudkjobing (28). For the s-continuum Rudkjobing used wave functions 
more satisfactory than those employed for the p-continuum (see discussion of 
Na 4s above). 

Mg* 3p->e’s and Si* 3p->e’s and 4p->es’. For these cases sensitivity is slight, 
cos¢(o) being almost unity. Our results are in satisfactory agreement with 
those of Biermann and Liibeck (29) (Fig. 10). 

Ca+ 4p-e’s. The good agreement with calculations by Green (31) (Fig. 10) 
may be largely fortuitous since Green neglected exchange in the calculation of 
continuum wave functions. For this case the sensitivity is very slight, 
(cos (0) = 0-999). 

(iii) mp—>e'd transitions 

Na 3p->e'd and 4p->e'd. There is no great sensitivity (cos¢(o0)~0-7). Fig. 11 
shows our results to be in good agreement with those of Rudkjebing (28). 

Mg* 3p->e'd and Si*+ 3p->e'd and 4p->e'd._ In Fig. 11 our results are compared 
with those of Biermann artd Liibeck (29). 
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Mg* 3p+e's 














Fic. 10.—mnp-— ¢’s transitions. 


- -- = Calculations of Rudjebing (a8) for Na, of Biermann, and Liibeck (ag) for 
Mg* and Si* and of Green (31) for Ca*. 
—- General formula. 
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No 3p+é’d Na 4p +6’d 





Si* 3p>e'd 

















Fic. 11.—np— ¢d transitions. 
Calculations of Rudjebing (a8) for Na, of Biermann and Liibeck (a9) for Mgt 
and Si* and of Green (31) for Ca*. 


General formula. 
General formula with yp’ =o. 
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The agreement for Mg* 3p->«'d is especially satisfactory in view of the fact 
that cos 6(0)= 0-16 and cos ¢(e’) goes through zero as «’ increases. 

The agreement for Si+ 3p->e«'d and 4p->«’d is poor despite the fact that 
sensitivity is only moderate (cos ¢(0)~ 0-6in both cases). This may be explained 
as follows. For the d-continuum Biermann and Liibeck used a central 
potential adjusted so as to give a correct energy for Sit 3d (see discussion 
of Sit 4s above) which is a perturbed member of the nd series. Fig. 12 shows 
quantum defects and y-defects (see QDM, Section 3) for Sit nd. This should 
be compared with Fig. 1 of QDM which shows quantum-defects and 7-defects 
for unperturbed d-series. It would appear that Si* 3s*3d *D is perturbed upwards 
by Sit 3s 3p? *D. The effect of the perturbation will be to give values of ,’ in 
the d-continuum functions of Biermann and Liibeck which will be too small. 
In Fig. 11 we give results calculated by our method with p’=o. It is seen 
that the results of Biermann and Liibeck lie between our results with p’=o 
and with p’ #0. 














~ ¥2 


Fic. 12.— Quantum defects (©) and n-defects (x) for Sit n'd. This figure should be compared 
with Fig. 1 of QDM, which shows quantum defects and y-defects for unperturbed d-series. 


Ca* 4p-e'd. The poor agreement between our results and those of Green (31) 
(Fig. 11) is probably due to neglect of exchange in the calculation of the continuum 
functions used by Green. It should be noted that the sensitivity (cos $(0) = 0°72) 
is greater than for the Cat 4p—’s case above. 


(iv) nd->«'f transitions 

Mg* 3d->e'f and Sit 3d->e'f. For these two cases the sensitivity is not great 
(cos¢(o)~ 0-7), and Fig. 13 shows our results to be in good agreement with 
those of Biermann and Liibeck (29). 

(v) nd—>e'p and nf->e'd transitions 

Mg* 3d->e'p and 4f-e'd. For 3d->«'p the sensitivity is not great (cos $(0) = 0°52), 
but Biermann and Liibeck (29) give values for df/de’ which are ten times greater 
than those obtained from the general formula. It does not appear that the values 
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Mgt 3d0f Si* 3d 











Fic. 13. 


(a) Mg* 3d+«f, (b) Si* 3d+ef. 
-— — — — Biermann and Liibeck (29). 
General formula. 


given by Biermann and Liibeck can be correct, as it would require a value of 
cos ¢(o) ~ 1-4 in the general formula to reproduce their threshold result. 

For 4f->nd, an accurate comparison cannot be made as Biermann and Liibecks’ 
results are given to only one significant figure; however their results agree with 
ours to within this accuracy. 

Si+.3d—e'p and 4f-e'd. Again Biermann and Liibeck give results to only one 
significant figure. For 3d->e’p their results and ours agree to within this accuracy. 
It is interesting to note however that the Biermann and Liibeck results for 


4f-<‘d lie between those obtained from the general formula with ,»’(e’d)=o0 and 
with p’(e'd) obtained by extrapolation of the nd series (see discussion of 


Si+3p—e'd and 4p-e'd above). 


(vi) Initial nl¢ states with q>1 


In the derivation of our formulae we considered radial equations of the type 

(18) and later identified the energy parameter E with the energy required 
to remove an electron. Although it is not easy to give a theoretical justification 
for this procedure for the case of initial states consisting of groups of equivalent 
electrons, the work of Bates and Damgaard suggests that useful results may 
nevertheless be obtained. For such cases (@u/d@v) may be quite large and careful 
consideration must be given to the factor [=(1+0/dv) required for the nor- 
malization of the bound-state wave functions. 
He 1s?. In Fig. 14 we plot quantum defects y against effective quantum numbers 
v for He 1s? 1S, 1s2s 4S and 153s 'S. For v=o the extrapolated quantum defect 
must have an integer value (Section 3 of QDM) and from the trend of the curve 
in Fig. 14, it would appear that this value should be unity. From the experimental 
data it is not easy to obtain a precise value of u/dv for 1s*. From the normaliza- 
tion of Hartree wave functions (Wilson and Lindsay, 32) we obtain 


(Ou/dv) = —0-167 
for 1s*. This gives the dashed tangent of Fig. 14 which is seen to be consistent 
with the other data. We adopt (du/dv)= — 0-167 for 157. 
Huang (19) has calculated (df/de’) using a six-parameter function for 15 
and undistorted Coulomb functions for the ejected electron. The sensitivity 
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Fic. 14.—Quantum defects for He 1s? 1S and 1s ns 'S against effective quantum numbers v. 
The tangent — — — is obtained from the normalization factor of Hartree wave functions. Note that the 
extrapolated quantum defect equals 1 for v=o. 


is slight (cos $(0) ~ 0°86), and Fig. 15 shows that the results are in good agreement 
with our results calculated with ~’=o. Improved calculations have recently 
been made by Stewart and Wilkinson (33). For 1s* they use the 6-parameter 
function employed by Huang and for the ejected electron they use Hartree-Fock 
‘ 








Fic. 15.—He 15* values of df/de’. 


Calculated by Huang (19) with undistorted Coulomb functions for the ejected electron. 
—————. General formula with p’ =o. 


functions. Let A(e’) be the correction to be added to (df/de’) calculated with 
neglect of distortion. In Table 1X values of A obtained by Stewart and Wilkinson 
are compared with those obtained by our method. It is seen that the corrections 
are small in both cases and that the agreement for small e’ is better than that 
for larger ¢’. 


TaBLe IX 


d 
(=) A(e’) 
de p=0 


Genera! formula General formula Stewart and Wilkinson 
1016 —0°043 —0°032 
0652 —o’o18 +0°O10 
0°407 —0°006 +0°028 
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Lassetre and Silverman (34) deduce oscillator strengths from experimental 
data on the scattering of high energy (~ 500 eV) electrons. When the energy is 
great enough for the Born approximation to be employed the scattering data 
may be used to obtain the generalized oscillator strengths defined by 


f(b, a K) = Ee 2) | (6) Pe" [a)|* (70) 


where K=k, —k, and where k,, k, are the initial and final propagation vectors. 
In the limit of K->o, f(6, alo)=f(b, a), where f(6, a) is the optical oscillator strength. 
The absolute values obtained by Lassetre and Silverman are obtained by cali- 
brating against f(2'P, 1s*| K) calculated accurately over a wide range of K values. 
Fig. 16 gives a comparison of optical oscillator strengths for transitions from 








0 


Fic. 16.—He 15s*: values of 4v"fy and dfjde’. 


— - — — Curve through experimental points (©) of Silverman and Lassetre (34). 
General formula (y’ #0). 
x Calculations using tables of Bates and Damgaard (1) and (8n/8v) = —0-167 for 1s*. 


He 1s? (values of 4v’%f, and of df/de’). In addition to the results of Lassetre 
and Silverman we give points for 1s*>2 1P and 3 'P calculated from the Bates and 
Damgaard tables with (du/dv)=—0-167 for 1s? and our values of (df/de’) 
calculated with allowance for distortion. In the range o<«’ <2:0 our results 
are in good agreement with the experimental results. At higher energies 
(e’ =2-60 and 2-86) the experimental results show peaks in oscillator strengths 
which are due to doubly-excited states of auto-ionization type. Such effects, 
which in other systems may occur at lower energies, are not properly taken into 
account in our method of calculation. 

Mg 3s”. For this case a plot of the 3sns 1S series quantum defects against v 
indicates that ,, ->2 as v>o, and the graph enables a fairly reliable estimate of 
€=(1+ép/dv) to be made since p,, does not vary as rapidly with v as in the 
previous case of the He m'S series (Fig. 14). In Fig. 17 our results for the cross 
section are compared with the experimental results of Ditchburn and Marr 
(35). Also shown are values of the cross section calculated using the Hartree- 
Fock 3s? function given by Biermann and Trefftz (36) and the continuum function 
Gy,(r) (equation (36)), with »’(e'p) extrapolated from observed 3s np *P quantum 
defects. The agreement obtained is probably as good as may be expected in 
view of the moderate sensitivity (cos ¢(o)=0-356) and the possible error of 
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+:20 per cent in the observed threshold cross section. The relative accuracy of 
the general formula results does not appear to be greatly inferior to that obtained 
using a Hartree-Fock bound function. 





0.0 : a 


0.00 0.05 010 4+ O45 





Fic. 17.—Cross section for photo-ionization of Mg 3s* 


— — — — Experimental results of Ditchburn and Marr (35). 
General formula. 
Calculations with Hartree-Fock 3s* function. 
Ions with outer 2p? configurations. In QD M (Section 4) it was shown that, for 


the 2p radial functions of atomic oxygen ions, a good approximation to the factor 


{ is! given by 
‘ a (v—1)(v+2)_ (71) 


v(v+r) , 


TaBLe X 





pd 





&o 
cos (0) cos $(0) 
General 
formula 














2p* °P ‘ 2°11 
2p? 'S 1°90 
2p* 3p « 2°63 


Na*, 2p%'S 2°47 
O* 2p°4S ; 3°35 


Mg**2p* 'S ‘ 3°30 
Ne** 2p* *P 3°77 
O** 2p? *P 4°16 


Ne** 2p? 4S ‘ ; 4°16 
O* 2p *P 4°22 


Ne** 2p? *P ‘ « 4°25 
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Fic. 18.—2p—'d and 2p->e's transitions in atomic oxygen ions. All calculations with p’ =o. 


- - — — Hartree—Fock calculations. 
——_ General formula. 
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This factor is adopted for all ions with outer 2p* configurations. From previous 
calculations (Bates and Seaton, (7); Seaton (37)) it is known that 2p->e’d transi- 
tions are more important than 2p->«’s transitions and that effects arising from 
distortion of the «’d continuum are fairly small. In order to assess the accuracy 
of the general formula for transitions from 2p configurations we compare our 
results with those obtained using Hartree-Fock 2p radial functions. In all cases 
we consider undistorted continuum functions for both ¢’s and «’d. Table X 
gives the results of threshold calculations for eleven different ions and Fig 18 
gives energy variations for atomic oxygen ions. The calculations with Hartree- 
Fock wave functions were made using formulae given by Bates (2). It should be 
emphasized that Table X and Fig. 15 are intended only to give comparison with 
Hartree-Fock calculations and that improved results, particularly for 2p->e’s, 
would be obtained by allowing for distortion in the general formula. 

It is seen that the overall agreement is good, but that our results are less 
satisfactory for systems, such as neutral atomic oxygen, for which the effective 
quantum number for 29 is particularly small. We may conclude that satisfactory 
results should be obtained on using the general formula for photo-ionization 
from 2p configurations and on using the Bates and Damgaard tables, with 
inclusion of the factor (71), for transitions 2p’>2p%—'!n's and n’d. 

8. Conclusions.—From the comparisons with results derived by other 
methods it may be concluded that our method will give good results whenever 
there is no great sensitivity to details of the atomic wave functions. 

Our main interest is in the application of our calculations to the calculation 
of rdcombination spectra and of total radiative recombination rates. For these 
problems one requires sums over large numbers of recombination rates to indi- 
vidual states. Whenever the calculations for any given states are sensitive to the 
wave functions the calculated rates will be small and will contribute little to the 
sums. Since there is no reason to suspect systematic errors in our method we 
should obtain particularly good accuracy in results involving sums over many 
States. 
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POLARIZATION MEASURES OF COMET AREND-ROLAND (1956h) 
AND COMET MRKOS (19574) 


M. K. Vainu Bappu and S. D. Sinvhal 


(Received 1959 June 15) 


Summary 

Measurements of polarization of Comet Arend-Roland and of Comet 
Mrkos made photoelectrically, using interference filters and conventional B, 
V filters are given. The interference filters were chosen to isolate the Na 
emission at 5890 A and the continuum at 4800 A and 4300 A. The polariz- 
ation at 4800 A in both comets is similar yielding a mean value of 20-1 per 
cent computed for a phase angle of 90°. This indicates that the agencies 
responsible for the continuum were identical in both comets. The polariz- 
ation of Na emission calculated for a phase angle of 90° for Comet Mrkos was 
37°4 per cent in agreement with the theoretical value for Na resonance 
radiation. 





1. Introduction.—The polarization of the light of comets was first detected by 
Aragp early in the nineteenth century. Subsequently, Secchi, Wright and others 
have made visual observations on many comets and have confirmed and extended 
the‘earlier results. These indicated fluctuations in amount of polarization similar 
to those observed in brightness and colour. 

The first systematic efforts to utilize polarization data on comets towards an 
understanding of the mechanism responsible for the origin of the light of comets 
were made by Ohman (x). Using a polarigraph in combination with an objective 
prism and telescope, he examined Comet Cunningham (1940c) and Comet 
Paraskevopoulos (1941c). A monochromatic study of this kind furnished data 
separately on the amount of polarized light in the continuum and in the light of the 
emission bands. In Comet Cunningham the emission spectrum was very intense 
as compared to the continuum, and the polarization measured in the light of the 
Swan bands agreed with what one would expect on the assumption that a 
fluorescence mechanism is responsible for the origin of the bands. Comet 
Paraskevopoulos on the other hand had a very strong continuous spectrum with 
hardly any bright bands present. Polarization measures indicated the dominant 
role played by scattering of light by small particles in the latter comet. 

The recent appearance of the two bright comets Comet Arend-Roland (1956h) 
and Comet Mrkos (1957d) provided a good opportunity at many observatories in 
the northern hemisphere for detailed studies in cometary physics. At Naini Tal 
we have carried out measures of the polarization of the nuclei of the two comets, as 
well as measuring their magnitudes and colours on a few nights. These measures 
were made photoelectrically, using interference filters for isolating selected 
portions of the spectrum in addition to the conventional filters employed in present- 
day photoelectric photometry. A preliminary account of some of our polarization 
data regarding Comet Arend-Roland was given earlier (2). Since the publication 
of these results Blackwell and Willstrop (3) have given an account of a similar 
study of Comet Arend-Roland employing photographic methods and combinations 
of broad passband filters to give near-monochromatic measures. For Comet 
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Mrkos, Hoag (4) has measured the polarization of the coma, as well as of selected 
portions of the tail, in integrated light, using a 1 P21 photomultiplier with no filters. 

2. Observing technique.—The 10-inch Cooke refractor of focal length 148 inches 
was used in conjunction with a photoelectric photometer. The light passed 
successively through a diaphragm in the focal plane, filters placed two inches behind 
it and a Fabry lens before reaching an unrefrigerated selected 1 P21 photomultiplier 
tube. A microscope inserted behind the diaphragm, at the will of the observer, 
enabled the checking of proper centring on the object. The diaphragm consisted 
of an eccentrically located circular six-position disk containing six different circular 
apertures, each of which could be brought into position centred on the optical axis 
by rotation. Three consecutive positions of the disk presented three apertures of 
increasing size. Of these, the smallest opening of diameter 2-15 mm was used in 
the measurement of magnitudes and colours. This opening admitted all of the 
coma and an unavoidable portion of the tail. In all measures of magnitudes and 
colours the coma was centred carefully in the aperture. The remaining three 
positions of the disk contained three nearly identical circular apertures of diameter 
210mm. These openings are exactly sixty degrees apart and notches provided 
on the outside of the disk ensure the shifting of each opening rapidly and precisely 
into position. Across these three successive openings, a sheet of polaroid was 
cemented on to the disk. Rotation of the disk, so as to bring each of the three 
openings in turn on to the beam, acts as the equivalent of a single polaroid turned 
successively in position angle through sixty degrees. It is well known that three 
measures of polarized light, transmitted by an analyzer shifted through three 
positions sixty degrees apart, give uniquely the plane of vibration of the incident 
light as well as the percentage polarization, provided the orientation of the 
polaroid in the photometer is known. 

The output of the photomultiplier is fed to a D.C. amplifier similar to the one 
described by Valley and Wallman (5) but with certain modifications introduced by 
Whitford and Johnson. The amplified signal is conveyed either to a Brown 
recording potentiometer or to an Esterline-Angus recorder. The former was 
used for Comet Arend-Roland, while the latter was used for measures on Comet 
Mrkos. Calibrations of coarse and fine sensitivity ratios of the amplifier were 
carried out frequently in combination with either of the recorders. 

The polarization data obtained directly from the recorder tracings need correc- 
tion before yielding final values of planes of polarization and percentage polarization. 
Two factors necessitate the use of such corrections. These are the non-identical 
areas of the three apertures behind the polaroid and the varying sensitivity of the 
photometer to different planes of vibration of the incident light. The first factor 
can be made very small, so as to be negligible. ‘The second can be compensated 
for by using proper correction factors. These correction factors were obtained 
by measuring bright nearby stars that had no polarization. Table I gives the 
correction factors used for Comet Arend-Roland and also those used for Comet 
Mrkos, along with their probable errors. The different correction factors for the 
two comets became necessary since a fresh polaroid was used for measures of 
Comet Mrkos. The correction factors used for Comet Arend-Roland are the 
weighted means of 15 values obtained with six different filters. Those for Comet 
Mrkos were obtained from measures through two filters on two different nights. 

The orientation of the polaroid in the photometer was determined by a labora- 
tory arrangement with the aid of an artificial light source and a Nicol prism. We 
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TaBie I 
Correction factors used for Comets Arend-Roland (1956h) and Mrkos (1957d) 
Correction factors 
Comet Polaroid Polaroid Polaroid 
position position position 
I II Ill 
Arend-Roland 0'9714 10000 09418 
+ + + 


00028 00000 o'o109 


Mrkos 0°9479 10000 1'0084 
+ = + 
00039 ©°0000 0°0038 


believe that the error in transferring these measures to the telescope does not 
exceed one degree. 

During our observing run on Comet Mrkos we carried out test measures of two 
stars that have noticeable polarization and whose polarization has previously been 
determined independently by Hiltner (6), Smith (7) and Hall and Mikesell (8). 
We planned these measures to act as checks on the instrumental performance from 
night to night. It is unfortunate that we do not have measures on a large number 
of stars for which polarization data exist. As such we cannot make a distinct 
comparative study of the efficiency of the different methods used. The results 
obtained on these stars are givenin Table II. Those for Rho Ophiuchi represent a 
mean value of the results obtained on three nights. HD 154445 was observed on 
only one night. The data are too meagre for any study of systematic errors, but it 


TaBLe II 
Values of polarization of Rho Ophiuchi and HD 154445 obtained by different investigators 
HD 147888 (p Oph) HD 154445 
Investigators P (per cent) @ (degrees) P (per cent) 0 (degrees) 


Hiltner 30 53 3°4 9! 
Hall 3°5 51 3°0 85 
Smith 2°8 52 3°5 88 
Bappuand Sinvhal 4:2+0°68 54 4°5 95 


does seem that the percentage polarization values obtained by us are systematically 
higher than those obtained by others by about one per cent. The position angles 
of the plane of vibration agree well within the errors of measurement. The errors 
of observation signified by the probable errors are very reasonable. In general, 
for the cometary measures we estimate the probable error of a value of percentage 
polarization to be less than +1 per cent, and that of a value of plane of vibration to 
be within + 5 degrees. 

3. Polarization of Comet Arend-Roland.—For obtaining measures of polari- 
zation at selected wavelengths we have employed interference filters. These are 
Bausch and Lomb second order filters having a peak transmission of 35 per cent 
and width at half intensity of 80 A. 

For Comet Arend-Roland we measured through such filters centred on 
5890 A, 5000 A, 4800 A and 4300A respectively. In an earlier publication (2) we 
mentioned having used an interference filter centred on 4700A to isolate the 
C,(1, 0) emission of the Swan system. We had depended on the manufacturer’s 
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calibration data and assumed that the filter was actually centred at 4700 A. Sub- 
sequent calibration in our laboratory indicated that the filter had its peak trans- 
mission at 4800 A. Assuch, much to our dismay, we have not been able to measure 
the polarization of any of the Swan bands in both comets. The 4800 A filter, 
however, transmits the continuum, without any contamination being caused by 
neighbouring emission. Fig. 1 shows the curve of intensity distribution in the 
cometary spectrum as on 1957 April 30-1. The transmission curves of two of the 
interference filters used (viz., those centred on 4300A and 4800A) for Comet 
Arend-Roland, drawn on an arbitrary scale for ordinates, are also given. These 
latter curves are meant to demonstrate only the wavelength regions isolated by the 
filters. The 5000 A filter picks up a portion of the C, (0,0) band. Similarly the 
5890 A filter isolates mostly the continuum and any sodium emission, if present. 
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Fic. 1.—IJntensity distribution in the spectrum of Comet Arend-Roland on 1957 April 30.1. 
The two lower curves represent the transmission characteristics of two of the interference filters used. 
The ordinates of these latter curves are also on an arbitrary scale, 


In the case of Comet Arend-Roland our measures started only after the sodium 
emission nearly ceased to exist and hence the values obtained through this filter 
indicate essentially the polarization of the continuum. We had no filter readily 
available for the elimination of any third order contamination of the radiation 
transmitted by this filter and hence it is likely that a small fraction of CN emission 
at 3883A has leaked through. The values of polarization obtained with the 
5890 A filter compare well with those of 5000 A and 4300A. As such the contri- 
bution of 3883 A to the measures of 5890 A must be considered to be negligible. 
The values of per cent polarization and position angle of plane of vibration 
obtained through interference filters are given in Table III along with phase angles 
11* 
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calculated from the ephemeris of Candy (9). By May 9 the comet had gone beyond 
the range of our telescope—polaroid—interference filter combination. The 
visual magnitude of the comet at this time was 7-3. The normal limiting 
magnitude for measures of precision with interference filters on the 10-inch 
refractor is 9-0. The absorption of light by the polaroid thus caused a loss of over 
one and a half magnitudes. 


Taste III 


Polarization measures of Comet Arend-Roland (1956h) through interference filters 
Amount of polarization Position 
Date Position angle of plane of vibration, « Phase angle of 
1957 r : rs Bs angle intensity 
May 5890 A 5000 A 4800 A 4300 A equator, 8 
4°658 17°5% eee 17°9% ee 76° 24’ 37° 48" 67° 
103° 107° 





5641 17°4% k 20°4% 73° 47 40° o1" 76° 
116° 114° 


6-648 19°7% 17°3% : 42° 12’ 74° 
119° 119 


7°655 12°9% y 16°8% 44° 16’ 79° 
121° 122° 


8-660 20°5% 14°5% 
123° 135° 


TasLe IV 


Py values of Comet Arend-Roland (1956h) 
Po 
(mean value) 
Wavelength per cent 
5890 A 20°17 + 2°34 
5000 A 19°53 42°17 
4800 A 20°30 + 0°62 
4300 A 18-66 + 0°78 


TABLE V 


Polarization measures of Comet Arend-Roland (1956h) through B, V filters 
V filter B filter 

Date Amount of Position Amountof Position Position 
1957 polariza- angle of polariza- angle of Phase angle of 
May tion plane of tion plane of intensity 

(percent) vibration (percent) vibration equator 

a a 

4°658 189 104° an mae 
5°641 18°3 111° 19°6 112° 
6-648 19°9 116° 15°9 116° 
7°655 17°4 118° 16°7 117° 
8-660 14°5 119° 16°5 122° 
11°656 12°8 128° wie ois 
14°662 14'0 147° 15‘2 145° 
15°662 15°3 146° 16°5 145° 
22°647 9°6 146° 14°4 145° 
26°650 6-2 148° 10°7 152° 


37° 48" 
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Assuming that the relation between amount of polarization and phase angle is 
given by 


_ Pq sin? 0 
* 1+ Py cos? 8’ 


where P, and P,, are the values of percentage polarization at phase angles 6° and go° 
respectively, we have calculated the mean values of P,, for the comet through the 
different filters. These are given in Table IV. The above relationship is valid 
essentially for fluorescence phenomena, and perhaps it is incorrect to use it for any 
of the wavelengths transmitted. However, it can be used as a convenient 
extrapolation formula. To this end, Table IV gives the most likely values of 
polarization at a phase angle of go° for the different wavelengths. 

For Comet Arend-Roland we have measures of polarization obtained through 
conventional B, V filters. These are given in Table V. We have used these 
filters so as to be able to follow the comet over a large range of phase angle. The 
broad pass-band filters, while unfit for a determination of the role of fluorescence 
in the origin of the emission bands, serve well in studies of the general polarization 
characteristics of comets. Our measures of the polarization through B and V 
filters serve to find a correlation between the change in position angle of the 
intensity equator and the position angle of the plane of vibration. The plane of 
vibration is nearly perpendicular to the intensity equator after May 14. Until 
May 14 it is tilted to the normal to the intensity equator by an angle larger than the 
probable error. We do not believe that moonlight has interfered with the mea- 
sures because the sky measures made during the dark of the Moon period are similar 
to those made while the Moon was past first quarter. This indicates that the 
contribution by moonlight to the overall intensity was negligible. In particular, 
the observations of May 4-658 through interference and B, V filters show the 
deviation of the plane of vibration from the position angle of the normal to the 
intensity equator to be largest. The Moon was only four days old at the time and 
the sky deflections were less than 0-3 of a Brown recorder division when measured 
through interference filters and polaroid. This value just exceeds the errors of 
estimation in reading of the Brown recorder charts. As such we are confident that 
the large deviation between the position angle of the plane of vibration and the 
normal to the intensity equator observed on certain days is real and is not something 
caused by any residual effects of uncorrected sky polarization. On the basis of this 
inference regarding the plane of vibration of May 4°658, we may conclude that the 
deviations marked until May 14 do not arise from spurious effects due to moonlifht. 
Using the relation mentioned previously, the Py, values calculated through the 
B, V filters are the following : 


Py (B) = 20°71 per cent + 1°33, 
Py (V) = 20-17 per cent + 1°84. 


Within the limits of measurement these compare well with P,, values obtained for 
the continuum regions centred on 5890 A, 4800 A and 4300 A. 

4. Polarization of Comet Mrkos.—News of the discovery of Comet Mrkos 
reached us during the peak of the monsoon period of cloudy weather. However, 
on four nights in August we had clear skies which enabled us to get some measures 
on the comet when its heliocentric distance was still less than 0-7 astronomical 
units. We had to limit ourselves to fewer filters for measurement because an 
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observing runonthe comet on any of these four nights lasted only forty-five minutes, 
by virtue of the low declination of this comet. On all four nights we had observa- 
tions through filters that transmit the regions centred on 5890 A, 4800 A and 4300 A. 
The 5890 A interference filter had along with it a Chance OY2 filter, which elimin- 
ated any possible contamination from the former’s third order. The remaining 
two filters centred on 4800 A and 4300 A were the same interference filters used in 
the study of Comet Arend-Roland. 

The polarization results of Comet Mrkos are given in Table VI along with 
values of the phase angle of the comet and the position angle of the intensity 
equator. The per cent polarizations at 4800 A and 4300A are almost identical. 
The per cent polarizations observed at 5890 A on August 19, 20 and 22 are systema- 
tically higher than those measured at 4800A and 4300A. This is due to the 
presence of Naemission. Spectrophotometric tracings obtained by Dr William 
Liller of the University of Michigan indicate the presence of Na emission even on 
August 23°04. Judging from our polarization data, we presume that Na emission 
must have been fairly strong on August 19. 


Tasie VI 


Polarization measures of Comet Mrkos (1957d) through interference filters 
Amount of polarization Position 
Date Position angle of plane of Phase angle of Mean 
1957 vibration, « angle intensity a—B 
August 5890 A 4800 A 4300 A equator, B 
19°618 186% 14°7% 14°5% 65° 21’ 60° 49’ 88° 
145 149 152° 


> 


20°618 16°7% 14°6% 13°1% 64° 40° 63° 35° go 
148 155 158 


21°620 15°7% 156%, 16°5%, 63° 25’ 66° o5’ g2° 
154° 156 163° 


22°618 17°1% 16°4°%, 15°5% 62° o9 68° 24’ 96° 
161° 166 166° 


Tasie VII 
Py values of Comet Mrkos (1957d) 
Pr 

(mean value) 

Wavelength per cent 
5890 A 22°02 + 0°83 
4800 A 19°78 + 1°00 
4300 A 19°21 + 1°35 


We find for Comet Mrkos small-scale fluctuations from night to night in all 
wavelengths, similar to those observed in the case of Comet Arend-Roland. 
These are, however, superimposed on the overall pattern of variation of percentage 
polarization with phase angle. 

The values of position angle of plane of vibration, for all three wavelengths, 
agree remarkably well within the probable errors of their determination. These 
position angles are perpendicular to those of the intensity equator. For Comet 
Arend-Roland we had found the angle between the two position angles to differ 
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considerably from ninety degrees for those phase angles at which Comet Mrkos 
was observed. 

On one night we were able to carry out observations through a V filter. This 
was on August 22-618 and yielded a polarization of 15-6 per cent and a plane of 
vibration having a position angle of 162°. 

The calculated mean values of P,, for Comet Mrkos are given in Table VII. 
Those for 4800A and 4300A are in agreement with those obtained for Comet 
Arend-Roland. The value of P,, for 58g0A is higher than those for the other 
two wavelengths. This is caused by Na emission. A valuc of Py for thi 
waveiength has intrinsically very little significance of its own, except to indicate the 
role played by Na resonance radiation. This is so because the percentage polariza- 
tion at 58go A, between phase angles of ninety degrees and sixty degrees for Comet 
Mrkos, is the weighted mean of the polarization of the continuum and the Na 
radiation. The intensity of the Na radiation is very much a function of the 
heliocentric distance of the comet. 

5. Magnitudes and Colours.—The comparison stars used for the two comets 
are givenin Table VIII. The magnitudes and colours were derived by individual 
comparisons and by using a mean extinction coefficient for the reductions, 
Because of the low altitude at which the observations were made, the probable 
error for a magnitude determination lies around +0™-o4 and for a colour around 

+0™-02. The derived magnitudes and colours are given in Table IX. 


‘ 


Tasie VIII 
Details of comparison stars used for determining magnitudes and colours of 
Comets Arend-Roland (1956h) and Mrkos (19574) 
Comet Star R.A. Dec. V B-V_ HD spectral 
(1900.0) (1900.0) type 
h m ° af m m 
Arend-Roland a Cam 4441 +66 10 4°45* +0°09* Bo 
B Cam 4 545 +60 18 3°92* +0°89* Gop 
8 Aur 5 51°3 +54 17 360t +0'99t Ko 
y U Ma 11 486 +54 15 2°43t +o0of Ao 


Mrkos A Ser 15 41°6 +07 40 4°43t +06of Go 
p Boo 14 27°95 +30 49 3:57t +120T Ko 
* Derived photoelectrically. + Taken from reference (12), 


Taste IX 
Magnitudes and colours of Comets Arend-Roland (1956h) and Mrkos (1957d) 
Comet Date 1957 V B-V 


m m 
Arend-Roland May 7.699 0°84 
8.620 °"90 
11°645 6°55 og! 
14.643 7°02 096 
22.628 7°87 0°90 
26°627 8-27 "90 
Mean 0’90 


Mrkos . 16.591 6°50 "90 
18600 6°65 088 
Mean 0-89 
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_ 6. Discussion.—In this investigation we have tried to determine the following 
characteristics of polarization of comets :— 
(1) the percentage polarization of radiation in the continuum and its variation 
with phase ; 
(2) the orientation of the planes of polarization for the different wavelengths 
studied and their variation with phase; and 
(3) the percentage polarization of Na emission. 


For a study of the first characteristic we have utilized the radiation transmitted 
through the 4300 A, 4800 A and 5890 A filters. The 5890 A values of polarization 
can be used for the continuum only when the sodium emission is totally absent 
oris negligible. ‘The transmission through the 5000 A filter is partly contaminated 
by the Swan emission at 5165A and hence the value of percentage polarization using 
this filter represents a lower limit to the value of possible continuum per cent 
polarization. The 4300 A and 4800 A filters yield values of the polarization of the 
continuum, though weak Swan emission can contaminate the result obtained 
through the 4300 A filter when the emission spectrum of the comet is very intense. 

Thesecond characteristic has been studied through both interference and wide- 
band filters, as can be seen from the data presented. 

The value of the polarization of the radiation transmitted by the 4800 A filter 
is 20°30 per cent+0-62 for Comet Arend-Roland and 19°78 per cent + 1-00 for 
Comet Mrkos. These values are extrapolated Py, values. They are remarkably 
close and demonstrate well the conclusion that in both comets the continuous 
spectrum was similar. The Py. values of the continuum at 4300 A are 18-66 per 

;gent + 0-78 for Comet Arend-Roland and 19-21 per cent + 1-35 for Comet Mrkos. 
The Py value of radiation transmitted by the 5000 A filter for Comet Arend- 
‘Roland is 19°53 per cent + 2°17 and is in agreement with the 4300 A values. We 
thus derive a value of 19°1 per cent for the mean continuum polarization at 5000 A 
and 4300A for both comets when the phase angle is ninety degrees. A similar 
value for the mean polarization at 4800A is 20-1 per cent. The soooA and 
4300 A values would necessarily be lower limits to the real continuum polarization. 

Our measures of a B, V colour indicate a colour excess of o™-26, in the coma. 
This value would tend to increase when the contribution of Swan emission is 
evaluated. The colour excess certainly originates as a result of scattering by 
particles in the cometary atmosphere that are larger than the wavelength of light. 
There would, no doubt, be a mixture of various particle sizes, and of differentspecies 
of particles. It should be possible, theoretically, to derive the particle size and 
refractive index from the observed reddening and polarization. _Liller (10) finds, 
for instance, by such a method that the particles in the tail of Comet Arend-Roland 
were mostly of iron with diameters around 0-6. The B—V value obtained by 
Liller in the tail is + 1™-18, thus giving a colour excess of o™-54. It is, of course, 
possible to match the variation of polarization of the continuum with those observed 
by reflection from different surfaces in the laboratory, and get a qualitative idea of 
the nature of the reflecting surfaces inthecoma. Sucha procedure would perhaps 
give an incomplete picture of the scattering processes in the head of the comet. 

To calculate the percentage polarization of Na emission we have evaluated from 
Liller’s tracings of August 23-04, after applying necessary corrections for extinc- 
tion, filter factor and photomultiplier response, the amounts of continuum and 
various emissions given out by Comet Mrkos in the region 5740 A-6040A. We 
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find that the contributions of the continuum, 5845A emission and Na 
emission are 94°33 per cent, 1°83 per cent and 3-84 per cent respectively. After 
passing through the 5890 A filter employed by us, the composition of the light 
reaching the photomultiplier becomes, continuum 88-94 per cent, 5845 A emission 
2°44 per cent and 5890A emission 8-63 per cent. This is diagrammatically 
represented in fig. 2. Assuming the last mentioned composition to hold for our 
observations of August 22.618, we find on using the observed P, value with the 
5890 A filter for the total polarization of the light passing through that filter, and 
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Fic. 2.—Composition of light around 5900 A transmitted by the interference filter in the 
spectrum of Comet Mrkos on 1957 August 23.04. 


the P, value for 4800 A for the polarization of the continuum, that the P, value for 
uncontaminated Na emission is 27-0 per cent, which corresponds to a Py. value of 
37°4 percent. We have assumed in the above calculation the polarization of the 
5845 A emission to correspond to a P,, value of 10 per cent. Since the intensity 
of Na emission relative to the continuum is small, the value of polarization derived 
for the emission is very sensitive to changes in the overall observed polarization 
of radiation transmitted by the 5890 A filter. 

The sensitiveness of the Py) value of pure sodium emission to slight changes in 
the observed values of polarization, and the fact that a difference of about 01-4 
exists between our observations and Liller’s, strongly favour our conclusion that 
the estimated polarization of pure sodium emission is in close agreement with the 
theoretical value of 40 per cent for sodium resonance radiation. 

Uttar Pradesh State Observatory, 

Naini Tal, India: 

1959 June 6. 





162 Polarization measures of comets 1956h and 1957d Vol. 120 


References 
(x) Y. Ohman, Stockholm Obs. Ann., 13, 11, 1941. 
(2) M. K. V. Bappu and S. D. Sinvhal, Nature, 180, 1410, 1957. 
" (3) D. E. Blackwell and R. V. Willstrop, M.N., 117, 590, 1957. 
(4) A. A. Hoag, P.A.S.P., 70, 203, 1958. 
(5) G. E. Valley and H. Wallman, Vacuum Tube Amplifiers, (New York : McGraw-Hill 
Book Co. Inc.), p. 480, 1948. 
(6) W. A. Hiltner, Ap. F., 114, 241, 1951. 
(7) E. v. P. Smith, Ap. 7., 124, 43, 1956. 
(8) J. S. Hall and A. H. Mikesell, Pub. U.S. Naval Obs., 17(1), 1950. 
(9) M. P. Candy, I.A.U, Circ. No. 1585 (1957, Feb. 20). 
(10) W. Liller, A. 7., 62, 245, 1957. 
(11) H. L. Johnson and W. W. Morgan, Ap. 7., 187, 313, 1953. 





THREE-COLOUR PHOTOMETRY IN THE SOUTHERN COALSACK 


A. W. Rodgers 


(Received 1959 June 30) 


Summary 


Photometric data in the UBV system of 1060 stars in seven uniformly 
obscured regions of the Southern Coalsack have been obtained. These data 
are analysed by the method proposed by Becker. The photographic 
measures are based on a UBV photoelectric sequence extending to V = 14:4. 
Absorption is determined as a function of distance in each of the regions. 
The distance of the Coalsack is found to be 174 + 18 parsecs, its photographic 
absorption ranging from o°7 to 2°4 magnitudes, and the minimum mass of 
the absorbing material is estimated to be 14 so!:r masses. Behind the 
nebula lies a clear region which extends to 800 parsecs. At greater distances 
obscuration again sets in, averaging o-85 mag./kpc in the region 1 kpc to 
6 kpe from the Sun. 





1. Introduction —The Southern Coalsack is an isolated, irregular cloud of 
obscuring matter; its geometrical centre has galactic coordinates /= 270°, b=o°. 
Previous studies of the nebula have been made by Unsdéld (1), Miiller (2) and 
Lindsay (3), who have analysed star counts in the nebula and in neighbouring 


comparison regions. Both Unséld and Miiller found the Coalsack to have an 
average photographic absorption of one magnitude and a distance of approximately 
150 parsecs. Lindsay made counts in the area near the star HD 311999, an exciting 
star of the Hi region, No. 46 in Gum’s survey of Southern H11 regions (4). 
Lindsay obtained a distance of 160 parsecs and a total photographic absorption 
of 1-6 magnitudes for this region of the dust cloud. A spectrophotometric study 
of stars in the Coalsack and in the neighbouring SA 194 was made by Briick (§) 
in 1934. His most important conclusion was that colour excesses found in and 
out of the Coalsack showed no systematic difference. This suggestion of neutral 
absorption was one of the primary incentives for the present study. More recently, 
Houck (6) has classified slit spectra and measured photoelectric colours of 
members of the association I Cru which lies in part behind the Coalsack and 
at a distance of 2-2 kpc. From a comparison of colour excesses in the obscured 
and unobscured parts of the association, Houck derived a mean visual absorption 
of 1-6 magnitudes to the eastern half of the nebula; this value, however, depends 
directly on the assumed space absorption between the Coalsack and I Cru. 

In 1938 Becker (7) proposed a two-colour photometric system with base 
lines chosen so that the difference between the colours was independent of 
interstellar reddening. Let us denote the three effective wave-lengths of a 
photometric system by A,, A, and A,, with A,<A,<A,. The observed colours 
on this system are then 


Ca= Tan t+ Rb apfA), and Ci. = Tne + Rby (A), 


where & is constant, 5(A) describes the wave-length dependence of the selective 
absorption, and /,,, and /,,. are intrinsic colours on the A,,, A,, A, system. In the 
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general case, where the ratio of colour excesses is m such that n= E,./E,,, we see 
that the quantity C,,—nC,, is independent of reddening. Becker chose a 
photometric system where n(A) was unity and the colour difference provides a 
good method of discriminating between various spectral types. This is achieved 
by measuring in one colour the variation of the Balmer discontinuity with type. 
Becker and his colleagues (8, 9, 10) have since applied the three-colour RGU 
system to studies of reddening and membership of galactic clusters with 
considerable success, and Tifft (11) has recently established a system of 
photoelectric standards suitable for standardization of RGU photometry. 

In 1953 Johnson (12) proposed a three-colour photometry system which is 
well calibrated in terms of standard stars of the MK spectral—luminosity 
classification. It is widely used in modern photoelectric photometry and has the 
advantage that its equatorial standards are available to southern hemisphere 
observers. In addition, in Johnson’s UBV system the ratio of the colour excesses, 
Ey »/E,-y and the ratio of total visual absorption A, to the colour excess E,,_, 
have been well determined (13). Thus it is possible to construct a colour 
difference diagram, i.e. a plot of the colour difference (U — B)—(B—V) against 
(B— V), for the bright, nearby stars selected by Johnson as photometric standards. 
This diagram is shown in Fig. 1. For slight reddening, E,,_ ,/E, = 0°76, i.e. 


Ew p(B y/ Evy. y) = —~ 0°24. 


The arrowed line in Fig. 1 represents this equation. Becker’s original proposal 
was to apply the three-colour data to a determination of interstellar reddening 
as 4 function of distance in galactic fields by comparing the colour difference 
diagram for stars at a given apparent magnitude with that formed from nearby 
stars. Considering stars of the same apparent magnitude in a restricted field, 
we see that these extend over a considerable range in depth. For example, stars 
of luminosity class V form a spectral and luminosity sequence, the colour difference 
diagram of which, in an unreddened field, would be identical with that of Fig. 1. 
In the case of a field which shows linearly increasing absorption with distance, 
the spectral sequence for the dwarfs would be displaced along the reddening line 
by an amount dependent on the appropriate colour excess for the distance being 
considered. This distance is dependent, of course, on the apparent magnitude, 
absolute magnitude, and absorption. It is to be noted, however, that for the 
dwarfs the displacement along the reddening line is not constant over all intrinsic 
colours. Thus for stars of earlier type and higher luminosity the colour excesses 
will be greater than for the nearby red dwarfs of the same apparent magnitude. 
The result is that the colour difference diagram for luminosity class V stars will 
be both displaced and distorted; for the giant and supergiant sequences similar 
remarks apply. In addition, the high luminosity stars will be more reddened 
than the dwarfs of similar type and hence the giant and dwarf sequences will be 
separated in the diagram by an amount depending on the absorption distance 
gradient. 

In the case of a single absorbing screen, i.e. an isolated dust cloud, there 
would be a displacement but no distortion for stars behind the screen, and for 
stars in front there would be no displacement of any kind. Thus, depending on 
the distance of the dust cloud, we would see a discontinuity in the spectral sequence 
of the colour difference diagram at a point determined by the intrinsic colours 





1960 Three-colour photometry in the Southern Coalsack 





| 


REDDENING UNE 
t/t, -0O7% 


CLASS Y - MAIN SEQUENCE Ome 
CLASS  - GianTs —— 
CLASS lad -SUPERGIANTS ome 


g 





(U-8)-(8-Vv) 

















i 
O06 





(8-v) ——- 
Fic. 1.—The intrinsic colour difference diagram for stars on the UBV system, 


of stars at. the same distance as the cloud itself. For example, consider a dust 
cloud at a distance of 251 parsecs, i.e. distance modulus 7-0, which has a 
photographic absorption of one magnitude. The colour difference diagram of 
gth magnitude stars in the field will then differ from that shown in Fig. 1 in the 
following manner: 

The main sequence stars fainter than M,= + 2:0, i.e. types A3-As, lie in 
front of the cloud and are unreddened, and are hence undisplaced in the diagram. 
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Stars brighter than this limit, however, lie behind the cloud and they will be 
displaced in the diagram along the reddening line. The ‘component along the 
B—V axis of this displacement is the colour excess, which in this case is o™-25. 
Because we are dealing here with an absorbing sheet, we find displacement but 
no distortion. All stars of luminosity class III are brighter than M,= + 2:0, 
so that displacement will also occur for the giant sequence, again by o™-25 along 
the B—V axis. The position of the discontinuity in the sequence in the colour 
difference diagram corresponding to the unreddened A3-As V stars enables us to 
determine the distance of the absorbing sheet. The displacement in the diagram 
gives us the colour excess and hence the absorption. 

The observations.—Seven regions of the Coalsack were selected for study. 
The regions were chosen so that the absorption in each was as uniform as possible 
and that, in absorption, they typified the observed variations in the nebula. 
The seven regions in the field of the Coalsack are shown in Plate 4. The 
photometric plates were obtained with the Uppsala Southern Station Schmidt 
camera 20/26-inch aperture, f/3°5, and a Zeiss astrograph, 5-inch aperture 
(diaphragmed to 4 inches during the course of this work) with focal ratio f/6-5. 
The plate-filter combinations used were the following: 


Ultra-violet: 103a-O+2mm Chance OX1; Zeiss and Uppsala. 
Blue: 103a-O + 2mm Chance OY10; Zeiss. 

Blue: 103a-O +2mm Schott GG13; Uppsala. 

Red: 103a-E + 2°5 mm Red Plexiglass; Zeiss and Uppsala. 


The in-focus images were measured with a Haffner-type iris-photometer at 
Mt Stromlo Observatory. ‘Table I lists the region centres and the number of 


stars (N) measured in each. The photographic measures were calibrated by 
cbmparison with a photoelectric sequence in the UBV system, established in 
the centre of the nebula. The photoelectric observations were made with an 
RCA 1P-21 photomultiplier mounted at the Cassegrain focus of the 30-inch 


TaB.e | 


Region R.A. (1900) Dec. (1900) l b 
h m e:: ° ° 

I 12 51°5 —61 08 271'1 +08 

II 12 3671 —6o 48 269°6 +11 

Ii! 12 49°! —62 38 271°0 —O'5 

IV 12 34°6 —62 10 269°3 —o'2 

V 12 33°3 —63 20 269°3 —1°3 

VI 12 43°! —64 08 270°! —2°1 

VII 12 45°0 —61 18 270°6 +06 179 


Reynolds reflector at Mt Stromlo. The photometer and filters used are those 
described by Eggen, Gascoigne and Burr (14). Transfer to the UBV system 
was made through observation of Johnson—Morgan equatorial standards, and of 
Harvard E-region stars measured at the Cape Observatory for which transfers 
to the B, V system were available. Average values of the extinction found at 
Mt Stromlo were used; they were: 


k, =0™-18/unit atmosphere, 
k, = 0™-37 —0°04(B— V)/unit atmosphere, 
k, =0™-52—0°04(B—V)/unit atmosphere. 
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The region of the Southern Coalsack. The two bright stars are x and B Crucis; the bright 
cluster is the « Crucis cluster. The regions studied are marked by dashed lines. The print is enlarged 
from a 5" Zeiss astrograph plate, 30 minutes’ exposure on 103a-O emulsion. 


A. W. Rodgers, Three-colour photometry in the Southern Coalsack 
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The region of the photoelectric sequence. The print is enlarged from a 103a—O plate taken with 
the Uppsala Schmidt camera, exposure 15 minutes. 


A. W. Rodgers, Three-colour photometry in the Southern Coalsack 
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Tasie II 


Star V B-V U-B Sp Type (MK/HD) 
E6-A 6-06 +0°49 F8 

D 7°35 —oror Ao 

33 6°80 +117 Ko 

40 5°51 +1°51 K2 
7 Vir 4°26 +0°10 A3V 
61 Vir 4°74 +0°70 G6V 
B Lib 2°61 —O1l B8V 
B Oph 2°78 +1°16 Kalll 
y Oph 3°72 +0°04 AoV 
Cin 2322 7°54 + 1°36 K7V 
e Aqr 3°77 +0°02 A1V 
7» Hya 4°31 —0'20 B3V 
Cin 1633 8-49 +1°41 Mo. 5V 
HR 3314 3°89 —oO'o! , AoV 


The comparison stars are listed in Table II. The relations between the 
instrumental system and the UBV system were found to be linear and can be 
represented by the following equations: 


V = Vp—0°083C p+ constant, 
B-—V=0-96C p+ constant, 
U—B=1-09(B— U),+ constant. 


The final magnitudes and colours for the sequence stars are given in Table III. 
Column 1 gives the current number, 2 the HD or HDE number, 3, 4 and 5 the 


photometric data, 6 the Harvard spectral type and 7 the number of nights on 
which each star was observed. From the internal and night-to-night consistency 
of the results the estimated mean errors for V, B—V and U—B are +0016, 
+o™-o11 and +0™-o20 respectively. The sequence stars are marked on an enlarged 
print in Plate 5. 

While the effective wave-lengths of the photographic blue and ultra-violet 
systems correspond closely to those of Johnson’s B and U, the photographic 
red system used has an effective wave-length of 6400 A. Thus the photographic 
red—blue colours have a longer baseline than does the B—V colour. It was 
found that a colour equation of the type Vp= V—m(B—V), where Vz is a red 
magnitude and m isa constant, enabled a colour-free calibration of the photographic 
red magnitudes to be made. The best values of m were o-51 for the Zeiss 
astrograph plates, and o-60 for the Schmidt plates. From observations of 
standard sequences, e.g. the Pleiades and IC 4665, the photographic probable 
errors are estimated as +0™-o55 averaged over all plates, i.e., the errors in the 
U-—B colours are +0™-078 and those in the B—V colours +0™-078/(1+m), 
i.e. +0™-052. 

The photographic measures were made to the limit of the sequence at V = 14:4, 
and the only stars excepted from measurement in each region were those which 
composed obvious unresolved multiple images. Extinction and field error 
corrections were applied to the final magnitudes. The extinction corrections 
were applied to each region as a whole using average values of the coefficients 
as given above. The field corrections were determined empirically for each 
instrument using standard techniques. 





A. W. Rodgers 


Tasxe III 
The photoelectric sequence 


Zz 


HD V B-V U-B Sp. Type (HD) 
111464 6°62 + 1°39 + 1°37 Ko 
LI11174 8-08 +0°39 +0°25 A2 
111024 9°00 +o°19 — 0°25 B8 
111124 9°30 +0°64 —o'4I Bo 
111557 9°57 +0°52 +0°08 F 
111125 9°81 +o°10 —0°45 B8 
IIIOQI 10°38 +0°28 +0°16 Ao 

10°61 + 0°93 — 0°05 
312050 1105 +0°25 +0°02 Ao 
11°25 +1°61 +091 
311888 11°76 +0°59 +0°20 A: 
12°40 + 1°57 +1°49 
12°43 +0°42 +o'10 
12°64 +0°50 +0°22 
13°02 + 1°23 +1°08 
13°05 +0°73 +0°27 
13°10 + 0°63 +0°23 
13°44 + 1°09 +0°16 
14°38 +0°76 +0°59 


~ 
cu Anewne £ 


Warner DOAUwWE ASSP OOMWMO 4 


The sequence stars are marked on an enlarged print in Plate 5. 


t Analysis.—Because of the relative paucity of stars brighter than 14™-4 in the 
‘Coalsack, the limit of apparent brightness used in construction of the colour 
difference diagrams was taken as one magnitude. Thus, for each region, stars 
in the intervals V <9, 9-9°99, 10-10°99, ..., were separately plotted. Some 
examples of these diagrams for regions II, III, and VII are shown in Fig. 2. 
Where available, HDE spectral types are given so that the manner in which the 
various spectral types are delineated on the diagram is shown. Lines were 
drawn through the plotted points which were interpreted as corresponding to 
reddened spectral sequences of Morgan’s luminosity classes V and III. For 
o™-o5 intervals of the intrinsic colour (B— V),, the displacement along the B— V 
axis of the observed sequences, from the unreddened position shown in Fig. 1, 
were measured. The corresponding point on the observed sequence was projected 
back along the reddening line (the slope of which is assumed independent of 
intrinsic colour or luminosity) to obtain the intrinsic spectral type and, from the 
tables of Morgan and Keenan (15), the absolute magnitude. The relation between 
colour excess and total visual absorption has been determined by Johnson (13) 
and taken in this paper to be A,=3-0 E, _,; the curvature in this relation for 
large value of reddening is neglected. Thus for a given region, we can construct, 
for each interval of apparent magnitude, a plot of colour excess or absorption 
against distance modulus corrected for absorption. For successively fainter 
intervals of apparent magnitude, the absorption curves overlap in distance and 
hence provide a check on the consistency of the reductions. Secondly, two 
groups of stars, the A type dwarfs and yellow giants, are unambiguously placed 
in the colour difference diagram for all values of reddening, i.e., a close analogy 
can be drawn between multicolour integral photometry and Morgan’s concept 
of natural groups in the classification of low-dispersion spectra. From these 
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two groups alone, determinations of absorption against distance can be made 
to supply a further check on the sequences drawn in the observed colour difference 
diagrams. 

The resulting plots of absorption against distance for the seven regions are 
shown in Fig. 3. The curves display the following features: (a) a discontinuity 
due to the Coalsack at a mean distance 174 pc, where the absorption rises to 
half its value attributable to the nebula; and (b) a clear region extending to 
800 parsecs at which distance obscuration is again evident. Table IV lists the 
absorption and distance of the Coalsack as determined in the individual regions. 


Tasie IV 


Distance modulus Distance Photographic 
(magnitudes) (parsecs) absorption 

(magnitudes) 

151 22 

178 10 

138 Il 

191 o7 

174 2°4 

166 1's 

219 12 


It is estimated that differences in the distances of the regions as great as 
+16 parsecs may be real. The mean photographic absorption weighted 
according to the area of each region is 1-45 magnitudes. This value agrees well 
with those obtained from star counts (1), (2), and (3). This result, that the 
absorption derived from colour excesses and that derived directly from counts 
are in agreement, negates the suggestion of neutral absorption in Briick’s work (§). 

The determination of the mass of the obscuring component of the Coalsack 
rests on assumptions concerning the optical properties and size distribution of 
the dust particles. At a distance of 174 parsecs, the Coalsack has a projected 
area on the celestial sphere of 176 pc. From the Mie theory, Greenstein (16) 
has computed that ice particles with refractive index 4/3 and an effective radius 
of 3x 10-5 cm in a column 1 cm? in cross-section containing one gramme of 
matter will absorb 8-91 x 10 magnitudes in the photographic region. We obtain 
from Greenstein’s result a mass for the Coalsack of 2-7 x 10@ gm or 14 solar 
masses. ‘The mass of the dust component of the nebula derived under such 
assumptions will be a minimum value, and the actual mass may include 
contributions from particles outside this restricted size range. For all reasonable 
assumptions, however, the minimum mass falls between 8 and 30 solar masses. 
If the nebula is taken to be spherical, the density of dust is found to be 
4°4x 107%gm/cm*. This value is approximately 30 times the mean density of 
dust in the galactic plane as given by van de Hulst (17); however, from one 
part of the nebula to another the density may vary by as much as a factor of 
three. 

Behind the clear region which extends to 800 parsecs, a steady increase in 
absorption is noted. It is not suggested that the present work gives evidence 
for general absorption rather than that produced by isolated dust clouds; at 
the larger distances the resolution of the method, depending as it does on 
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photometric data, falls rapidly. In addition the colour excesses are derived from 
early type stars which have a large intrinsic dispersion in luminosity. The presence 
of the clear region behind the Coalsack together with the increase in obscuration 
up to the I Cru association at 2-2 kpc (6) suggest that at /= 270° we are looking 
through the Coalsack to an inter-arm region, behind which lies an extension of 
Morgan, Whitford and Code’s (18) Sagittarius arm. The star-poor galactic 
cluster NGC 4609 in region IV has been studied by the author. Magnitudes 
and colours were measured for approximately 40 stars within a radius of 
6 minutes of arc of the centre of the cluster. Photoelectric measurements based 
on the Coalsack sequence were made on 6 of these stars as zero point checks 
on the photographic measures. From the colour difference diagrams it was found 
that only 14 stars brighter than V = 14-4 were probable members of the cluster. 
The derived visual absorption was 0-95 magnitudes. The distance of the cluster 
was determined by classifying the brightest probable members on objective 
prism plates taken with the Uppsala Schmidt telescope and assuming that these 
stars were of luminosity class V. The spectral types ranged from Bs to A3. 
The derived distance is 1520 pc. The absorption agrees well with the value 
shown in Fig. 4 (d) corresponding to this distance. 

In conclusion it should be stressed that the basic limitations to the application 
of three-colour photometry to absorption studies are set by the assumptions of 
uniform obscuration over the field and the restricted number of stars in a field 
in which this assumption is valid. 
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ON THE MAGNETIC FIELD OF THE GALAXY 
F. Hoyle and }. G. Ireland 


(Received 1959 May 27) 


Summary 


The view that the spiral arms of the Galaxy are tubes of force winding 
into, and out of, the galactic nucleus is considered. It is shown that the 
winding of such a tube in the differential rotational field of the Galaxy leads to 
a degree of wrapping of about the observed order in a time interval as short as 
400 million years. 

It appears likely that a further winding of the field must lead to the 
outbreak of instabilities in which loops of the magnetic field emerge out of the 
plane of the Galaxy. It is suggested that the magnetic field of the galactic halo 
is built from a multitude of such loops. 

The rate at which energy passes from the orbital motion of gas in the 
Galaxy, first into magnetic energy, and then into random dynamical turbulence 
is computed. The result is 10~*’ erg cm=* sec™, a value in close agreement 
with the estimated rate of dissipation of the turbulence of the interstellar gas. 





1. Introduction.—Although it is still only ten years since Fermi (1) proposed 
the existence of a magnetic field within the Galaxy, this concept is now widely 
accepted. The discovery of the polarization of the light of distant stars by Hall 
and Hiltner added strong support to Fermi’s suggestion. 

The theory of the polarization effect proposed by Davis and Greenstein (2) 
requires the magnetic field to be mainly orientated along the spiral arms of the 
Galaxy, an idea that has been used by Chandrasekhar and Fermi in a discussion 
of the radial equilibrium of a spiral arm (3). 

Neglecting external fields, the difference of gravitational potential between 
the central axis and the surface of an arm is ~mG, where m is the mass per unit 
length and G is the gravitational constant. The cross-section of the outer arms 
of the Galaxy (Orion arm or Perseus arm) may be taken as ~1°5 x 10" cm’, 
while the total average mass density (gas + stars) can reasonably be estimated as 
~8x 10-“ gm cm-, giving m= 1-2 x 10%gm, and Gm=8x10"%erggm—. This 
may be compared with the random dynamical energy of the gaseous material 
within the arm, which, for random velocities of the material of order 5 kmsec~, 
is ~1:25x10%erggm-". If now the radial extension of the arm were entirely 
maintained by a balance between the dynamical motions and gravitation, the 
dynamical energy would have to be about 0-5 Gm (the exact factor in the virial 
theorem depending on the degree of isotropy of the motions) and this does not 
seem to be so. In fact the gravitational energy appears to be about six times 
greater than the dynamical energy. Hence Chandrasekhar and Fermi were led 
to postulate the existence of a magnetic pressure opposing the gravitational 
contraction of the arm. To keep the arm in balance the quantity H?/8mp must 
be of order 5 x 10! erg gm~", where H is the magnetic intensity and p is now the 
gas density. With p now of order 3x 10-“gmcm~ for the gas alone, the 
required value of H is ~6 x 10~* gauss. 
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This argument can hardly be said to be very decisive, however, since it derives 
its force from quantities that are not determined with precision—the estimates 
for the gravitational potential and for the dynamical energy are not very certain, 
nor was the shearing effect of the main gravitational field of the Galaxy taken into 
account in the calculation. The argument does show, however, that the magnetic 
intensity cannot be much greater than ~ 10~° gauss, otherwise the arms would 
expand. Particularly would this be the case in the direction perpendicular to 
the galactic plane. In this connection it may be recalled that the arms are not 
of circular cross-section, but appear to have a diameter of only ~ 200 parsecs 
in the direction normal to the galactic plane, as compared with ~ 500 parsecs 
or more in the plane itself (4). 

Support for the importance of magnetic pressure within the arms comes from 

an independent consideration of Davis (§), who pointed out that the magnetic 
pressure must be adequate to control the general gas motions, otherwise the 
overall direction of the field along the arms (as required by the polarization of 
starlight) would be destroyed. It can, however, be questioned whether Davis 
was right in suggesting that the field must be strong enough to control even the 
motions of dense clouds (density ~ 10-** gmcm~*, velocities up to 20 km sec~*) ; 
for this would require a magnetic intensity ~ 10~* gauss—and at such an intensity 
the whole arm would expand. Hence it is likely that the field will be locally 
distorted by such clouds. Indeed, radio observations by Mills (6) suggest that 
the field must possess an appreciable mean square component perpendicular to 
the arms. 
{| Inaddition to the field within the spiral arms, a field appears to exist throughout 
ia large halo surrounding the Galaxy. Until very recently, the diameter of this 
halo was set at ~30kpc, but observations by Hanbury-Brown and Hazard (7) 
show that the corresponding halo for M31 has a diameter of order 1ookpc. The 
halo surrounding the Galaxy probably has dimensions of a similar order. Accept- 
ing the synchroton theory of the origin of the radio emission from the galactic 
halo, together with a reasonable estimate for the energy density of the highly 
relativistic electrons (0-01 eV to o-1eV percm®), it can readily be shown that 
the magnetic intensity in the halo cannot be much less than 10~* gauss. The 
existence of such a field adds plausibility to the view that the spiral arm field has 
intensity ~ 10~° gauss; for plainly the field intensity is likely to be greatest within 
the comparatively dense material near the galactic plane. 

2. The origin of the galactic magnetic field.—While the above remarks corres- 
pond to views that are widely held by astrophysicists, there are other questions 
concerning the galactic magnetic field on which there is no agreement at all, and 
indeed which have scarcely begun to receive serious discussion. The origin of 
the field, and the nature of its large scale orientation, are two such questions. It 
would probably be generally agreed that the lines of magnetic force wind inwards 
to the centre of the Galaxy along one set of spiral arms, and then, after traversing 
the central regions, wind outwards along a second set of arms. But how are the 
lines directed at the outermost part of the galactic disk, say at distance 15~z0 kpc 
from the centre? 

A possibility is that the lines pass from one set of arms into the halo, and that 
after traversing the halo in some complex way, they are directed back into the 
other set of spiralarms. This possibility is especially attractive to those physicists 
who believe that cosmic rays originate within the Galaxy (8, 9), since it provides 
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a closed magnetic bottle within which the cosmic rays can be confined (except 
possibly at the highest energies). There are serious objections to this point of 
view, however. ‘These objections will be reviewed briefly. 

In the first place there is empirical optical evidence that the outer arms of the 
galaxies do not join with each other at great distance from the centre. The nearby 
Andromeda nebula is an Sb galaxy of very similar structural form to the Galaxy. 
Although a great deal of detail can be observed in the Andromeda nebula there is 
no sign of any such joining up of the outer arms. Indeed it is a general rule 
(except in some barred structures) that, down to the telescopic limit, the spiral 
arms of the galaxies appear always to head outwards into intergalactic space. In 
extreme cases, bridges actually exist that connect the spiral arms of one galaxy 
with those of a quite different galaxy. 

To this we may add a severe difficulty concerning the origin of the field. The 
dissipation of a field of galactic dimensions is very slow, so slow that the field cannot 
have the form of a closed bottle unless the field had a purely local origin. Yet a 
purely local origin demands the existence of an enormous impressed E.M.F. 
It can readily be shown (10, for example) that the self-inductance of the Galaxy 
is so large that an E.M.F. of ~10'*eV would have to operate for 10!’ sec in 
order to build a galactic field of the required intensity. It is emphasised that this 
E.M.F. must be ‘ seen’ by the material—it is not one that simply transforms to 
zero in the rest frame of the material. No suggestion has ever been made that 
would provide for an E.M.F. remotely approaching this magnitude. 

These remarks are thrown into relief by recalling older ideas on the origin of 
the galactic field. At one time it was believed that the field was highly irregular 
in structure, and that the intensity had gradually been built up through the effect 
of turbulence operating on a small initial ‘seed field’. It was suggested that the 
seed field owed its origin to an E.M.F. arising from thermal pressure gradients. 
The magnitude of such an E.M.F. is ~ 1 eV, ludicrously small compared with the 
10'*eV required to produce a large scale ordered field. The notion of turbulence 
had of course to be abandoned once it was realised that the field possesses large 
scale order. 

This apparently insuperable theoretical difficulty suggests that the galactic 
magnetic field is probably primaeval, and that the lines of force are not closed 
within the Galaxy—a view that has the advantage of providing the Galaxy with a 
direction, namely the initial direction of the magnetic vector. Subsequent 
winding of this direction leads to spiral structure. 

Although most questions are simplified by this picture, there is one problem 
that requires special consideration. If the lines of force in the spiral arms go 
out into intergalactic space, what is the nature of the field within the halo? This 
important problem will be discussed at a later stage. 


3. Orders of magnitude.—Taking gravitational forces first: 


(i) At distance 8kpc from the galactic centre: 
Radial acceleration due to whole Galaxy ~ 2 x 10-* cmsec~*, 
Acceleration normal to galactic plane~ 1-5 x 10-*scmsec*, z<0-2, 
where 2 is the distance from the galactic plane measured in kpc. (For z>0°2, 
the acceleration increases more slowly, rising to ~ 6 x 10-*cmsec™* at z= 2). 
Acceleration due to a local spiral arm ~ 10~*cmsec~*. 
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" (ii) At distance 4 kpc from the galactic centre: 

Radial acceleration due to Galaxy ~ 4 x 10-*cmsec~*. 
Acceleration normal to the plane ~ 10~’z cmsec™*, z< 0-2. 
Acceleration due to a local spiral arm ~ 5 x 10-*cmsec~*. 

Magnetic forces act directly only on ionized gas. The frictional coupling 
between the ionized component and the neutral component of the interstellar 
gas may be sufficiently strong, however, for the magnetic forces to be taken as 
acting effectively on the whole of the gas. We consider first the ionized component 
alone. The effect of frictional coupling will then be discussed. 


(a) Effect of the curvature of the magnetic lines of force 

Curvature of the lines of force introduces an acceleration towards the galactic 
centre, H*/(47rpion), where r is the radius of curvature and pion is the density 
of the ionized gas. According to Westerhout (11) we may put 

Pion = 3 x 10-6 gm cm-* 
at a distance of 8 kpc from the centre, and pion = 1°5 x 10-* gmcm~ at a central 
distance of 4kpc. Taking H=10-gauss at 8kpc and H=2x 10~*gauss at 
4kpc, the respective accelerations are 10-*cmsec~* and 1-6 x 10-* cmsec~. 
(6) Effect of the gradient of H* in the plane of the Galaxy 

Assuming that the magnetic vector changes its orientation (from ‘inwards’ 
to ‘ outwards’) as we pass from one spiral arm to the next along a radial direction, 
the value of H? presumably falls to zero between the spiral arms. Hence the 
magnitude of the gradient of H? in the plane of the Galaxy ~2H?/Ar, 
whert Ar is the radial separation of adjacent arms. The resulting radial accelera- 
tion ip the plane of the Galaxy is H?/(47Arpion). For Ar=1 kpc, the accelerations 
are ~6x10-*cmsec™ at 4kpc, and about 8 x 10-*cmsec~ at 8 kpc. 

(c) Effect of the gradient of H® normal to the galactic plane 

The gradients normal to the plane may well be three of four times greater 
than they are in the plane, corresponding to a significant decline of H? in the first 
hundred parsecs of height above the plane. 

Evidently then, the magnetic forces acting directly on the ionized component 
of the gas are considerably greater even than the main radial component of the 
gravitational force. Plainly, the ionized material will not obey anything like the 
usual system of dynamics unless there is a strong coupling with the neutral gas. 
This question will now be investigated. 

Only the friction between protons and the neutral atoms need be considered. 
Write n for the number of neutral atoms per cm*, and F for the magnitude of 
the magnetic acceleration. The mean free path of a proton between successive 
collisions with neutral atoms ~10%n-!cm, and the time between successive 
collisions is 10%"~'y~! sec, where v is the thermal speed of the proton. Under 
an acceleration F, the proton acquires a drift speed 10%n-v—'F cmsec~ in this 
time interval (the drift speed being small compared to v). The time required 
for systematic drift of distance Dcm between the ionized and the neutral gas is 
hence given by 

10-5” vF-'D sec. (1) 


With n=1, v=10'cmsec™!, F-!=10’cm—'sec?, D=3x10"%cm the time of 
drift is thus of order 3 x 10’ sec, which is close to the age of the Galaxy. Since 
D=3 x 10” cm is small compared to the scale of the Galaxy it follows that the 
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ionized and the neutral components of the gas can be regarded as bound together 
provided n is not much less than unity, as indeed it is not. In this connection it 
may be noted that the expression (1) is linear in nm. Hence it is not necessary to 
consider the detailed irregularities that exist in the density of the neutral gas (pro- 
vided that D is taken as a distance large compared with the mean size and spacing of 
the irregularities—as indeed it was). Moreover F was set high and v low in 
arriving at the time of drift. We conclude therefore that the two components of 
the gas are stuck effectively together by their mutual friction, and on this basis 
more realistic values for the magnetic accelerations will now be given. 

We simply replace pion by the total gas density. Taking this to be 
~ 3x 10-*4 9m cm~ at the galactic plane, and within the spiral arms, the resulting 
values are: 


(i) from the gradient of H? between adjacent arms, 
~ 3x 10~*°cmsec™ at central distance 4 kpc, 
~8x 10~!’cm sec at central distance 8 kpc. 


(ii) from the gradient of H? normal to the galactic plane, 
~ 10~*cmsec~ at central distance 4 kpc, 
~2x10~*cmsec~ at central distance 8 kpc. 


These accelerations amount to only about 10 per cent of the main radial gravita- 
tional acceleration of the Galaxy. Hence the rotational dynamics of the gas is 
largely unaffected by the magnetic forces. The magnetic accelerations are, 
however, comparable with the gravitational effects of a local spiral arm. In fact, 
this is just a restatement of the conclusions of Chandrasekhar and Fermi. The 
magnetic accelerations are also closely comparable with the total gravitational 
accelerations in a direction normal to the galactic plane. 

Before continuing the main discussion it is worth computing the rate at which 
heat is supplied to the interstellar medium through friction between ionized and 
neutral gas. It was seen above that the mean drift velocity of the ionized com- 
ponent relative to the neutral component is ~104%n-!v-1F. The kinetic energy 
of a proton at this speed ~5 x 10%~*y-*F*M,, where M, is the mass of the 
proton. This kinetic energy passes into heat in a time ~10%n~'v~'sec (the 
time between successive collisions) so that the mean rate of conversion into heat 
is ~5 x 10'4n—1y-1F2M,. Writing xn for the number of ionized atoms per cm’, 
the rate of heat production per unit volume is ~5 x 10“xu-'F*M,. Putting 
*x=0'02, v=10°cmsec~, and F=10~’cmsec~, the heat production 


~ 10-“erg cm~* sec, 


This is less by a factor ~ 10* than the rate of heating of the gas due to the dissi- 
pation of random mass motions (cf. Section 10). 

4. The winding of the galactic field—-The considerations of the previous 
section carry certain remarkable implications which, although well-known, do not 
seem to have received much serious attention in the literature. The gas in the 
Galaxy meves about the centre in nearly circular orbits, the dependence of orbital 
velocity on central distance having been given by van de Hulst, Raimond, and 
van Woerden (12). Now as the gas moves, the magnetic lines of force must 
certainly move with the ionized component to a very high degree of approximation 
indeed. But since the ionized and the neutral components are effectively bound 
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together, this implies that the lines of force move with the gas as a whole. That 


is to say, the same particles that lie along a line of force at one moment of time 
also lie along a line of force at a later moment. 

Suppose now that the magnetic vector is everywhere parallel to a particular 
direction at an initial moment, and suppose that the lines of force all traverse 
the central bulge of the Galaxy, as in Fig. 1. (In accordance with the remarks of 
Section 2, the direction of the magnetic vector may be thought of as the direction 
of the primaeval field from which the field of the Galaxy is derived.) The 
subsequent behaviour of the lines of force can then be determined a simple way. 
Consider a chain of particles lying initially along any particular line of force. 
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By using the rotational velocity curve it is easy to determine the positions of the 
same particles at any subsequent epoch, and the line joining these positions must 
be a line of force. 

The procedure is illustrated in the following figures for the case of the line 
of force that lies initially along the exact diameter (i.e. that passes through the 
very centre of the Galaxy): Fig. 2 gives the rotation curve used in the calculations, 


100 MILLION YEARS 


200 MILLION YEARS 














Fic. 4. 


| 400 MILLION YEARS 
300 MILLION YEARS 











Fic. 5. Fie. 6. 


and Figs. 3, 4, 5 and 6 give the shapes of the line to a distance of ~ 15 kpc from 
the centre for time intervals of 100, 200, 300, and 400 million years respectively. 
The magnetic vector may be taken as ingoing for the solid part of the line and as 
outgoing for the dotted part (or of course vice versa). A comparison between 
Figs. 4 and 6 shows that at twice the time there are twice as many wrappings of 
the field. This result applies generally, the number of wrappings rising linearly 
with time. 
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It is most remarkable that the degree of wrapping shown in Fig. 6 corresponds 
quite closely to the degree of wrapping found by Mills (6) to be actually present 
in the Galaxy. Yet the age of the Galaxy is some 25 times greater than the 400 
million year interval used in the construction of Fig. 6. Indeed, if the computa- 
tion were carried through for the whole age of the Galaxy, the number of wrappings 
would be 25 times too numerous. This conclusion would appear to follow 
inescapably from the rotation curve of Fig. 2. It seems therefore that either the 
rotation curve of Fig. 2 is grossly wrong—which appears very unlikely—or that 
some important new process must be at work. In the following sections we shall 
consider what process this might possibly be. 

We conclude the present section with the remark that, although the initial 
magnetic configuration shown in Fig. 1 is no doubt somewhat idealised in form, a 
similar problem would still arise for a more realistic initial configuration—i.e. 
one in which the lines of force possessed some initial degree of curvature. 

5. An energy source in the galactic nucleus.—Reference to Fig. 6 shows that 
in the topological sense there are no wrappings at all, since they can simply be 
deformed away at the centre! But since a line of force is a locus of particles such 
a deformation process is impossible so long as the particles move in circular orbits 
in accordance with Fig. 2. To permit deformation the particles would all have 
to move into the centre—a possibility too remote to command serious 
attention. 

By a similar token, the wrappings cannot be lost at the outer edge of the Galaxy 
unless the whole interstellar gas is also lost at the outer edge. Such a resolution 
of the problem would therefore require the whole present-day interstellar gas 
to be lost in a time ~5 x 10° years. A gaseous medium could not therefore be 
maintained unless it were resupplied by an outflow of material from the nucleus. 
The gas would have to flow through a distance ~3 x 10cm in a time of 
+ 1°5 x 10sec and hence would have to maintain an average outward velocity 
~ 20 km per sec throughout its passage from the nucleus to the outer edge of the 
Galaxy. 

In support of such a picture, gas is in fact streaming out of the nucleus, and is 
doing so at an adequate speed (13). Whether the rate of flow of gas is sufficient 
is somewhat doubtful, however—present indications are that the amount of gas 
may be inadequate by a factor ~5 to resupply the whole interstellar medium in 
a time of order 5 x 10% years. And a serious difficulty is that the gas would have 
to pick up angular momentum in some fashion as it moved outwards. The 
main radial gravitational acceleration of the Galaxy must always be closely com- 
pensated by centrifugal acceleration and this evidently requires that there shall 
be a marked increase of angular momentum as the gas moves outwards. This 
requirement can in fact be seen directly from an inspection of the rotational 
velocity curve. 

It is particularly to be noticed that a spirally wound magnetic field, of the sort 
considered above, does not supply net angular momentum to the gas, although 
it can reapportion angular momentum from one part of the gas to another. 
Indeed in order that net angular momentum be derived through the magnetic 
field it would be necessary for the field configuration to be dictated by some agency 
other than the gas itself, the field serving as a coupling device between the gas 
and this other agency. In view of our present inability to identify any such agency 
we accordingly think it worthwhile to consider other possibilities for resolving 
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the wrapping problem. Before we go on to consider these possibilities, there are 
certain features of the field configuration near the galactic centre that are worth 
mentioning. 

The curve of Fig. 2 has no relevance at small distances from the centre, where 
the true curve may possess irregularities (depending on the precise distribution 
of stars within the nucleus). Such irregularities could produce a very much 
quicker winding of the field than the curve of Fig. 2 does. Partly because of 
strong winding and partly because the gas density near the centre is low (~0°3 
atoms per cm*) the magnetic accelerations could be very large, in which case a 
compression towards the centre and a removal of the magnetic wrappings by 
deformation could well take place. In this way the very centre of the Galaxy 
emerges as a possible singular point where the energy of wrapping of the field 
can be partially released. 

Certainly some strong energy source in the central regions of the galactic 
nucleus is required to explain the remarkable observations of van Woerden, 
Rougoor, and Oort (13), although the actual situation seems to be still more 
complicated than that described above. (Rather does it seem as if the rotation 
curve may permit counterwinding at certain distances from the centre, such as 
would occur if the curve possesses a maximum and minimum.) 

6. The spiral arms as a transient phenomenon.—It might be suggested that the 
whole interstellar gas is circulated through stars in a turn-over time of about 400 
million years, and that the spiral arms are entirely renewed in a period of the 
same order. But this suggestion cannot work unless the rotation curve is dras- 
tically changed at some epoch during each circulation period. The lines of force 
must unwrap themselves in order that a configuration similar to that of Fig. 1 
be reinstated every 400 million years or so, and this can only happen if the magnetic 
accelerations become larger than the main radial gravitational acceleration of the 
Galaxy. In this case, however, the lines of force must expand far out from the 
galactic plane, so that the magnetic field becomes effectively lost from the immediate 
neighbourhood of the plane. The lost lines of force could only be retrieved by 
a widespread condensation of gas from the halo on to the galactic plane—and this 
invokes a second process apparently quite independent of the circulation between 
gas and stars. 

In any case the balance that apparently exists between magnetic energy, the 
dynamical energy of the interstellar gas, the size of the spiral arms, and the thick- 
ness of the Galaxy, would seem to be so delicate as to make any discontinuous 
resolution of the problem appear most improbable. Moreover, if a discontinuous 
explanation were correct, it is a little surprising that no hint of the explanation 
is given by an empirical examination of the structural forms of extragalactic 
systems. 

7. The origin of the halo field.—It is of interest to inquire what would happen 
if the winding of the field¢were to continue beyond the stage shown in Fig. 6. 
As the field becomes more and more wrapped, the magnetic intensity rises. But 
long before the magnetic acceleration becomes comparable with the main radial 
gravitational acceleration of the Galaxy the spiral arms cease to be in balance, in 
the sense of Chandrasekhar and Fermi. The balance fails not only because of 
the increasing magnetic intensity, but also because the mass per unit length of 
the arms decreases with winding—since their total length becomes steadily 
greater with time. 





F. Hoyle and }. G. Ireland, On the Vol. 120 


Still more important, the thickness of the distribution of gas and magnetic 
field in the direction normal to the galactic plane must increase. Reference to 
‘Section 3 shows that the magnetic and gravitational accelerations in this direction 
are already closely comparable under present-day conditions. A further winding 
of the field would hence produce a considerable increase in the thickness of the 
whole distribution. In fact the thickness must increase roughly linearly with 
time, as the following argument shows. 

Suppose first that the thickness z of the gas were to remain constant. Then 
the magnetic intensity H wouid increase proportionately to the number of wrap- 
pings—i.e. proportionately to the time ¢. Hence the normal magnetic 
acceleration would increase proportionately to ¢* (the magnetic acceleration is 
~ H?/(8mpz), where p is the gas density). Now let z change. Then Hatz-! 
and pocz-!, Hence H?/(8mpz) changes as f#z-*. The proportionality zoct 
therefore gives a constant normal magnetic acceleration, and this is approximately 
the correct condition, because the balancing inward gravitational acceleration is 
not much changed by the winding process, since it arises in a considerable measure 
from stars rather than from the gas. 

Now we have seen that the present-day value of z is closely consistent with 
t~5x 10% years. Plainly then, the thickness of the distribution of gas and 
magnetic field will swell to about twice the present value in only a further time 
of ~5 x 10° years, unless some other process intervenes. To understand what 
this process may be, we note that any swelling of the gas could hardly be extremely 
uniform over the whole face of the Galaxy. There must be particular places 
where the gas density is abnormally small, and in these places the outward magnetic 
acceleration must be abnormally high. Clearly then, bubbles must form. The 
bubbles, moreover, are unstable, since the gaseous material of the bubble must 

énd to slide sideways and downwards towards the galactic plane—much as the 

material of a solar prominence is observed to fall down into the chromosphere. 
As the material thus falls away, the gas density in the bubble declines, and the 
consequent rise of the outward magnetic acceleration causes the bubble to swell 
still further. It appears likely that the instability can only end by a complete 
breach of a portion of the face of the Galaxy, a breach through which a loop of 
the magnetic field emerges into the halo—again much as loops of the solar magnetic 
field emerge from the photosphere into the chromosphere and corona. _It is an 
attractive possibility that the magnetic field of the halo originates in this way. 
On this basis the field of the halo is an extraordinary magnetic complex emerging 
from, and returning to, the plane of the Galaxy at a multitude of points. 

8. Ingoing and outgoing spiral arms.—Although the discussion has now reached 
its most interesting stage, it seems worth inserting the remarks of the present 
section, even though the flow of the main discussion is thereby interrupted. 

In drawing up Figs. 3-6 it was assumed that the ingoing and outgoing arms 
interlace each other in the plane of the Galaxy (the solid line and the dotted line 
in the figures). If the Chandrasekhar—Fermi condition for the arms should fail, 
adjacent arms must begin to press on each other. It then becomes a serious 
question to decide whether the arms would continue to interlace each other in 
this way. An alternative possibility is that the ingoing and outgoing windings 
of the field would form two separate flat spirals, parallel to each other and on 
opposite sides of the galactic plane. Instead of individual arms pressing on each 
other, the two complete spirals would then be pressed together. For stability, 
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a sheet of material would have to exist along the galactic plane, and the gas pressure 
within the material would have to be comparable with the magnetic pressure in 
the spirals—i.e. with H?/87. For H=10~*gauss, and a gas temperature of 
~ 100 °K, the density in the sheet would thus be of order 10? atoms per cm*. 
Although the 21cm hydrogen line data would seem to exclude this possibility, 
the point is worth making—particularly as the whole question of diffusion between 
the ionized and neutral components of the gas (cf. Section 3) would have to be 
reopened. In Section 3, significant diffusion distances were taken to be 
100—1,00opc. Here however, diffusion through only 10 pc could be important, 
since the sheet of material along the galactic plane might be very thin. 

9. The nature of the spiral arms.—The argument of Section 7 gives encourage- 
ment to the suggestion that the winding problem is to be resolved by the emergence 
of the magnetic field into the halo. The magnetic field in the disk is wound 
continuously in a manner consistent with the rotational velocity curve. But 
after a certain degree of winding the magnetic intensity in the plane rises sufficiently 
for bubbles to break out of the galactic plane. Loops of the field then emerge 
into the halo with the result that portions of turns, or even whole turns, are lost 
from the plane of the disk. In this way, the wrappings of the field in the 
plane can be prevented from increasing steadily in number with the passage of 
time. 

We are now forced, however, to abandon the simple picture of the spiral 
arms as tubes of magnetic force winding smoothly into (and out of) the nucleus 
of the Galaxy. There is now no continuity of the lines of force. Instead, the 
arms become a statistical property of a collection of spiral arcs. The question 
evidently arises as to whether this less regular picture is consistent with the 
observational data. 

The most complete evidence on the spiral structure of the Galaxy comes from 
21 cm hydrogen line observations. The reconstruction of the hydrogen distribu- 
tion given by Schmidt (14) presents just such a collection of arcs as the present 
considerations require. Indeed it is fair to say that Schmidt’s picture of the hydro- 
gen distribution can hardly be interpreted as anything else but a collection of arcs. 

It is also well-known that while the eye readily picks out apparent spiral arms 
on optical photographs of M31 and of other Sb galaxies, the arms cannot really 
be traced continuously. (The present remark is not intended to apply to Sc 
galaxies. The whole question of the winding of Sc galaxies would require a 
separate discussion, in which the rotation curve, the mass distribution, and the 
magnetic intensity were all known. Such complete information is not at present 
available for any Sc galaxy.) ‘Taking accord of the well-known tendency of the eye 
to connect arcs and dots into a continuous line, the interpretation of the arms as 
a statistical phenomenon does not appear to us to be in discordance with the optical 
data. In this connection it may be noted that a statistical interpretation does not 
require all the arcs included in a particular arm to be wholly ingoing or wholly out- 
going—the two senses of the magnetic field could be mixed together in the same 
arm 


10. The dynamical energy input into the interstellar gas.—We come now to a 
strong point in favour of the present scheme. The process described above leads 
to a dynamical energy input into the interstellar gas, This input turns out to 
be precisely of the amount required to compensate for collisional dissipation 
within the gas. 
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Oort has estimated (15) a time of ~ 2 x 10sec for the decay of the random 
motions of the interstellar clouds of the solar neighbourhood—the decay arising 
from collisions between the clouds. At a mean cloud velocity of 5 kmsec™ 
and a mean density of 1 atom cm~* (obtained by smoothing the clouds through 
space) the mean dissipation is ~ 10~* ergcm~*sec™, 

Magnetic energy is lost in the expansion of a bubble. As the magnetic 
pressure lifts gaseous material from the galactic plane the magnetic energy is 
first converted into gravitational energy. The gravitational energy is then 
recovered by conversion into dynamical energy as the material slides out of the 
bubble back towards the galactic plane. Thus the energy input into the inter- 
stellar gas is dynamic rather than thermal and it comes fro:n gas falling inwards 
to the galactic plane. 

Once a steady state is reached, the energy input becomes equal to the rate at 
which work is done in winding the magnetic field. This can readily be estimated. 
For H=10~* gauss in the solar neighbourhood the present magnetic energy 
density H?/8x is ~4x 10-"ergcm~*. We have seen that H would be doubled 
by the winding in a time of about 400 million years—i.e. the magnetic energy 
density would be increased to 1-6x10-“ergcm~* in 400 million years. 
This implies an energy increase of 1:2 x 10~'ergcm~* in 400 million years 
(1-2 x 10sec), which gives precisely 10~*’ ergcm~*sec™!. 

It is emphasized that this significant coincidence does not depend at all on 
theory—except so far as the estimate H = 10~° gauss depends ontheory. Granted 
thié! field, and granted the rotational velocity curve, energy must be fed into the 
magnetic field at a rate ~10-" ergcm~*sec-4. The close coincidence of this 
value with the dynamical dissipation rate strongly suggests that a process must 
exist whereby energy fed into the magnetic field ultimately appears as dynamical 
energy. The mechanism of bubble formation supplies just such a process. 

The orbital motion of the gas in the Galaxy represents the primary source of 
the energy—it is this motion that supplies the mechanical forces required to 
wind the magnetic field. 

It was suggested many years ago by Weizsicker that an energy transfer from 
orbital motion to turbulent motion might take place. This appears to be so, 
but the emphasis is different from that of Weizsicker who used a Reynolds number 
argument to support the idea of an immediate origin of the turbulent motions. 
The present process is highly indirect. Orbital energy is first converted to 
magnetic energy. The magnetic energy is then converted into gravitational 
potential energy and only with a fall in the gravitational field does the final 
dynamical energy appear. 

11. Remarks on star formation.—The argument of the previous section 
implies that the dynamical energy of the interstellar gas is not derived mainly 
from the heating effect of early type stars. This agrees with an estimate by Oort 
(x5), who found that the observed dynamical input from expanding heated masses 
of gas amounts to a spatial average of ~ 10-* ergcm~* sec. 

The total emission of radiation by stars beyond the Lyman limit corresponds 
to a spatial average of ~ 10-" erg cem~*sec~! (cf. 16, for example). Of this about 
10 per cent, ~ 10~* ergcm~* sec~!, appears as heat in the interstellar gas. (‘The 
bulk of the energy is taken up in the ionization of neutral hydrogen atoms). The 
crucial question therefore is whether the conversion of heat to mass motions can 
have an efficiency as high as 10 per cent. Some investigators have concluded 
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that the efficiency might be as high as this, while others have arrived at a contrary 
result. Oort’s empirical estimate suggests a 1 per cent efficiency. 

We wish simply to point out that, if the heating and the gross dynamics of the 
interstellar gas were really controlled by the stars themselves, a distinctly odd 
situation would arise. 

The physical process responsible for the condensation of gas is almost certainly 
one of cooling and compression—cooled gas compressed by hot surrounding gas. 
The heating and the dynamics of the interstellar medium are thus intimately 
connected with cloud formation, and therefore with star formation. Hence, if 
the heating and the dynamics of the interstellar medium were to arise from the 
stars alone, star formation would be dependent on existing stars—i.e. star 
formation would possess feed-back. 

Write x for the density of the recently formed early-type stars (the ones 
supposed responsible for the dynamics of the medium). Then for a feed-back 
process, x varies with time ¢ according to an equation of the form 


>=; = xXe. 


dt 
With «>o and constant, x increases exponentially, With «<o and constant, 
x decreases exponentially. It is likely, however, that « varies with x. For 
a(x) increasing slowly with x, and with a zero at x, say, star-formation increases 
very greatly if x > x9, but decreases to zei aif x<x». Exactly this form of behaviour 
seems to be operative in the Magellanic Clouds—the former case in the Large 
Cloud and the latter in the Small Cloud. 

But the acceptance of a similar equation and a similar functional form for 
a(x) leads to serious difficulty in the case of the Galaxy, since star-formation in 
the Galaxy seems to be a steady process, not a luxurious effloresence like the 
state of affairs in the Large Magellanic Cloud, nor a miserable trickle like the 


Small Cloud. An equation for dx/dt with a stable solution is clearly required, 
such as 


. =xa+B, (8 a positive constant), 


which possesses a stable solution if «(x) and B are such that xa(x) + B = 0 possesses 
a solution. 


The likely existence of an extra term in the equation for the galactic case 
implies the existence of some new physical process—which might well be inter- 
preted-as arising from the steady input of dynamical energy into the gas due to 
winding of the magnetic field. 


St John’s College, Emmanuel College, 
Cambridge : Cambridge. 
1959 May. 
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AN ATTEMPT TO DETECT THE GALACTIC MAGNETIC 
FIELD USING ZEEMAN SPLITTING OF THE HYDROGEN 
LINE 


Jj. A. Galt,* C. H. Slater and W. L. H. Shuter 
(Received 1959 July 1) 


Summary 


An attempt has been made to determine the strength of the galactic 
magnetic field by observing the inverse Zeeman effect on the 21 cm absorption 
line of neutral hydrogen. Preliminary measurements using the Cassiopeia A 
radio source have shown no detectable Zeeman effect. This indicates that 
the magnetic field component in the line of sight is less than 5 x 10~* oersted 
at the point in the Orion spiral arm where the absorption occurs. 





1. Introduction.—A general magnetic field can be postulated to explain inter- 
stellar polarization of starlight, the cosmic ray spectrum, and the stability of the 
spiral arm structure of the galaxy. According to Chandrasekhar and Fermi (1), 
a magnetic field of the order of 7 x 10~* oersted may be expected although Davis 
and Greenstein (2) suggest fields up to 10~* oersted. 

Bolton and Wild (3) have suggested that the galactic magnetic field may be 
measured by observing the inverse Zeeman effect in the hyperfine structure of the 


21cm absorption spectrum of strong radio sources, using the radio analogue of 
Babcock’s (4) method of measuring weak solar magnetic fields. The present 
paper reports an attempt to make this measurement. 

2. Theory.—In the presence of a weak magnetic field the 21 cm line is split 


into three components. ‘The frequency separation of the extreme components 
is given by 


eH 
2mmc (t) 


where H is the strength of the magnetic field. A field of 10-* oersted will thus 
produce a Zeeman splitting, Av, of only 28c/s. As this is much smaller than the 
observed line-widths, the splitting cannot be observed directly. When the 
magnetic field is in the direction of observation, however, the components are 
circularly polarized in opposite directions and the small displacements may be 
measured with a system sensitive to changes in the sense of circular polarization. 
The technique used in the present experiment was to scan in frequency through 
an absorption line using a narrow band receiver which switched rapidly between 
opposite senses of circular polarization. It can be shown that for splitting of a 
line of gaussian shape the maximum temperature difference, AT,,, is related to 
the frequency separation Av by 


Av= 


Av=o-7 re ; (Av<p) (2) 


where p and 7°, are the half-width and depth of the line respectively. (See Fig. 1.) 


* Now at the Dominion Radio Astrophysical Observatory, Penticton, British Columbia, Canada. 
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Maximum sensitivity will be obtained with narrow intense lines such as those 
observed in the absorption spectrum of Cassiopeia A (23N5A) by Hagen, Lilley 
and McClain (5) and by Muller (6). The absorption lines originate in indi- 
vidual neutral hydrogen clouds of the galaxy. The narrowest of these lines, which 
has a half width of about 18 kc/s and is associated with the Orion spiral arm, was 
studied in this experiment. 























v 


Fic. 1.—(A) Theoretical profile of absorption line. 
(B) Displaced polarized components of a Zeeman triplet (splitting 
exaggerated). 
(C) Zeeman effect pattern obtained by subtracting the components in (B). 
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3. Equipment.—The 250-ft steerable telescope at Jodrell Bank was used with 
a triple superheterodyne hydrogen line spectrometer. A simplified block diagram 
of the receiver is shown in Fig. 2. The aerial feed was a wave-guide of square 
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Fic. 2.— Simplified block diagram of the receiver. 


cross-section flanged at the open end. The feed and its support were sym- 
metrical about a diagonal plane in the wave-guide to minimize spurious polarization 
effects. Two mutually perpendicular coplanar probes extracted orthogonal 
components from the radiation in the guide. These probes were matched to 
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separate crystal mixers fed in phase from a common local oscillator. The two 
' 27 Mc/s intermediate frequency signals were amplified by cascade preamplifiers 
and passed down equal lengths of cable to the laboratory suspended beneath the 
paraboloid. Here an extra quarter wave-length of cable was inserted in either 
lead to introduce a phase shift of go° between the two signals. A phase inverter 
which switched the phase by 180° 75 times a second was inserted in one of the 
leads. The signals were then added. The system was therefore sensitive 
alternately to left and right circular polarization every 1/75 of a second. After 
further amplification and detection the output was monitored on three recorders 
displaying : 


(i) Total power in a 600kc/s bandwidth, averaged between the opposite 
senses of circular polarization. 


(ii) Power difference between signals in a 5-2kc/s bandwidth and a 600kc/s 
bandwidth for observing the line profile. This again was the average between 
the opposite senses of circular polarization. 

(iii) Power difference between the opposite senses of circular polarization in 
a 7°7 kc/s bandwidth. 

The overall noise factor for reception at the hydrogen line frequency was 
estimated as 11; the output time constant of the 7-7kc/s polarization channel 
was 10 seconds and the frequency scanning rate about 2kc/s per minute. 

4. Calibration.—Noise measurements using an argon discharge tube showed 
that the set noise was approximately 3000°K. The total power received from 
Cassiopeia A (either in the broad band channel or in a narrow band channel away 
from an absorption line) was found to be 20 per cent of set noise. Since there 
was no primary image rejection this implies an aerial temperature of about 300°K 
at the hydrogen line frequency. The depth of the absorption line T,, as observed 
on the profile recorder was taken to be 300°K since absorption is almost complete 
(6). The polarization recorder was calibrated by removing the extra quarter 
wavelength of cable and injecting a c.w. signal into the two amplifiers in phase. 
By observing the ratio of the deflections thus produced on the profile and polari- 
zation records it was possible to assign aerial temperatures to deflections on the 
polarization recorder. It should be noted that the absolute value of 7, need not 
be known accurately since only the ratio of AT, to T,, is required to estimate the 
Zeeman splitting. The sense of the polarization was checked by radiating a 
circularly polarized signal at 1420 Mc/s from a helix into the main lobe of the 
aerial beam. 

5. Sensitivity of the equipment.—The r.m.s. noise on a single polarization 
record was about 10°K. Substituting the parameters AT,,=10°K, 7.,=300°K 
and »=18kc/s into equation (2) gives the minimum Zeeman splitting which 
would be detectable in a single scan as + 420c/s. By integrating 16 records this 
limit is reduced to +105c/s corresponding to a magnetic field of 3-8 x 10-% 
oersted. 

In practice this limit was not attained because of a spurious deflection of 
about 30°K which appeared on all polarization records of Cassiopeia A; its shape 
corresponding to that of the absorption line. The origin of this spurious 
deflection was not definitely established, but it is thought to have been caused 
either by cross-coupling between the probes in the wave-guide or by a systematic 
difference in gain of the paraboloid to left- and right-handed circular polarization. 
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To detect a Zeeman effect in the presence of this spurious deflection it was 
necessary to subtract the line profile from the polarization records. Any mis- 
alignment of the line centres in this process would produce a spurious Zeeman 
effect which was estimated to be +70¢/s. 

Thus in the present experiment it was not possible, in 16 scans, to detect a 
Zeeman separation of less than about + 125 c/s corresponding to magnetic fields 
of less than 4°5 x 10-5 oersted. 


6. Experimental procedure.—Before each series of observations precautions 
were taken to ensure that the receiver was balanced. These include equalizing 
the local oscillator injection into each mixer and balancing the gain and phase in 
each channel so that no square wave component appeared at the input to the main 
intermediate frequency section of the receiver. The telescope was then set to 
follow Cassiopeia A, the total power record being used to ensure that the source 
remained in the aerial beam. The second oscillator was set to scan from about 
35 kc/s above the line centre to 35 kc/s below. Coincidence marks were put on 
all records at regular intervals. After every two scans on Cassiopeia A a control 
scan was made on a reference region at /=82-2°, b= —2-38°. This region was 
chosen from the survey by Muller and Westerhout (7) because the 21 cm emission 
spectrum is similar to that in the immediate vicinity of Cassiopeia A. 


7. Analysis of results —Observations were carried cut during January 1959, 
and covered a total observing time of 269 hours during which approximately 200 
frequency scans were made on Cassiopeia A and the reference region. The first 
series of scans disclosed the presence of the spurious deflection referred to in 
Section 5, and to remove the systematic errors introduced by this effect it was 
necessary to improve the technique of aligning the profile and polarization 
records. ‘This series was therefore rejected and a new series was made with an 
improved technique of alignment. These latter records, taken on January 28 
and 29, have been analysed. They consisted of two sets of eight scans on 
Cassiopeia A and eight control scans on the reference region. For each scan on 
Cassiopeia A the position of the line centre was found from the profile record and 
the corresponding polarization record was divided into equally spaced intervals of 
about 3 kc/s symmetrically placed about the line centre. The average deflection 
from an arbitrary base line was then measured for each interval and the sum of 
the deflections was taken for each set of eight scans. The spurious deflection was 
eliminated by subtracting the line profile (which had been normalized and 
corrected for the slightly different bandwidths in the two channels). The 
observations were then analysed for the presence of a Zeeman component by 
calculating the cross-correlation between the observed curves and the expected 
shape of the Zeeman pattern. This cross-correlation when suitably normalized 
is a measure of the magnitude of the Zeeman component. The control scans 
were similarly analysed omitting the correction for the absorption line profile. 
The results for each set of eight integrated scans expressed in terms of the apparent 


Tasie I 
Observed Av 
(with r.m.s. error) 
c/s 
Cassiopeia A 1959 January 28 — 24+ 60 
1959 January 29 +10+60 
Reference region +40+ 60 
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Zeeman splitting Av are shown in Table I. (A positive sign indicates a magnetic 
field component directed away from the observer). 

The results for Cassiopeia A, when combined, indicate that the Zeeman 
splitting was Av= —7+42¢c/s which is clearly not significant. Taking into 
account that there is a much larger possible error of +125 c/s as explained in 
Section 5, we conclude that the longitudinal component of the magnetic field is 
less than 5 x 10-5 oersted. 


8. Conclusions.—No significant Zeeman effect has been detected and it is 
probable that the longitudinal component of the magnetic field in the clouds of 
neutral hydrogen which produce absorption is less than 5 x 10-5 oersted. It 
should be noted that the line of sight in the direction of Cassiopeia A is inclined 
at an angle of about 45° to the Orion spiral arm so that, if the general magnetic 
field is aligned with the arm, then the corresponding upper limit to the field must 
be raised by a factor of 4/2 over that quoted. 
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